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ABSTRACT 

Methods  for  obtaining  ionization  rate  contours,  fractionation 
numbers,  atom-concentration  intensity  ratios,  and  particle  size  distribution 
from  the  Miller  fallout  model  are  given.  Included  are  complete  computer 
programs  for  estimating  sublimation  pressures,  ionization  rate  contours, 
and  particle  size  parameters. 

Graphical  integration  of  activity  from  fallout  over  the  watersheds 
serving  Houston,  Texas  and  New  York  City  is  described.  Maximum  levels  of 
selected  contaminants  at  the  water  intake  for  these  cities  under  most  adverse 
wind  conditions  were  calcuilated  and  reported  in  ^lc/ml  as  follows:  New  York 
City;  Sr-89,  6.8  X  lO"^;  Sr-90,  5-9  x  lO"^;  Ru-106,  7.1  x  lO"^;  I-I3I, 

7.6  X  10"^;  Cs-137,  X  10'^;  Ba-l40,  5.8  x  10 Houston;  Sr-89,  I.9  x  lO'^; 
Sr-90,  1.9  X  10’*^;  Ru-106,  2.5  X  10'3;  I-I3I,  2.5  x  10“^;  Cs-137,  1.0  x  lO’^; 
3a-l40,  1.9  X  10'^. 

It  is  concluded  that  the  contribution  of  induced  radioactivity 
to  the  contamination  of  water  supplies  wovild  not  be  significant. 

Many  commercially  available  instruments  of  the  survey  type  can 
be  used  to  detect  adequately  the  presence  of  substantial  concentrations 
of  radioisotopes  in  water  under  emergency  conditions.  For  more  accurate 
measurement  of  such  concentrations  down  to  safe  levels  more  sensitive 
equipment  or  methods  of  preconcentration  of  the  samples  are  needed.  Such 
methods  or  instruments  should  extend  the  range  of  sensitivity  by  a  factor 
of  the  order  of  100  and  should  be  readily  adaptable  for  field  use.  Outlines 
of  procedures  for  the  analysis  of  radioactivity  in  water  are  given.  Efficiency 
of  water  decontamination  by  various  methods  is  discussed. 

Body  burdens  of  individual  radioisotopes  may  be  detemlned  from 
a  consideration  of  the  rate  of  ingestion,  the  effective  decay  rate,  and  the 
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affinity  for  the  isotope  of  the  critical  organ.  It  is  "believed  that 
uptake  of  radioactive  isotopes  may  he  reduced  or  prevented  hy  selective 
blocking  of  critical  organs  with  stable  isotopes. 
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I .  lOT  RODUCTIOK 


The  purpose  of  this  study  was  to  evaluate  aad  sunmiarize  availatle 
information  on  the  prohlem  of  water  contamination  hy  radioactive  fallout 
in  the  event  of  nuclear  war.  Consideration  has  teen  given  to  the  current 
theories  of  the  formation  and  distribution  of  fallout  and  the  level  of 
fallout  that  might  result  from  a  possible  nuclear  attack* 

Of  the  varioiis  attack  models  that  have  teen  recently  descrited'^’^ 
the  one  prepared  hy  Technical  Operations,  Inc.'  has  been  selected  as 
most  appropriate  for  the  present  study,  liiis  model  appears  to  have  teen 
well  thought  out  and  carefully  developed,  giving  due  consideration  to  the 
relative  importance  of  various  military,  industrial,  governmental,  and 
power  resource  targets.  This  model  assumes  all  detonations  to  he  .surface 
b'orsts,  although  perhaps  a  more  real. Stic  situation  woul,d  include  many 
air  hursts,  which  would  he  more  effective  against  an  '.’jahardened  militaiy 
or  industrial  target.  As  the  type  of  hurst  has  a  considerable  influence 
on  the  fallout  produced,  it  must  he  carefully  evaluated  in  a  study  of 
water  contamination. 

r  4 ) 

A  thorough  analysis  of  the  fallout  model  of  Miller'  has  heen 
made.  Several  important  functions  derived  from  this  model  have  been 
utilized  in  this  study  and  are  discussed  in  detail  in  Section  II,  "3vmiiiar;y 
of  Analysis  of  the  Fallout  Model". 

Information  from  the  Miller  model  has  also  been  used  in  the 
study  of  possible  levels  of  contamination  by  selected  isotopes  in  the  water 
supply  systems  of  Houston,  Texas  and  New  York  City.  With  appropriate  con¬ 
sideration  to  the  assumptions  made,  these  results  may  be  applied  to  the 
watersheds  of  municipalities  other  than  those  shown  here. 
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T.iiis  ,r  i;  ^dy  it  conceded  with  the  erfect  of  indlvidos.!  iiaclear 
ar-d  does  nc^  con.-ider  the  effect;--  of  overlapping  patters. 

T.he  ccirip.i.ex.  protiem  &f  overlap  has  teen  treated  by  Technical  Cperaticns, 
ihc.'"  E-nd  tut  E.d'lltlorL5l  studies  on  this  subject  appear  desirable. 

■.'■r.-jnspc.-:''.  of  iallcvt  cy  soriace  wat-er  has  teen  studied,  Et/alua- 
r -.on  o1  the  effect.:  of  vater.thed  eharacteriarics  or.  the  concentration  of 
r-,d,io'iC'’ ivit./  in  .t  ..rr's.ue  vater  vould  he  .-aluahle. 

ft eontafliirjr..tlcn  proced'ures,  both  standard  and  emergency,  have 
teen  described  previously. '"'v See  Appendix  A)  These  procedures  are  e-'-alnated 
uccs  ."ilng  to  their  efficacies  for  the  removal  of  various  radioisotopes  and 
-also  iccorli-'ig  to  the  feaeibility  of  their  use  after  a  nuclear  attack. 

purpoces  of  this  study,  it  has  teen  assumed  that  vatervtorks  facilities 
have  p-'iffered  only  nunimal  damage,  and  t-hat  i>erFonnel  will  be  available 
for  operation. 

A  study  of  published  methods  of  ra.diochemical  analysis  and 
existing  instrument at  Ion  he.s  teen  made.  Due  to  the  question  of  availability 
o."  personnel  and  facilities  it  appears  that  this  Information  will  be  of  more 
value  in  recognizing  long-*erm  hazards  than  in  determining  potability  in 
immediate  post-attack  period. 

The  biclcgical  iptake  an,!  resultant  body  burdens  in  man  have  been 
/S-'ed  in  rei.atlon  to  the  levels  of  water  contamination  considered  likely, 
f  t;,.- i:  g  .d-.tu  u-ss  been  put  into  a  -workable-  form  so  that  the  body  burden 
c:'  ir:.  of  the  se..ected  biclogically  important  radioisotopes  may  be  determined 
a--*-  -.nz  .Tpcci'. :c  irat  alter  the  init-.Lation  of  Intake, 

Ihe  internai  radiation  hazard  resulting  from  nuclear  attack  is 
cl  iiHport.^.^e  -uhiv  to  those  persons  who  sur’/ive  the  catastrophe  for  a 
perio'l  of  several  year.-.  As  a  primary  annlhllator,  it  is  of  no  consequence. 
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as  the  external  dose  concomitant  with  a  lethal  internal  dose  is  more  deadly 
"by  several  orders  of  Magnitude.  In  generalj  radiation  from  sub-lethal 
concentrations  of  internally  dejxssited  radionuclides  acts  slowly,  and 
its  i-esults  do  not  become  apparent  for  years.  Therefore,  the  internal 
radiation  hazard  is  nob  among  the  forces  debilitating  the  popi^lation  and 
reducing  its  inanediate  ability  to  cope  with  a  radically  altered  environ¬ 
ment.  Bather  it  has  long-term  effects  which  c\ilminate  in  producing  diseases 
such  as  cancer,  shortening  of  the  life  span,  and  possibly  affecting  the 
germ  plasm  of  mankind  to  a  small  degree. 


II.  SUMMARI  CF  ANALYSIS  OF  TIG!  O.C.D.  FALLOUT  MOUEL 


A.  Basic  Assumptions  of  the  O.C.D.  Model 

During  the  initial  phase  of  this  project,  the  Postattach  Resecrch 
Division,  Office  of  Civil  Defense,  requested  that  all  studies  be  based  on  the 
fallout  model  as  developed  in  "Fallout  and  Radiological  Countermeasures’'^^^ 

A  thorough  and  quantitative  study  of  the  model  was  conducted  because  the 
evaluations  of  water  contamination  depend  directly  on  previously  established 
Information  concerning  the  distribution  and  properties  of  local  fallout.  The 
results  of  the  detailed  study  are  now  presented  in  summary  form. 

The  model  was  based  on  both  established  physico-ohemicr.l  theories 
and  actual  test  data,  brought  together  by  the  scaling  method.  The  rationality 
and  simplicity  of  the  basic  assumptions  provide  this  model  with  many  advan¬ 
tages.  The  basic  assumptions  include: 

(1)  The  model  applies  generally  to  a  "ground  surface’’  burst,  but 
it  can  be  extended  to  other  tj’pes  of  detonations  by  suitable  transformations, 

(2)  The  required  input  data  consist  of  the  weapon  yield  and  the 
wind  velocity. 

(3)  The  radioactive  decay  schemes,  chemical  properties  and  rela¬ 
tive  abundances  of  the  fission  products  are  based  on  the  data  developed  by 
Bolles  &  Ballou'' Katcoff^^^,  and  Miller  and  Loeb^^\ 

(4-)  The  increase  of  the  fireball  volume  with  time  is  based  on 
the  adiabatic  expansion  of  an  ideal  gas  with  the  thermodynamic  equations 
modified  to  include  a  teim  for  the  change  in  free  eneigy  with  altitude  and 
with  the  external  pressure  proportional  to  exp  (-mgz/kt)  . 


>- 

For  explanation  of  symbols,  see  Glossary. 
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(5)  The  radioactive  cloud  is  assiamed  to  have  the  shape  of  an 
ob].ate  spheroid;  hence  the  ionization  rate  contours  or  fallout  patterns 
under  a  constant  wind  velocity  are  generally  cigar-shaped. 

(6)  All  values  for  the  dimensions  and  properties  of  the  fireball 
and  fallout  patterns  are  scaled  from  test  data  by  the  scaling  method  to  be 
a  function  of  weapon  yield  and  wind  velocity. 

(7)  The  condensation  process  of  fallout  is  divided  into  tvro 
time  periods:  (a)  the  first  period  is  characterized  by  the  presence  of 
gas  and  liquid  phases  and  ends  when  the  bulk  carrier  or  substrate  material 
of  the  particles  solidifies,  and  (b)  the  second  period  is  characterized 
hy  the  existence  of  gas  and  solid  phases.  The  temperature  of  i673°K  is 
used  to  divide  these  two  periods  of  condensation.  The  radioactive  material 
which  condenses  in  the  first  period  will  be  fused  inside  the  fallout  parti¬ 
cles  and  generally  becomes  insoluble  in  water.  The  material  which  condenses 
during  the  second  period  will  be  adsorbed  on  the  outside  of  the  fallout 
particle  and  becomes  readily  soluble  in  water.  With  these  basic  assump¬ 
tions,  the  model  will  yield  the  intensity,  the  activity  anu  atom-concentra¬ 
tion  for  a  specific  damwind  area  and  other,  related  information. 

B.  Formation  and  Geometry  of  Stem  and  Cloud 

After  a  surface  burst,  the  fireball  forms  the  shape  of  a  sphere, 
and  its  radii,  both  horizontal  and  vertical,  expand  exponentially  with 
altitude  as  it  rises.  Since  the  expansion  coefficients  are  different  for 
the  two  radii,  the  fireball  grows  to  the  form  of  an  oblate  spheroid,  with 
a  circular  top-view  and  an  elliptical  side-view.  When  the  fireball  reaches 
its  final  stabilized  height,  it  is  commonly  known  as  the  "cloud".  A  stem 
is  shaped  by  the  trajector;/-  of  a  continuously  expanding  fireball  and 
assumes  the  form  of  an  inverted  exponential  horn. 
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The  initial  spherical  radius  of  the  flrehall^  the  final 
horizontal  semi-axis,  a,  the  final  vertical  semi-axis,  b,  and  the  final 
height  of  the  center,  h,  of  the  cloud  are  related  to  the  weapon  yield, 
W,  through  empirical  data  as  shown  in  the  following  scaling  functions; 


Rg  =  2.09  X  10^  ft.  W  =  1  to  10^  KT  (l) 

a  =  2.45  X  10^  ft.  W  =  1  to  10^  ET  (2) 

b  =  1.40  X  xO^  ft.  W  =  1  to  10^  ET  (3) 

h  =  0.66  X  10^  ft.  W  =  1  to  28  ET  (^) 

h  =  1.68  X  10^  ft.  W  =28  to  10^  KT  (^a) 


At  a  given  altitude,  Z,  the  horizontal  and  vertical  serai -axis 
of  the  fireball  can  be  obtained  by 

k  Z 
a 

a„  =»  a  e 
Z  o 


(5) 


(6) 


where  a  ,  b  ,  k  ,  h,  can  be  solved  from  the  known  information  of  R  ,  a,  b,  h. 
o  o’  a'  tj  s’’ 

Ttie  geometry"  of  both  stem  and  cloud  is  shown  in  Figure  1. 

C .  Particle  Size  Parameter  and  Ionization  Rate  Contours 

As  mentioned  previously  under  the  basic  assumptions,  the  fission 
products  are  either  fused  Inside  or  adsorbed  on  the  outside  of  the  dust 
particles  and  are  inhaled  into  the  fireball.  The  ultimate  ionization  rate 
contours,  or  fallout  patterrs,  on  the  ground  are  related  to  the  size  of  the 
failing  particles.  The  particle  size  parsuneter,  ct,  is  defined  as  the  ratio 
of  wind  ve3.oclty,  V^,  to  the  average  falling  velocity  of  the  particle,  V^, 
since  it  is  more  concerned  with  the  falling  rate  than  the  actual  physical 


size  of  a  particle. 


(COORDINATE  ORIGIN  AT  GROUND  ZERO)  (COORDINATE  ORIGIN  AT  CLOUD  CENTER) 


Figure  1.  Geometry  of  Stem  aod  Cloud 
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For  the  stem,  it  is  assumed  that  the  particles  at  a  given  altitude  in 

the  rising  fireball  are  concentrated  on  a  horizontal  plane  through  the  center 

of  the  fireball.  Therefore,  all  particles  which  have  equal  values  of  a  fall 

out  of  the  stem  starting  at  the  same  altitude.  Any  particl*'  starts  to  fall 

when  its  instantaneous  falling  rate,  V^,  equals  the  rate  of  rise  of  the  sur- 
dZ 

rounding  air  mass,  There  is  an  empirical  relation  between  V„,  V„.  and 

the  fireball  altitude  Z; 

V2/V^  =  p  +  qZ  (7) 

where  p  =  0.95>  q  =  1.02  x  10  where  Z  =  5j000  to  50>000  ft 

with  particle  diameter  d  =  200  to  1,200  microns 
p  =  0.58,  q  *  l.?'^  X  10'^,  when  Z  =  50,000  to  110,000  ft 

with  particle  diameter  d  =  3OO  to  1,000  microns 


There  is  also  '’.n  empirical  approximation  for  Z: 


Z  -  (1  -  e  ^  ), 


for  t  =  20  to  480  seconds  (8) 


where  Z^  is  a  yield  dependent  multiplier 


Thus, 


is  0.011  sec~^. 


(9) 


For  values  of  Z^,  see  Figure  5. 
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Therefore,  the  particle  size  parameter,  CL,  for  the  particle  group 
falling  from  the  altitude  Z  is 


a  = 


V^(p  +  qz) 


(10) 


Moreover,  it  can  he  detemnined  from  geometry  of  the  stem  that 


a  Z  =  X  *  a„  (  +  for  a  .  ; 
m  Z  min 


for  CL  ) 
max 


(11) 


Therefore,  a  and  a  .  may  he  obtained  for  the  particles  falling 
’  max  min 

at  a  distance  X  from  ground  zero.  In  addition,  the  rising  and  falling  time 
for  particles  of  a  given  value  of  ct  can  he  determined  from  combining  Equa¬ 
tions  (8)  and  (lO),  as  shown  below 


,  ak„Z  +  V  qZ 

^  /...\  1  ,_fZo  v^oi* 


(12) 


(sec) 


aCca^z^  -  v„p) 


(13) 


The  outer  dimensions  of  the  cloud  are  defined  by 


2  2  2 

2  ^2  ^ 

a  b 


(ll) 


Consequently,  ct  and  a  ,  for  the  downwind  distance  X  can  be 
^  max  min 


derived,  as  shown  by  tl>*  aqwttlon 

*  (1 . 


GL 


hX 


)  +  (a^  -  y^)  [h^  -  b^(l  -  ^  /a^)] 


b^(l  -  y  lef) 


(15) 


This  calculation  usually  Includes  the  180  sec.  average  delay  time  due 
to  fireball  toroidal  circulation,  which  should  be  subtracted  in  order 
to  obtain  the  actual  rising  time. 
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(where  the  sign  is  +  for  a  ,  and  -  for  ct  .  ) 

'  max’  min' 

Since  the  cloud  is  assumed  to  be  stabilized,  all  particles  with 
the  same  values  of  a  fall  along  the  same  slope.  Hence,  the  times  of  arrival 
and  cessation  for  particles  of  a  given  oc  can  be  determined  from  the  georaetry 
of  the  cloud 


cc  (h  +  z  )  /  V 
a  w 

(16) 

a  (h  +  z  )  /  V 
c  '  w 

(1?) 

vhere  z  and  z  are  the  smallest  and  largest  intercepts  of  the  elope  to 

8*  C 

the'  outer  dimension  of  the  cloud,  respectively. 


a(X  -  Qh)  b^  i  ab  H 


o^b^) 


(1  -  y^7b^)  -  (X  -  Ohr 


2 

a  + 


a^b^ 


(18) 
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(wiiere  the  sign  is  +  for  and  -  for 

Assuming  that  the  ionization  rate,  or  fallout  intensity, 
varies  exponentially  with  downwind  distance,  and  its  distribution  is  closely 
related  to  a  and  W,  also  confirmed  by  empirical  test  data,  a  complete  set 
of  scaling  functions  which  defines  the  variation  of  intensity  versus  down¬ 
wind  distance  and  the  bi-elllptical  outline  of  the  ionization  rate  contours 
is  tabulated  in  Chapter  3  of  "Fallout  and  Radiological  Countermeasures''^^^  . 
Those  scaling  functions  are  not  repeated  in  this  report,  but  a  corresponding 
computer  program  was  established,  as  presented  in  Figure  5-  A  typical 
ionization  rate  contour  for  a  ID  MT  weapon  yield  and  a  15  mph  wind  velocity 
is  shown  in  Figure  6 .  Since  the  only  inputs  are  weapon  yield  and  wind 
velocity,  this  computer  program  is  relatively  simple  to  understand  and 
to  use  hy  others  to  obtain  the  information  on  ionization  rate  contours. 

It  should  he  pointed  out  that  these  contours  have  been  corrected 
to  a  standard  reference  time,  H  +  1  hour,  by  a  decay  correction  factor 

I  (t;  actual  time) 

"  I  (l;  H  +  1  hour}  * 

XL 


This  decay  correction  factor  can  be  obtained  directly  from  a 
typical  close-in  fai],out  decay  curve  as  shown  in  Figure  2  . 

Knowing  the  ionization  rate,  the  exposure  dose,  in  the  cloud 
fallout  area,  which  is  of  greater  interest  than  the  stem  fallout  area,  can 
be  estimated  by 


Ix(t)  (it 


(20) 


(id  os/Noissid  y3d  yH/a  J^_ol)  diva-NoiivziNoi 


Typical  Ionization  Rate  Decay  Cuive  for  Veiy  Close-In  Fallout 


D.  Condensation  of  Fallout  and  Fractionation  Numbers 

During  a  nuclear  detonation,  a  tremendous  amount  of  heat  is  created 
by  nuclear  fission  and  fusion.  This  heat  melts  and  also  vaporizes  eveiy- 
thing  inlialed  into  the  air  mass  from  the  ground  and  forms  the  familiar 
fireball.  When  the  fireball  rises,  its  volume  begins  to  expand,  and  its 
temperature  is  assumed  to  drop  according  to  the  ideal  gas  law.  The  material 
swept  inSide  the  fireball  starts  to  liquidify  and  finally  solidify  to 
particles,  in  or  on  which  the  fission  products  begin  to  condense.  The 
condensation  process  is  generally  divided  into  two  time  periods  as  mentioned 
in  the  basic  assumptions.  The  dividing  temperature  depends  largely  on  the 
composition  of  the  fireball.  For  a  model  surface  burst  on  a  soil  consisting 
of  silicate  minerals  this  temperature  is  assumed  to  be  l673°K,  as  most  ele¬ 
ments  will  have  started  to  solidify  at  this  temperature.  Since  the  two 
periods  are  characterized  by  the  coexistence  of  either  gas  and  liquid,  or 
gas  and  solid,  the  vaporization  or  sublimation  pressure  plays  an  Important 
role  in  the  condensation  process.  The  following  formula  for  the  evalua¬ 
tion  of  the  pressure  is  presented  in  "Fallout  and  Radiological  Counter- 
measures  '  ' : 


log  Pj  =  I  +  B  +  C 


(21) 


where  p  is  the  vaporization  or  sublimation  pressure  for  element  j, 

<J 

gm/sq.  cm., 

T  is  temperature  of  the  fireball,  °K, 

A,  B,  C  are  the  empirical  constants  for  vaporization  or  sublimation 

reactions  of  elements  contained  in  the  fireball  which  were 

(k) 

tabulated  in  Miller's  manuscript  . 

The  established  computer  program  for  estimating  is  shown  in 


Figure  3 • 
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Figure  3.  Computer  Program  for  Estimating  Sublimation  Pressure 

B*n-2?n  sT*NtJ*«n  version 

COMMENT  SUBttMiTION  PRESSURL-  ,*•  PtJ  RY  T.  E.  JiMFS.  % 

INTEGER  J.If ABC»DEE.JK^4.  *  i^RAY  Q(7<>  »  TI70>  S 

READitiEAA)  S 

WRITEIitEAdiEAC)  S 

START..  PEAOlJtOATA)  1 

IF  mOD(ASC»27>  EQL  0  S 

HRITEI  fTEAU.EACJ . .  .  - 


02^0  FOR  I  a  (l.l.OEF)  t 

0241  (l<n  •  EXP<  <2.302965)(  A/T<  1 )  ♦  R  *  C)>  ..I 

02M  .  uRITE<FSE4U*6A6»  _  .  .  ..  * 

0,2A8 _  ,  GO  START  $  craft 

B2.A.3. .  UPUT  0ATA(A9C  ,  a  ,  B  ,  ,C  JK A ,  JKp, JKo , J<0  S 

02S1 _ INPUT  I  *  (1»1»DEF)  $  T«.I1.»  .  _  .  % 

0315  OUTPUT  E*®‘F0«  .I..»  .(t».l*.0£F_>  .$  JI  I.>  >  ....  * 

03‘'6  output  EAD<  JKA.  JK9.JK0...I.KD,.A.8.C.0(1)  )  * 

03''4..  . _ FORMAT  p  Ar  (  HO .  •  lOENT  1 F  I  CATl  ON».  BTl.  t  A  t ,  81  Q  .  »B  a.  BU  .  »e«  .  H9  .  _ 

03*4  _•  PJ  FOR  1*73  .PfG..*».I_T0_._W3L _  * 

03A3  format  FAe<4A5.X12.0*XH1,3.X13..6,F?O.S.W4»  t 

0391  . .  FINISH  _ _ _  _ L.CRAFT 


compiled  PROGRAM  ENnS  AT  0392 
PROGRAM  variables  regin  AT  4230 


0200 

0200 

0200 

0200 

0204 

0212 

0216 


Figure  3  (Continued).  Computer  Program  for  Estimating  Sublimation  Pressure 
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Figure  3  (Continued).  Computer  Program  for  Estimsting  Sublimation  Pressure 
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Figure  3  (Continued).  Computer  Program  for  Estimating  Sublimation  Pressure 
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Pigure  3  (Continued).  Computer  Program  for  Estimating  Sublimation  Pressure 
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The  elements  in  the  fireball  are  radioactive  and  hence  the  number 
of  moles  of  each  species  is  constantly  changing  vith  time.  If  y..(t) 
denotes  the  amount  of  radionuclides  of  element  J  and  mass  number  A  present 
at  time  t  after  fission,  then  the  total  amount  of  element  J  present  at  time  t 
is 


r.(t) 

J 


atoms  or  moles 


(22) 


Correspondingly,  the  total  amount  for  the  entire  chain  1  of 
mass  number  A  is 


h  -- 


atoms  or  moles 


(23) 


in  which  is  constant,  except  for  mass  chains  containing  neutron  emitters. 

For  the  first  period  of  condensation,  the  material  balance  for 
element  J  between  gas  and  liquid  phase  is 


{2k) 


where  n..°(t)  is  the  number  of  moles  of  element  J  with  mass  number  A, 

J  A 

which  is  mixed  with  moles  from  other  mass  chains  to 


form  n  moles  of  vapor,  and  is  abbreviated  to  n  ,  and 

J 

n  (t)  is  the  number  of  moles  of  eleir.ent  j  wita  mass  number  A 

JA 

dissolvea  in  the  ri(-^)  moles  of  liquid  carrier,  which 
is  the  paroicle  in  the  liquid  phase,  prior  to  solidi¬ 
fication,  and  is  abbreviated  to  n.. 

J 

For  a  perfect  gas,  the  ratio  between  the  mole  fraction  of  element  i, 


in  the  vapor  pnase  and  that  in  the  liquid  phase  is 
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11°  n?/a  n°  ET/p.V  k, 

.1  0  d  J 

iij  =  -  -aynCir  == 


(25) 


where 


is 

the 

is 

the 

n 

is 

the 

n(£) 

is 

the 

R 

is 

the 

1 

X 

is 

the 

is 

the 

V 

is 

the 

is 

the 

mole  fraction  of  element  j  in  the  vanor  phase 

mole  fraction  of  element  j  in  the  liquid  phase 

number  of  moles  of  vapor 

total  moles  of  liquid  carrier 

molar  gas  constant 

absolute  temperature 

total  vaporization  pressure 

molar  volume ,  and 

Henry's  Law  constant 


From  Equation  (25),  the  following  relation  between  n°  and  n^  can  be 

obtained; 


o 

n. 

J 


n.k. 


[d(4/V 


V]  RT 


where  k°  =  k  /[n(i)/V]  RT 
J  J 


(26) 


Substituting  Equation  (26)  into  Equation  (2^4-)  yields 


(27) 


The  fractionation  number  of  the  first  period  of  condensation, 

r  (A,t),  is  defined  to  be  the  fraction  of  the  mass  chain  condensed,  or  in 
o 

mathematical  form 


0 


(28) 
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Agaln  according  to  Henry's  Law 

J  “j 


where 


is  the  fugaclty  of  the  element  in  the  liquid  phase 

is  the  mole  fraction  of  the  element  in  the  liquid  phase ^  and 

is  the  partial  pressure  of  the  liquid. 


However,  by  Raoult's  Law  for  an  ideal  solution 


Therefore 


The  method  for  determining  p.  has  been  presented  in  Equation  ( 21) . 

J 

The  [n(.<i)/V']  ratio  is  given  hy  a  scaling  function 

[n(f)/v]  =  5«67  X  10"'^exp  [-0.510  £12i|£  (;■ 


where  t,  in  seconds,  is  related  to  the  fireball  temperature,  T,  by 

T  =  1.66  X  10^  exp( -0.546  (33) 

and  the  data  of  y..(t)  were  calculated  and  tabulated  as  N.(A,  t)  by  Holies 
and  Ballou^  .  Hence,  the  fractionation  number,  f  can  be  evaluated  for 

the  elements  contained  In  the  fireball. 

During  the  second  period  of  condensation,  if  an  excess  of  solid 
surface  area  is  present,-  the  amount  condensed  at  any  time  after  solidifica¬ 
tion  of  the  carrier,  is  given  by 
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n!  =  -  n" 

i  3  j 


{3h) 


vjhere 


is  the  amount  of  element  j  condensed  on  the  surface  of 
the  solid  particles 

is  the  amount  of  element  .j  that  has  not  condensed  In  the 
liquid  carrier,  and 

is  the  amount  of  element  j  in  the  vapor  phase. 


For  a  perfect  gas 


n. 


KC 


(35) 


wnere 


is  the  sublimation  pressure  and  may  be  evaluated  according 
to  Equation  (2l). 


The  material  balance  for  element  j  is 


yj(t)  rr:  nj(t)  +  n^(t)  +  n^(t) 


(36) 


Tnerefore,  employing  Equations  (2T)  and  (35) 


n. 

J 


i  +  k 


o 


(37) 


or 


k° 


VjA  ‘  .  \o  VjA  '  RT  VjA 

^  ‘"j 


(38) 


In  this  equation  the  partial  pressure  of  each  nuclide  is  propor¬ 
tional  to  its  abundance  at  any  given  time,  and  is  given  by 


(39) 
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where 


.  ..o  „  o 

=  VoA 


The  fractionation  number  of  the  second  period  of  condensation 
t)  is  defined  to  be  the  fraction  of  the  mass  chain  condensed  up  to 
a  gi'ven  time  during  the  second  period,  or  expressed  mathematically 


r  ■  '  A ,  t ) 

c  ■  ’ 


'A 


V  ^  s  1 

0.23  :EWF.T 


(40) 


or 


[A,t)  =  1 


(A,  t)  - 


V 


0.23  hWET 


4^) 


vhere  Y 


is  the  fireball  volume  and  is  given  by  scaling  the  function 
V  7.62  X  10^^  W  exp  (O.^iO  cra^ 


B  is  the  ratio  of  fission  to  total  yield,  and 

0,23  i&  a  conversion  constant  changing  weapon  yield,  W,  to 

number  of  fission-moles  per  atom. 


Physically  in  Equation  (hi),  v7rT  gives  the  number  of  moles  of  element  j 

t3 

in  the  vapor  phase  and  0.2}  tVTi  ?  gives  the  total  number  of  moles  of  the 

j 

element  produced  .and  in  existence  at  a  given  time  in  the  second  period, 
i'he  combined  ti-rai  in  Eq.^ation  f'hi)  gii-eE  the  fraction  still  remaining  In 
the  vapot  Elate  during  the  second  period  of  condensation.  The  1  - 
t-f-TTn  gi'.'e.’  the  fraction  w'hich  has  not  condensed  in  the  liquid  carriers  during 
their  existence  in  the  first  period  of  condensation. 
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The  fractionation  nimbers  yield  the  information  ahout  the  fractions 
of  mdioactivity  condensed  in,  or  on  the  particles,  at  any  time  after  a  hurst. 
Therefore,  the  fractionation  numbers  not  only  make  the  conversion  between 
intensity  and  activity  possible,  but  also  lead  the  way  to  the  evaluation  of 
soluble  and  insoluble  activity  which  is  importaiit  to  biological  availability 
and  radiological  counts measxires  research. 

E.  Conversion  between  Intensity  and  Activity 

The  conversion  between  intensity  (measured  as  ionisation  -rate)  and 
activity  density  in  fissions  per  tinit  area  is  presented  in  "Fallout  and 
Radiological  Countermeasures'',^^^  based  on  the  simple  ass’jmption  that 


I^(t)  I^(t)(ideal) 
“A^  “  A^(ideal)  " 


(U2) 


where  the  conversion  factor,  K^(t),  is  further  broken  down  into  the  fol¬ 
lowing  components 

Kj^(t)  =  (^3) 

in  which 

D  is  an  Instrument  response  factor,  usually  assigned  a 

value  of  0.75 

q  is  the  terrain  shielding  factor,  usually  assigned  a 

A 

value  of  0.75 

r^Ct)  is  the  gross  fission  product  fraction  number,  sometimes 
also  designated  as  rQ,(t) 

i^^(t)  is  the  air  ionization  rate  per  fission  at  3  ft  above 

an  infinite,  ideal  plane  for  a  uniform  distribution  of 
the  normal  fission  product  mixture,  and 
i^(t)  is  the  same  unit  for  neutron  induced  activities. 
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Details  about  the  neutron  induced  ionisation  rate,  i^(t),  are  almost  entirely 

classified,  hut  generally  assigned^^^  as  1-9^  of  i„  (t)  The  other  components 

in  equation  (^3)  can  be  readily  derived. 

Detailed  information  about  the  activities  of  the  U-235  fission 

( 7) 

products  has  been  provided  by  Bolles  and  Ballou  .  These  data  can  be 

/  Q  N 

correlated  to  the  data  supplied  by  Katcoff'  and  applied  to  other  types 
of  fission  elements,  such  as  U-238,  P\.i-239  ty  taking  proportions  of  the 
corresponding  mass  number.  This  conversion  has  been  systematized  and 
tabulated  in  "Fallout  and  Radiological  Countermeasures".  Moreover,  Miller^^^ 
translated  the  activity  data  from  Bolles  and  Ballou  into  ionization  rate  by 
a  disintegration  mialtiplier,  m.  Both  tables  of  correlation  and  of  disinte¬ 
gration  multipliers  are  reproduced  from  Miller's  manuscript  as  Tables  I 
and  II.  Therefore,  the  ideal  or  normal  ionization  rate  of  fission  product 
mixture  may  be  computed  according  to  the  following  formula: 

i^p(t)  =  (44) 

4 

where  is  the  activity  of  the  nuclide  at  time,  t,  per  10  fission. 

In  accordance  with  the  earlier  discussion,  the  fission  products 
do  condense  and  fractionate  with  time  instead  of  condensing  completely  at 
the  very  beginning  as  in  the  ideal  situation.  Therefore,  the  real  ioniza¬ 
tion  rate,  l^*(t),  should  be  multiplied  individually  for  each  radionuclide 
by  its  fractionation  number  at  any  specified  time.  Mathematically,  this 
may  be  expressed  by 

l^*(t)  =  Bm  A^[r^(A,  t)  +  r^(A,  t)]  (4?) 

The  r^  (A,  t)  term  vanishes  in  the  first  period  of  condensation. 

As  time  increases,  the  sum  of  the  two  fractionation  numbers  approaches  unity. 
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TABLE  1 

(4) 

Cumulative  Mass-Chain  Yields  of  Fission  Products^ 
(Values  are  in  per  cent  of  fissions) 


Mass 

u235 

u238 

Pu239 

Wumber 

Thermal 

Neutrons* 

Fission 

Neutrons 

Fission 

Neutrons 

8-Mev 

Neutrons 

Thermal 

Neutrons 

Fission 

Neutron; 

72 

1.6x10"^ 

4.6x10“^ 

5.0x10"^ 

_ 

-4* 

1.2x10 

73 

-k 

I.LkIO 

0.0012 

3.7.x10“5 

- 

2.2x10"^ 

- 

7k 

(3.2x10"^)°' 

0.0034 

-4 

1.1x10 

0.001 

1  -V 

4.1x10 

0.0011 

75 

(8. 8x10“^) 

0.0062 

8.3x10"^ 

0.0040 

-4 

T.xlO 

0.0023 

76 

(O.CO29) 

0.012 

0.0012 

0.0078 

0.0014 

0.0051 

77 

0.0083 

0.023 

0,0038* 

0.014 

0.0026 

0.011 

78 

0.021 

0.048 

0.0095 

0.026 

0.0049 

0.025 

79 

(0.04l) 

0.096 

0.019 

0.053 

0.0090 

0.043 

3o 

(0.077) 

0.19 

0.045 

0.096 

0.016 

0.075 

81 

0.l4 

0.21 

0.088 

0.18 

0.030 

0.14 

82 

(0.29) 

0.50 

0.20 

0.35 

0.056 

0.2,3 

83 

0.5H 

0.80 

0.4o* 

0 . 66 

0.10 

0.37 

81^ 

1.00 

1.3 

0.85 

1.02 

0.17 

0.60 

85 

1.30 

1.85 

0.80 

1.45 

0.28 

0.92 

86 

2.02 

2.5 

1.38* 

1.9 

0,45 

1.15 

87 

(2.9h) 

3.3 

1.90 

2.25 

0.73 

1.5 

88 

(3.92) 

4.2 

2.45 

2.7 

1.2 

1.9 

89 

1^.79 

5.1 

2.9* 

3-17 

1.9* 

2.4 

Continued 
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TABtE  I.  vcontinaed) 

Cumulative  Mass-Chain  Yields  of  Fj  ssion  Products 


(j'alues  are 

in  per  cent 

of  fissions) 

Mass 

.^35 

.238 

Numher 

Thermal 

neutrons^ 

Fission 

Meutrons 

Fission 

ISeutrons 

8-Mev 

Oieutrons 

Thermal 

Neutrons 

Fission 

Neutrons 

90 

5^77 

5.8 

3.2* 

3 .7 

2.1- 

3.0 

91 

5.8h 

5.85 

3.6 

4.3 

3,0 

3.7 

92 

6.0'i 

6.0 

i|-.l 

k.6 

3.7 

4.4 

93 

6.^5 

6 .1- 

i^.85 

5.2 

4.6 

5.0 

9k 

6.k0 

6.1 

5.3 

5.1-5 

5.5 

5.4 

95 

6.27 

6.3 

5.7^ 

5.6 

5..9*' 

5.6 

96 

6.33 

6.3 

5.8 

5.7 

5.7 

5.3 

97 

6.09 

6.1 

5.7 

5.61 

5.6* 

5.2* 

98 

9.78 

5.8 

5.7 

5.6 

5.4 

5-4 

99 

6.06 

6 . 

6. 

6.2** 

5-9* 

5.9"' 

100 

6.30 

6.7 

6.1 

6.1 

6.0 

6.4 

101 

5.0 

5.5 

5.5 

6.5 

6.0 

5.9 

102 

1.1 

2.9 

5.6 

5.9 

5.9 

5.3 

103 

5.0 

1.7 

6.6 

5.0 

5.8* 

4.6 

lot 

1.8 

0.95 

5.1 

3-2 

5.0 

3.5 

105 

o.go 

0.54 

3-9 

2.2 

3.9" 

3.2 

l06 

0.38 

0.30 

2.7* 

1.5 

5.0* 

3.6 

lOT 

0.19 

0.17 

1.35 

1.0 

1,0 

3-1 

i08 

70. 085) 

0.095 

0.67 

0.70 

3.0 

2.6 

109 

(0.039, 

0.053*^^^ 

0.32* 

0.18 

1.5* 

1.9* 

no 

(0.020) 

0.030 

0.15 

0.33 

0,65 

0.8i 

111 

(0.015; 

0.022-*»* 

0.073* 

0.23*-* 

0.27* 

o.3i. 

112 

fO.Olj/ 

0 ,020*** 

0,016* 

0.19 

0.10* 

0.14* 

111 

1  0.012) 

0.018 

0.0i4  3 

0.17 

0.055 

0 . 090 
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TABLE  I  (continued) 

Ciomulative  Mass-Chain  Yields  of  Fission  Products 


(Values  are 

in  per  cent 

of  fissions) 

Mass 

u235 

^J238 

Pu239 

Kumher 

Thermal^ 

Neutrons 

Fission 

Neutrons 

Fission 

Neutrons 

8-Mev 

Neutrons 

Thermal 

Neutrons 

Fission 

-Teutrons 

11^4- 

(0.011) 

0.017 

0.04l 

0.16 

0.046 

0.075 

115 

0.0104 

0.017*** 

o.o4o* 

0 . 15*** 

0.041* 

0.069* 

ll6 

(0.010)^ 

0.017^ 

0.039 

0.l4 

0.039 

0.065 

U7 

(0.010) 

0.017 

0.039 

0.14^ 

0.038 

0,065 

118 

(0.010) 

0.017 

o.4q^ 

0.14 

0.038^ 

0.064^ 

119 

(0.011) 

0.017 

0.041 

0.l4 

0.039 

0.064 

120 

(0.011) 

0.018 

0.042 

0.15 

0.041 

0.065 

121 

(0.012) 

0.020 

0.044 

0.16 

0.044* 

0.066 

122 

(0.013) 

0.022 

0.046 

0.17 

0.04t 

0.069 

123 

(0.015) 

0.030 

0.050 

0.19 

0.052 

0.076 

121 

(0.017) 

0.053 

0.055 

0.23 

0.058 

0.082 

125 

0.021 

0.095 

0.072 

0.33 

0.072* 

0.14 

126 

(0.058) 

0.17 

0.175 

0.48 

0.175 

0.35 

127 

(0.145) 

0.30 

0.39 

0.70 

0.39^ 

0.60 

128 

0.37 

0.54 

0.77 

1.0 

0.77 

1.9 

129 

0.90 

0.95 

1.45 

1.5 

1.45 

2.5 

130 

2.0 

1-7 

2.5 

2.2 

2.5 

3.2 

131 

(2.88) 

2.9 

3.2* 

3.2 

3.8* 

3.8 

132 

(4.31) 

4.7* 

4.4 

5.0 

4.6 

133 

(6.48) 

6.1 

5.5^ 

-.4 

✓  * 

5.27- 

4.9 

13I 

(7.80) 

7-3 

6.6* 

6.5 

5.69> 

5.2 

135 

(6.40) 

6.3 

6.0* 

5.9 

5.53^ 

5.1 

136 

(6.36) 

6.4 

5.9* 

5.8 

5.06* 

5.3 

137 

(6.05) 

6.0 

6.2 

5.85 

5.24* 

6.4*t 

Continued 


-30- 


TABLE  I  (co£LtlE\^ec3  ; 

Cumai.ative  Maes-Chain  Yields  of  Fission  Products 
(/alues  a, re  in  per  cent  of  fieeions) 


Mai 


..,235 


^,238 


Pu 


239 


i'T'uitiber 

Thermal 

Neutrons. 

Fission 

Neutrons 

Fission 

Neutrons 

8-Mev 

Neutrons 

Thermal 

INe'utrone 

Fission 

Neutrons 

138  '. 

5.74 

5.7 

6.4 

5.9 

5-5 

5-4 

139 

(6.34) 

6.4 

6.5 

6.0 

5.7* 

5.2 

II.O 

6.1+4 

6.4 

5.7^ 

5.6 

5.68'* 

5.0* 

II.I 

(6.30) 

6.3 

5.7 

5.5 

5.2* 

4.7 

1.12, 

(5.85) 

5^9 

5.7 

5.4 

6.69* 

4.9 

ii^3 

(5.87) 

5.8 

5.5 

4.97 

5.4* 

5.0 

5.67 

5.1^* 

4.9^. 

4,3-** 

5.29* 

4.8 

145 

.3.95 

4  .2 

3.7 

3.7 

4.24* 

4.4 

146 

3.07 

3.3 

3-1 

3.17 

3.53^ 

3.7 

1“  7 

2,38 

2.5*^* 

2 . 6*^ 

2.7** 

2.92* 

3.0 

l48 

1.70 

1.85 

2.0 

2.27 

2.28* 

2.36 

l'-9 

1.13 

1.3**-^' 

1.45 

1.9^* 

1.75 

1.86 

150 

0.67 

0.80 

1.05 

1.45 

1.38* 

1.48 

151 

0.1+5 

0.50 

0.74 

1.02 

1,08 

1.16 

i;-2 

0.265 

0.31 

0.50 

0.66 

0.83* 

0.92 

153 

0.15 

0.19*-^ 

0.32 

0.41** 

0.52 

0.60 

154 

0.077 

0.096 

0.19 

0.25 

0.32* 

0.37 

155 

0.033 

0.048 

0.11 

0.15 

0,20 

0.23 

^56 

0.014 

0.02  3^** 

0,066* 

0.092** 

0.12* 

0.l4 

157 

0.0073 

0.012 

0.034 

0.057 

0.064 

0.075 

0.002 

0.0062 

0.016 

0.032 

0.034 

0.043 

159 

0.00107 

0.0034*'^' 

0.0090** 

0.017** 

0.020**** 

0.025 

Cor.'*,  ir.uei 
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TABLE  I  (continued) 

Cumulative  Mass-Chain  Yields  of  Fission  Products 
(Values  are  in  per  cent  of  fissions) 


Mass 

u235 

y38 

Pu239 

Number 

Thermal  Fission 

Neutrons*  Neutrons 

Fission 

Neutrons 

8-Mev 

Neutrons 

Thermal 

Neutrons 

Fission 

Neutrons 

160 

3.5x10-4 

0.0012 

0.0036 

0.0085 

0.0092 

0.011 

161 

7.6x10-5 

4 . 6x10-4** 

9.4x10-4 

0.0044** 

0,0038**** 

0.0051 

*Seymour  Katcoff,  Fission- Product  Yields  From  U,  Th  and  Pu,  Nucleonics, 

Vol.  16,  No.  !+,  p  7^5-85  (195B) 

**■1.  R.  Bunney,  E.  M.  Scadden,  J.  0.  Abriam  and  N.  E.  Ballou,  Radiochemical 
Studies  of  the  Fast  Neutron  Fission  of  and  Second  'JN  Inter¬ 

national  Conference  on  the  Peaceful  Uses  of  Atomic  Energy,  A/Conf.  15/F/6i4-3, 
USA,  June  1958 

**+G.  P.  Ford,  J,  S.  Gilmore,  et  al.  Fission  Yields,  LADC-3083,  1958 

R.  Bunney,  E.  M.  Scadden,  J.  0.  Abriam,  N.  E.  Ballou,  Fission  Yields 
in  Neutron  Fission  of  Pu^39^  USNEDL-TR-268,  1953,  Unci. 

a.  Parentheses  indicate  estimated  values  or  where  Katcoff's  value  was  altered 
in  order  to  adjust  the  yields  to  a  gross  sum  of  100  in  each  neak. 


b .  Line  indicates  division  of  two  peaks  that  was  used  for  individual  peak 
sums . 
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a?ABLE  II 


Air  lonigatlon  Disintegration  Mixltlpliers  for., 
the  Fission  Product  and  Other  Eadionuclides 


Nuclide 

10"^  r/hr-ft^ 

Nuclide 

10"^  r/hr-ft' 

(dis/sec ) 

(dis/sec ) 

Zn-72 

6.75 

Nb  -95111 

1.4l 

Zn-7it- 

4,22 

Nbp-95 

4,56 

Nb^-97m 

4.42 

Ga-l^ 

14.3 

Nfc,-97 

4.02 

Ga-73 

2.11 

Nb-98 

11.7 

Ga-74 

12.4 

Mo -99 

0.772 

Ge-75 

0.216 

Mo-101 

6,85 

Ge-77 

13.1 

Mo-102 

0.0 

Ge-78 

2.74 

Tc-99ni 

0.782 

As -77 

0.0735 

Tc-101 

2.03 

As-78 

2.51 

Tc-102 

2.45 

A3-79 

0.0 

Ru-103 

2.99 

Se^ -8lm 

0.117 

Ru-105 

4.36 

Se^-8l 

0.0 

Ru-106 

0.0 

Se-83 

12.8 

Rh-103in 

0.150 

Br-83 

0.017^ 

Rh  -105m 

0.283 

Br-84 

8.12 

Rh;-105 

0.559 

Rh^lOe 

1.81 

Kr-83ru 

0.0959 

Rh-lO? 

5.08 

Kr  -85m 

0.973 

Krp-85 

0.0316 

Pd -109 

0.0 

Kr-87 

7.03 

Pd-111 

5.00 

Kr-88 

7.40 

Pd-112 

0.0225 

Rh-88 

4.30 

Ag-109ni 

0.172 

Rh-89 

12.0 

Ag-111 

0.175 

Rh-91 

6.67 

Ag-112 

4.59 

Ag-113 

0.158 

Sr-89 

0.0 

Ag-115 

1.28 

Sr-90 

0.0 

Sr-91 

5.03 

cd  .115 

0.187 

Sr-92 

6.93 

cd;:-ii5 

1.11 

Sr-93 

1.20 

Cd-117m 

8.11 

Cd-118 

0.0 

Y-90 

0.008 

Cd-120 

3.08 

Y  -91in 

3.24 

Yi-91 

0.0135 

In-115 

1.12 

Y-92 

6,56 

in-117 

1.36 

Y-93 

0.817 

In-ll8 

13.5 

Y-9^ 

6.67 

In-119 

0.0371 

Zr-95 

4.32 

Sn-121 

0.0 

Zr-97 

0.346 

Sn.123 

0.0 

Sn-125 

6.19 

Sn-126 

0.0 

Sn-127 

3.08 
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Nuclide 

Table  II  (continued) 

Air  Ionization  Disintegration  Multipliers  for 

m 

10"^  r/hr-fi 

the  Fission  Product 

and  Other  Radionuclides 

m 

10" 5  r/hr-ft^ 

Nuclide 

(di s/sec ) 

(dis/sec 

Sb-125 

2.78 

Ce-l4l 

0.470 

Sb-126 

16.1 

Ce-l43 

2.34 

Sb-327 

2.70 

Ce-144 

0.188 

S'b-128 

13.0 

Ce-l45 

4.22 

Sb-129 

6.10 

Ce-l46 

1.60 

Sb-131 

3.66 

Pr-l43 

0.0 

Te, -125m 

0.0525 

Pr-l44 

0.163 

Ter>127 

O.lOl 

Pr-l45 

0.0 

TeJ-127 

0.0249 

Pr-l46 

6.42 

16^-12 9m 

0.398 

Tep-129 

1.52 

Nd-l47 

0.968 

Te^-131 

8.42 

Nd-l49 

2.31 

Te;-131 

2.08 

Nd-151 

6.44 

Te-132 

1.52 

^e, -133m 

0.790 

Pm-147 

0.0 

Tei-133 

9.12 

Pm-l49 

7.25 

Te-13i+ 

7-89 

Pm-150 

4.99 

Pm-151 

2.21 

1-131 

2.38 

Pm-152 

3.79 

1-132 

13.0 

Pm-153 

4.76 

1-133 

3-59 

1-13^^ 

9.63 

Sm-151 

0.0544 

1-135 

9.60 

Sni-153 

0.455 

Sm-155 

1.80 

X3-,  -131m 

0.207 

Sm-156 

1.20 

Xer-133ra 

0.413 

Sin-158 

3.08 

XeJ-133 

C.304 

Xe^-135m 

2.69 

Eu-155 

0.314 

Xe^-lSS 

1.55 

Eu-156 

5.55 

Xe-138 

7.89 

Eu-157 

3.78 

Eu-158 

7.21 

Cs-137 

0.0 

Cs-138 

11.2 

Cki-159 

0.430 

Cs-139 

4.78 

Tb-161 

0.123 

Ba-137m 

3.65 

Ba-139 

0.888 

U-237 

1.13 

Ba-lUb 

1.10 

U-239 

0.388 

Ba-l4l 

4.22 

U-240 

0.0 

Ba -1)4.2 

5.76 

Np-239 

1.11 

La-l)+.0 

13.0 

Np-240 

2.07 

La-li+l 

0.400 

1&-Ik2 

11  6 

Mn-56 

9.40 

La-ll4-3 

6.67 
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Ihe  gross  ioni:;ati'^n  fractionation  number,  is  defined  as 

the  ratio  between  the  ionization  rates  at  H  +  1  hour  of  condensed  and  normal 
fission  products.  It  may  be  expressed  as 


rjj(l) 


Thus,  all  the  components  for  the  coi.i  :  ~ion  factor,  ICjj(t),  have  been  cited 
and  derived.  The  intensity-activity  conversion,  shown  as  Equation  (42), 
is  very  useful  in  fallout  contamination  studies  and  the  determination  of 
other  prope.rties  of  fallout . 

F.  Biological  Availability  and  Atom-Concentration  Intensity  Ratio 

The  fraction  of  radionuclides  condensed  on  the  outside  of  fallout 
particles  Is  potentially  available  for  biological  uptake.  By  an  analogy 
to  the  intensity-activity  conversion  presented  in  the  foregoing  section, 
the  ratio  of  the  number  of  atoms  condensed  on  the  exterior  of  the 
particle  to  the  standard  intensity  at  H  +  1  hour  is  given  by  the  following 
foimula 


riTj  ■  K^(x} 


(47) 


where  nutaber  of  atoms  (of  the  end  member  of  the  mass 

chain  and  condensed  on  the  outside  of  the  particle) 
which  land  per  square  foot  of  ground. 

The  other  terms  in  Equation  (47)  are  the  same  as  defined  previously. 
The  ratio  N^I  is  called  the  atom-concentration  intensity  ratio. 

As  shown  in  Equation  (l2)  the  fireball  rising  time,  t.  Is  a  function 
of  the  particle-size  parameter,  a.  Because  of  time  dependence  this  ratio  is 
also  dependent  on  a.  However,  the  a  can  be  approximated  and  simplified  by  an 
average  value  Oi^,  which  is 
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«  =  I 

o  h 

Since  cn^  does  not  vary  with  cross-wind  distance,  the  atom-concentration 

intensity  ratio  holds  constant  as  a  first  approximation  for  all  values  of  cross- 

wind  distance  at  a  given  value  of  downwind  distance  Moreover,  as  shall 

he  seen  in  the  following  section,  a  given  value  of  a  represents  a  group  of 

*  / 

particles  with  a  small  range  in  diameters.  Therefore,  the  ratio  is  essen¬ 

tially  associated  with  a  given  paiticle  size. 

The  values  of  (A)/l(t)  for  six  hiologically  important  isotopes  at 
different  downwind  distances  for  a  typical  10  MT  hurst  have  heen  calci0.ated 
as  shown  in  Table  VI. 

This  ratio  can  he  applied  to  water  contamination  studies  to  yield 
the  activity  of  soluble  radionuclides  in  water  supplies.  If  r^(A,  t)  is 
replaced  by  r^(A,  t),  the  ratio  N^/E  will  provide  the  information  for  the  activity 
of  insoluble  radionuclides,  since  those  radionuclides  fased  into  the  silicate 
part, ides  are  assumed  to  be  not  readily  soluble  in  water. 

G.  Ionization  Rate  Contour  Ratios  and  Particle  Size  Distribution 

Several  ionization  rate  contour  ratios  have  been  defined  to  deter¬ 
mine  special  properties  of  fallout.  They  are: 

(a)  The  mass  contour  ratio,  M^(t),  which  is  the  ratio  of  the 
fallout  mass  per  unit  area  at  any  downwind  distance  to  the  ideal  ionization 
rate  at  that  location,  or  mathematically: 

M^(t)  - 

is  the  mass  of  fallout  per  unit  area  at  any  downwind 
distance,  X 

is  the  ideal  ionization- rate  at  3  ft  above  an  extended  open 


where 


"x 


yt) 


area  covered  with  fallout . 


rBBiB.-WKiiiy««g^^grfVTT«-yr!- 
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The  ma&s  of  fallout;  m,  empirically  determined  to  be  a  function 
of  both  a  and  W;  is 

m  =  f(a)  mg  (50) 


where 

f(a)  =  5-33  X  10"^  mg/fission. 

a 

=  0.1  to  0.5 

=  5-64  X  10"^  mg/fission. 

a 

=  0.9  to  2.0 

=  7.3.4  X  10“^^  mg/fission 

a 

>2.0 

According  to  Equation  {^2),  it  is  found  that 


1^(1)  X  Area  =  D(l)q^  [r^(l)  +  i^Cl)!  BW  (5l) 

Therefore,  the  mass  conto\ar  ratio  associated  with  H  +  1  hour,  M^(l), 
for  a  ground  surface  burst  will  be 


w  _  _ f(Q)  _ mg/sg  ft 

”  D(l)  qyB[r^(l)  ij.p(l)  +  r/hr  at  1  hr 


(52) 


(b)  The  fission-product  contour  ratio,  FP^(t),  which  is  the  ratio 
of  the  number  of  atoms,  or  moles,  of  fission  products  per  unit  area  at  any 
downwind  distance  to  the  ideal  ionization  rate  at  that  location,  or  mathe¬ 
matic  al].y; 

FP^(t)  = 


where 


is  the  number  of  atoms, 
per  unit  area. 

Prom  the  discussion  of  Section  6., 


or  moles,  of  fission  products 


this  contour  ratio,  corrected 


to  H  +  1  hour,  for  a  ground  surface  burst  may  be  expressed  as 
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1.16 

‘  — D(i')  q”tr„(ij'i";'n:r'i'iiiiT]" 

(c)  The  fraction  of  device  contour  ratio,  FD^(t),  vhich  ig  the 
ratio  of  the  fraction  of  weapon  device  per  unit  area  at  any  downwind  distance 
to  the  ideal  ionization  rate  at  that  location,  lliis  ratio  when  corrected  to 
H  +  1  hour  may  he  expressed  as 

-  D(1)  “*  " 

These  contour  ratios  predict  various  characteristics  of'  local 
fallout,  as  they  yield  information  about  the  mass,  the  number  of  atoms  of 
fission  prodacts  per  unit  area,  the  fraction  of  weapon  yield  per  unit  area 
and  other  properties  of  fallout  at  a  specific  location.  This  information 
is  valuable  for  the  design  of  fallout  shelters  and  other  radiological 
countermeasures.  Moreover,  the  use  of  these  contour  ratios  and  their 
scaling  functions  makes  possible  the  extrapolation  of  a  limited  amount 
of  experimental  data  to  cover  a  number  of  operational  cases  since  the  con¬ 
tour  ratios  are  not  constant  but  are  point  functions  whose  values  depend 
on  many  variables.  Furthermore,  the  average  density  of  the  fission  products 

can  be  estimated  from  the  ratio  of  FP  (l)  to  M  (l). 

r  r 

Another  advantage  of  the  O.C.D.  fallout  model  Is  that  particle 
size  groups  for  any  given  downwind  location  can  be  estimated  by  a  graphical 
method,  provided  that  the  particle  faxl-rate  is  known.  Since  the  fall-rate 
for  a  given  size  parameter  ct  does  denote  a  group  c  particles  with  a  range 
of  diameters,  with  the  fall-rate  for  a  given  particle  size  parameter  known, 
the  range  of  diameters  for  a  group  of  particles  can  be  determined,  depending 


moles  fp/sq  ft 
r'/hr  at  1  hr 


on  the  thickness  of  the  cloud  and  the  altitude  from  which  the  group  falls. 
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A  family  of  curves  which  shows  the  size  of  ideal,  spherical  particles  falling 
through  a  standard  atmosphere  from  various  altitudes  as  a  function  of  fall- 
rate  was  prepared  by  D.  E.  Clark''  and  appears  as  a  sketch  in  Figure  3. 

When  estimating  the  particle  size  range  at  a  downwind  distance,  it 
should  be  first  decided  whether  the  source  of  the  particle  group  originates 
from  the  stem  or  the  cloud .  The  origin  of  the  particle  group  can  be  generally 
judged  by  comparing  the  downwind  distance  with  the  downwind  range  of  the 
stem  fallout.  In  case  the  distance  falls  within  this  range,  the  particle 
group  predominantly  originates  from  the  stem;  otherwise  it  is  essentially 
derived  from  the  cloud.  Within  the  stem  range,  the  maximum  and  minimum  a 
values  may  be  calc'ilated  by  Equation  (H),  which  is  a  laborious  procedure, 

A  computer  program,  as  shown  in  Figure  7,  has  been  established  to  solve  this 
equation.  The  height  of  fall  for  the  stem  fallout  is  derived  from  Equation  (lO) 


which  is 


1.56) 


For  particle  groups  falling  inside  the  cloud  range  at  any  down¬ 
wind  location  (X,  x),  the  crosswind  distance  X  should  first  be  converted 
to  the  cloud  coordinate,  y,  according  to 


y  = 


(57) 


where  Y 

m 


a 


is  the  maxlm’um  half -width  of  the  crosswind  distance  on  a 
1  r/hr  contour,  and 
is  the  cloud  radius . 


Following  this  coordinate  transfer,  use  equations  (l5  )  (18)  ho 

ietermine  0:  and  the  two  extreme  heights  of  fail,  z  and  z  .  For  the  actual 
m  '  a  c 
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dlstance  of  fall,  the  height  of  the  cloud,  h,  should  he  added.  The  fall-rate 
is  eval-uated  from  Equation  (lO),  which  is  the  defining  equation  of  a  and 
which  yields 


V  = 

f  a 


Therefore,  with  weapon  yield  and  downwind  distance  known,  and  with 
the  help  of  the  fall-rate  curves,  the  approximate  particle  size  distrihutlon 
at  (X,  Y)  can  be  estimated.  Ttiis  information  on  particle  size  distribution 
is  usefxlL  for  studies  of  transport  phenomena  in  fallout  investigations. 

H.  Conclusion  -  The  Reality  and  Applications  of  the  O.C.D.  Model 

The  O.C.D.  fallout  model  mainly  uses  the  technique  of  the  scaling 
method,  hence  most  of  its  formulae  are  derived  from  realistic  experimental 
data.  Therefore  it  serves  as  a  very  good  bridge  between  the  ideal  theoretical 
hypothesis  and  practical  empirical  data.  Moreover,  this  model  has  been 
applied  to  several  real  test  shots,  and  was  found  to  yield  very  satisfactory 
results.  On  the  other  hand,  the  fundamental  thermodynamic  theory  on  which 
this  model  is  based  is  quite  simple  and  convincing,  and  all  the  computations 
involved  are  relatively  easy.  Another  advantage  of  this  model  is  its  wide 
application,  since  it  gives  not  only  information  on  fallout  patterns  but 
also  that  of  activity,  hence  it  is  very  valuable  to  radiological  research. 
Although  there  are  still  a  few  defects,  such  as  constant  wind  velocity, 
the  sole  dependence  on  weapon  yield,  and  the  cigar- shape  contour  assump¬ 
tions,  etc.,  this  model  is  generally  considered  very  rsalistic  and  applicable 
to  the  problems  of  fallout  research.  Its  main  applications  will  be  radio¬ 


activity  prediction,  countermeasure  design  and  biological  protection. 
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III .  COMPUTATIONAL  METHODS  OF  THE  FALLOUT  MODEL 


A.  lonlzatior  Rate  Contours 

Although  the  physical  characteristics  and  the  mathematical  approxi¬ 
mation  for  the  rising  stem  and  the  cloud  have  been  discussed  in  detail, 
little  has  been  said  about  how  the  fallout  contours  are  determined.  Even 
though  the  available  data  from  actml  land  s'orface  shots  is  small  and  limited 
primarily  to  low  weapon  yields,  some  general  conclusions  have  been  drawn 
from  the  test  data. 

Because  most  of  the  observed  properties  of  the  fallout  patterns 
are  in  terms  of  intensity,  the  fallout  contours  are  usually  given  in  units 
of  roentgens/hour  (where  the  time  was  corrected  to  a  standard  time  of 
H  +  1  hour).  The  general  shape  of  the  fallout  pattern  can  be  approximated 
by  the  overlapping  of  ellipses,  for  both  the  stem  and  cloud  fallout. 

The  dimensions  necessary  for  constructing  the  ellipses  are  the 
distances  from  ground  zero  along  the  downwind  axis  where  the  lateral  dis¬ 
placement  is  zero,  (Y  =  O).  First,  to  construct  the  intensity  vs.  distance 
curve  (see  Figure  4)  the  downwind  distances  of  interest  are  the  highest 
intensity  and  a  preselected  low  intensity,  which  varies  with  wind  velocity, 
but  for  a  15  mph  wind  is  1  r/hr  at  1  hr. 

For  the  stem  these  locations  have  been  labelled  as  follows: 
and  are  the  upwind  and  downwind  distances,  respectively,  from 

ground  zero  to  the  1  r/hr  at  1  hr  intensity. 

Xg  and  x^  are  the  distances  to  the  upwind  and  downwind  shoulders, 
respectively,  of  the  intensity  ridge. 

For  the  cloud  the  following  notation  is  used: 

X^  and  X^  are  the  arrival  location  for  the  group  of  particles  of  the 
the  selected  low  Intensity  from  the  rear  and  front  of  the 


cloud,  respectively. 


000  01 


<D 


•H 


is  the  distance  from  grovurid  zero  to  an  intermediate  point 


7 


^8 


between  X^.  and  X„. 

5  7 

X„  is  the  downwind  distance  from  ground  zero  to  the  peak  of 

the  intensity  contour. 

is  the  downwind  distance  from  ground  zero  where  the  conto^lr 
intervals  has  a  ma:x;imum  half -width. 

The  values  for  these  quantities  are  determined  by  equations 
appearing  in  the  Ionization  or  Intensity  rate  computer  program.  Because 
of  the  limited  characters  available  on  the  computer  printer,  the  notation, 
used  in  the  analysis  of  the  model  had  to  be  changed.  The  corresponding 
notation  between  the  model  and  the  computer  program  is  given  in  Table  III. 

The  following  discussion  describes  the  method  for  calc\ilating  the 
X's  and  how  the  scaling  functions  can  be  found  in  the  computer  program  in 
Figure  5-  First,  yield-dependent  scaling  functions  for  the  size  parameters 
a  =  0^  a^,  a^,  ,  .  .  >  have  been  determined  from  data  that  has  an 

effective,  or  average  wind  speed  of  15  mph.  These  scaling  functions  are 


then  multiplied  by 


W 

- r-  to  correct  for  the  wind  velocity  V,,.  ■^j'the 

15  mph  W  2,3 


ionization  rate  between  X^  and  Xy  is  calculated  by  two  formulas .  For  yields 
W  :>  lOOKT,  ^  3^^)  ^f^2  3  vhere  the  values  for  ^(l) 

and  n  were  generated  by  fitting  curves  through  the  values  of  a  Table  in 
Chapter  Three  of  Carl  F.  Miller's  manuscript  "Fallout  and  Radiological 
Countermeasures" .  For  W  <  100  a  curve  was  fitted  through  the  first  three 
values  of  ^  given  in  another  table  from  Miller  and  adjusted  for  wind 

‘-}J  y 

/  W\-0.79 

velocity  by  multiplying  by  ■  The  resulting  equations  appear  on 

lines  O7A1  through  0802  in  the  computer  program  appearing  in  Figure  5- 


In  all  the  equations  that  follow  the  scaling  function  formulas 
and  constants  are  given  in  Table  IV.  For  any  scaling  function  dependent 


o^ou 
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0400 

0400 

0400 

0400 

0404 

0408 
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0422 

0436 
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0484 
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0621 

0628 
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Pigore  5*  Computer  Program  for  Ionization  Rate  Contours 

BkC-2Z0  STANDARD  VERSION  3/1/62 

DUMP  J  S  T6S 

eOMHENT  RICHARD  E.  UAMES  POH  A  P  633.  S 

integer  F  .  G  .  J  >  P  *  N  >EDC  .  799...  i 

array  ALP(9)>I(9).IC(49I<KAA<9).L(49>.LI(9>,4(49)>PHI<9>,X(9>  S 

PROCEDURE  SEMHXl#  I1*X2. 12  %  H,B)  . . S 

SEGIN  C  >  LOG( 111  i 

D  >  L0Q(I2)  1 

H  ■  4X2  -  XI)  /  (D  -  C)  S 

B  ■  (XI  .  0  -  X2  .  C»  /  <D  -  C)  S 

RETURN  END  SEHIO  $ 

FUNCTION  LPF(N.W)  •  2.3025851  t(N>  4.  LM  .  H(N>  S 

FUNCTION  Fr<N.H)  •  EXP  <  I.FF  ( N»  N  )  )  % 

FUNCTION  IFF(X,MM»8B)  •  EXP((X-8e)/HM)  $ 

FUNCTION  SLL<1*M.8»  «  M.LOOdl  ♦  8  S 

READdSQAZ)  S  INPUT  UAZ(E0C>F0R  UXl.l.EDC)  S  (F  .K  F  )  >  M<F  k )  )S 
START..  READdAEBWV)  S  INPUT  E8HV(H,VU>  i 

tW  «  lOG(W)  » 

F  ■  (  W  QTR  28  )  1 

Q  XH  UTR  9.0)  A 

V  ■  VM/15.IJ  S 

AH  ■  FF(39  ♦  F  ,  H)  S 

A  a  (2.45a«3)Ha0.43l  % 

B  a  (l,40»*3)Ha0,30'’  $ 

LA  a  LFF(3>M)  S 

LARS  a  LFF(2*H)  1 

Las  a  LFFIB.W)  | 

HS  a  FF(4>H)  S 

H  a  FF(’54F.U)  a 

KA  ■  LARS  /  (N  -  RS)  t 

LAO  a  LA  -  KA.H  S 
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Figure  5  (Continued).  Computer  Program  for  Ionization  Rate  Contours 


0632 
0637 
064!> 
0651 
0669 
0683 
0685 
0709 
0713 
0717 
0721 
0733 
0741 
0741 
0749 
0754 
0761 
0761 
0762 
0767 
0776 
0781 
0786 
0793 
0  796 
0802 
0802 
0811 
0819 


2S  ■  (LAS  -  LAO)  /  KA  S 

ALP23  »  V.KFt7,W)  S 

ltl>  ■  1(4)  c  {(9)  c  1.0/V  % 

FOR  J  «  4.8«6.7«6  i 

ALP(J)  «  V.FF<21*J>M>  i 

ALPS  ■  ALP(5)  S 

20  ■  Q  (  1900.0  *  (ALP23  «  0 . 020) 2S ) 7 A.P23  *  (l-6){H-6)  t 

AS  c  exP(LAS)  $ 

X(2>  •  ALP23.26  'AS  % 

X(3)  ■  ALP23.2S  «  AS  S 

XC4)  ■  ALP«4)  1  ALP<4)  .  20  -  19U0.0  )  7  «  ALP(4)  +  0,020  >  i 
K12  ■  FF<41  F,W)  I 

either  if  W  LSS  100  S 

BEUIN 

1023  ■  EXP(t-0.20i7e297LH  ♦  1. 5887143). H  ♦  7.4558767)  S 

123  •  1023  .  VA-0.79  I 

ENLi  1 

OTHERWISE  S  BEGIN 

V'^’N  •  -EXP(  U  (  (-7.94235B6«*-4)LW  ♦  0 . 071809740  )LH 

-  0.22171897)LH  ♦  0 . 97745048  )  LH  -  1.82115694)  $ 


VHSN  ■  VH»VUN 

1023  ■  EXP<  (((-0.0017555634LH  *  O . 046499101 ) LH  - 
0.427651U4>LH  ♦  1.05351e7)LW  ♦  13.053662  ) 
123  ■  <023.VWSN 
1 023  ■  1023  15»VHN 

END 

X«l)  •  K(2)  -  L0B(123)7K12 
ALP5P  •  V.FFdl.W) 

ALP<9)  ■  V,FF(23»F,W) 


0828 


either  if  ALP(5)  oeq  ah 


-k6- 


Pigure  5  (continued).  Computer  Program  for  Ionization  Rate  Contours 


0637 

0841 

0845 

0849 

0851 

0687 

0883 

0890 

0897 

0907 

0907 

0918 

0918 

0923 

0937 

0951 

0961 

0974 

0974 

0981 

0981 

0984 

0995 

1021 

1029 

1036 

1038 

1054 

1063 

1068 


KAA(5»  ■  rF<30*F*H>  * 

OTHERWISE  _  ,S 

BEGIN  KAA(5)  »  Pr<32,H)  ,  « 

ALP<5)  «  AI.P5P  END  I 

FOR  J  ■  6.7  . . .  . .  . . 

KAA(J>  •  FF<22<-2J*<ALP<J)  LSS  AH>  ,  W)  ^  .  * 

KAA«9»  ■  77138. HI  I 

ABSQ  ■  FF(l.H)  S 

FOR  J  •  6.7.5  * 

BEGIN  either  if  aLPIJ)  LSS  AH  .  S 

OC  ■  ALP23  * 

OTHERWISE  * 

00  ■  ALPIJ)  -  ah  S 

PHKJI  ■  <(AUP(J)  ♦  AH)  ♦  80RT(ABS8  ♦  (  ALP<  J)  +  AH ) .  ( AtP<  J  I  +  AH)  )  > 
/((QQ)  -»  SGRTIABSO  *  QU.UO))  END  S 

FOR  J  ■  6.7.5  $ 

KJ)  •  2.A.KAA(J).L0QIPHI(J))  S 

EITHER  IF  ALP(9)  GEO  AH  S 

AUP<5)  ■  ALP5P  S 

OTHERWISE  » 

ALP«»)  ■  ALP5  S 

FOR  W  •  (8.1.9J  S 

X(J)  ■  <6600  *  1020UFmW*<0.445  -  U  .  2BlF  ) )  ALPIJ )  S 

IF  ALP(5)  GTR  ah  t 

X<5)  ■  X<9l  -  8  .  SORHABSO  ♦  ALP5P  .  ALP5P  )  i 

ALP(5)  ■  ALPS  S 

SEKI(^(7I.1<7).X(9).1(9)  S  8789, 0769)  t 

1(8)  ■  IFFIXIO) ,N789,B789)  S 

Y815  ■  GXP(((((  9.9684B14**-'<)LW  -  0 . 027025999 1 L  W  * 

0.22433092)LH  -  0 . 12450 012 > L H  *  6.7992249)  6 


Figure  5  (Continued).  Computer  Program  for  Ionization  Bate  Contours 


1075 

1089 

_109g 

1095 

1098 

1109 

1117 

ll!5 

1162 

1192 

1192 

1213 

1216 

1220 

1233 

1236 

1245 

1251 

1267 

1283 

1299 

1315 

1331 

1347 

1347 

1347 

1377 

1377 

1377 


SVW  •  <(3.U892VU.VU/(9»«5)  •  3.6099299/9000  *  0.0174592)VW 


•  0.346408)99  *  3.44944  S 

va  ■  7615  .  SV9  S 

112)  >  !<3i  «  123  % 

FOR  J  »  (1.1.9)  S 

LKJ)  •  L0<a(nJ>)  % 

WRITgdSECl.ECFl)  $ 

OUTPUT  ECKFOR  J  ■  (1.1.9)  t  J  .  FOR  J  a  (1,1, 9>  S  X(J)/52Bg  , 

FOR  J  a  (1.1.9)  5  MU)  ,  W  ,  VR  ,  76/5280  )  S 

FORMAT  ECFl(Il2,8113.M3.aX  •. 9S13. 6, M4 , . (  « , 9813. 6. W6, X3S • 0 , •  KT*, 

XIO.O.*  HPHa.645.*y6  a».X6.4.U6)  S 


COMMENT  computation  CF  VALUES  RHICH  WILL  BE  OONSTANT  INNEXT  SECTIONS 

YOS  a  FFiao.W)  _ _  _  I 

YS  »  YOS  .  LOQ(123)/LOStl_02J)^ _ _ _  S 

Xe7  a  X(8>  “  1(7)  _  S 

0  a  S0RT(Y8.Y8  ♦  X87.X87)  $ 

LI7  a  L06  (  M7)  /  1(9)  )  % 

SEMI(X(1).1.0.X(2).I23  S  M12.B12>  S 

SEMI(X(3).123.X(4).I(4>  S  M34.B34)  S 

$EMl(X(5).I(9).X(6>.M6)  S  M56.B56)  I 

SEMI(X(6).I(6).X(7>.1«7>  S  M67,B67)  S 

SEHKO. 0.123. YS. 1.0  S  MYS.BYS)  % 

SEHI(0.0,M7).D.1.0  t  HYB.SYB)  S 

MRITEdSECFS)  i  FORMAT  E0F5(*  l».BB.»Xl  •  X2  X3  P  X4*. 

B4,*X5  t  X6  -  K7  t  X6  -  X9* , b5. *Y-ST6M*, 65 . *Y-CLOUD*. B9 , 
•C9*.Bt5,«F0CIa.N4)  S 

COMMENT  ME  NON  FIND  THE  ENDPOINTS  OF  THE  SEVERAL  ELLIPSES  AND 

THE  FOCI  OF  THE  UPHIND  CLOUD  ELLIPSES  S 

READ(SSEAIC)  S  input  EAICIEOC.FOR  Ja(l.l.EDJ)  S  IC(U>)  S 
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Figure  5  (Continued).  Computer  Progi'am  for  Ionization  Rate  Contours 


1416 

begin 

IJ  ■  !C(Ui 

S 

3  -tSO 

XCl  ■  SULllJ  •  M12  .  S12> 

s 

1489 

XC3  •  SLL( lU  •  H34  ,  834) 

t 

1438 

XC4  •  SLUt 1 J  «  M56  t  856) 

% 

1447 

IF  IJ  QTR  1(6) 

i 

1455 

XC4  ■  SLL< IJ>M67.B67) 

i 

1461 

XC6  ■  SLLCiU  •  M789  ,  8789) 

% 

1470 

Ties  ■  SLLdJ  •  MYS  .  BYS) 

6 

1479 

DPR  ■  SLL( IJ  >  Mye  >  BY8) 

S 

1488 

YlCC  «  DPR  .  Y£  7  D 

i 

1493 

XPRS  ■  X(2)  •  LOG  (  183  /  iJ  )  /  K12 

.6 

1504 

XXa  X(3)  *  (X(4)  .  X(3))LO&(123/  IJ  >/L0S(  1 83/ 1 ( 4 > ) 

% 

1585 

SA  •  3.14l5927<<X<2)*XPfiS)<X<8>-X'»RS)»YS(XX4XPRS-8X<8))/8) 

% 

1542 

XXYY  ■  L00(I<7)/|J)  /  LI7 

% 

1590 

YY  ■  Y0.XXYV 

1 

1593 

XXO  ■X(7)  ♦  <X(9)  -  X(7))XXYY 

S 

1598 

P  ■  6  -  (  U  L&J  1(6)) 

S 

1570 

XPRC  ■  X(PI  ♦  (X  P+1)  -  X(P))LDG(  IJ  /  -Ot  :(P+1)/1(P)) 

t 

1594 

CA  »  1.5708  YY  (  XXC  -  XPRC  > 

s 

1599 

CB  •  (X(7>  ♦  X07,Y1CC  /  Y8  ) 

I 

1605 

ELLC  ■  SQRT((CB  -  XC4)(C9  -  xC4)  -  YlC:.YlCC) 

s 

1618 

WRITE(SSEC2.ECF8) 

s 

1686 

OUTPUT 

EC8I lJiXCl/52e0,xr3/5280.Xg47528j,XC6/?280, YICS/6P80. 

1696 

Y I CC/52eu,Ce/5260.(C8-ELLC) 75280 1 (CB+ELLC) 75280  ) 

% 

1680 

format 

ECF2(X9.2« 2 02. 2X11. 4), 63.2X11.4,83.3X11.4.94) 

s 

1691 

END 

s 

1692 

Z99W1NDV6L0CITYMPM  •  VH  S  Z99YIELUKT  W 

% 

1700 

COMMENT 

here  wi  find  the  constants  for  a  fdrhul*  for  finding 

1700 

INTENSITY  IN  The  DOmnHino  CLOUd  REGION 

s 

1704 

X?  •  X(7)  /  5280  5  XB  •  X(B)  /  5280  S  X9  >  X(9)  /  5280 

i 

Figure  5  (Continued)-  Computer  Program  for  Ionization  Hate  Contours 


1712 

1718 

1722 

1729 

1733' 

1741 

1770 

1770 

1770 

1770 

1611 

leie 

1619 


787  ■  X6  •  X7  S  796  •  79  •  XS 

HYB  •  re  /  9260 

744  ■  798  .  X96  -  787  .  767 

Cl  «  L17  /  744 

WRIT£(SSEC6.EC^6) 

OUTPUT  EC6(LI7«>L0a<I<9})<X87.Cl<X9B.Cl>X44/(Hre.Mre).X7 
FORHaT  6CF6(W,*F0H  OOHNKINO  CLOUD  INTENSITIES.  *.W4.*|  ■  gXP  (*, 

910. 3»*  ♦  (•.X12.10»*J .7X7  -  ( *. Xll. 7, * ) .SORT < XX7 . XX7  ♦  <- 
S10.8.*)Y.r  )  )««W4  «  B5;<»WNERE  7X7  ■  7  •  •«S1Q.8.B5. 

*.  (  7  4f.O  Y  IN  HILES 

IF  U  EQL  10000.0  S  STOP  4464$457 
60  START  % 


FINISH 

compiled  program  ends  AT  1820 
PRQSRAM  variables  BEGIN  AT  3739 


Figiire  5  (Continued).  Computer  Program  for  Ionization  Rate  Contours 


EXP  «  9.111X921  ♦  I  . 001X390636). XX7  -  (  .  020635561)  .S0RTUX7.XX7  ♦  <  X5.i7X538)Y 
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Pigyre  5  (Continued).  Computer  Program  for  Ionization  Rate  Contours 


LMT  UAMi.  rtttSD  mm 

U  PftOfiftAH 

_«T*nr 

VAHUILE  IN  PMOHAM 


atari  uoaii 

NUHRCH  OtJlHES  EXSCUT£D 
00.01 

VAI.UE 


A  _ 

At«Q 

AH 

ALf93 

A(.P9 

4LP9P 

aA 

i 

BYt 

BTO 

.%M 

•94 

•94 

•97 

t7BB 

CA 

CA 

Cl 

0 

OPR 

9&C 

ELLC 

F 

0 

H 

N 

[029 

U9 

KA 

R12 

LA 

LARI 

LAI 

LAO 

Lit 

LH 

NYI 

HY« 

Nil 

NS4 

Klft 

NA7 

NtW* 

P 

QQ 

RB 

•A 

lYW 

V 

VM 

VHN 

YHIH 

B 

XCl 

XC9 

IC4 

ACA 

RR«C 

APR! 

IX 

xxe 

XKVV 

K44 

X7 

Al 

1B7 

t% 

X9« 

TlCC 

_  ..tlCI .  - 

T» 

tt 

Yl 

'  Y^19 

-  -  IK  -  -  -  . 

It 

mnwYtL. 

IPtTliLDRT 


,12079794.  06  .. 
,)4107990.  08 
a70219B4*  01 
.62979402.  00 
,91991296.  00 
.21199911*  PI 
.1B709B19.  09 
.22198903,  09 
.9690;292*  09 
.37991424,  09 
•>. 60421978,  09 
.19264999,  99 
.14990999,  09 
29397919*  07 
.27991092.  07 
.92707093*  11 
,92726409*  06 
.47416109,-04 

.  ,97091449,  06 

.46097390*  09 
OO0OOOUOO9 
.21919324,  06 
OO0OUOOOO1 

oooouooooi 
.76807999*  09 
.90000000*  04 
.99499941*  04 
,39499941*  04 
0000001)009 

.46930770.-04 

.1‘271979*-03 
.11771117,  02 
.51693793*  HI 
.90366434,  Dl 
,91969230*  01 
.91114931*  01 
.94109404,  0i 
-.•1394(197,  04 
.69219246*  02 
,98719999*  04 
-.14943980,  09 
.11616057,  09 
.39999764*  06 
-.27221409*  09 
DOODOOOOD9 
.4t620?70*  00 
.4<771329*  04 
.40229969*  10 
.lOOOUCDOi  01 
.lODOUOOO,  01 
•19OQ000O*  02 
-.73994731*  00 
.13492961*  00 
.10000000*  09 
.169T3960*  04 
.43000790,  09 
.l(i70«509*  06 
,60i99«20*  09 
.19709309*  09 
,19073919*  04 
,43004799*  01 
.60999921*  06 
,t2Ui9U«i  06 
,19819216*  66 
<§•300099,  02 
.•B94«S01i  02 
.30349402.  02 
,92»69o24.  os 
^43940094,  03 
,41796029,  09 

_ a39..0X7ta*.ji- 

.69927293,  09 
.  vUtKiUt  -•*  - 
.96917299,  09 
,IAA99761«  09 
.'34439768*  89 
.,44iimujLi- 


.39312209,  09 

itmmu— 

0000010096 


52 


Figlire  5  (Continued).  Computer  Progi*am  for  loniKation  Pa.te  Contours 

ARRAY  ALR 


tOOOOOOOO* 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•42299999 • 

01 

•  al49C23«i  00 

•14927949. 

01 

•40497992 • 

01 

•9t9«79a4» 

01 

•99449790i 

02 

ARRAY  t 

alOOOOOOO* 

01 

•99496941 • 

C4 

•99499941* 

04 

•  lOOOOOOO* 

01 

i21972«72*  01 

•92472Z92I 

04 

•909a9020* 

04 

aoeaosont 

04 

•looooeoot 

01 

ARRAY  tC 

aXOOOOOOOa 

01 

• 30000000 i 

01 

•  lOOOOOOOi 

02 

•lOOOOOOOi 

02 

•  lOOOOOOOi  01 

•  lOOOOOOOi 

01 

•  lOOOOOOOi 

04 

•iOOOOOOOt 

04 

•lOOOOOOOi 

09 

•200000001 

09 

•  lOOOOOOOi 

03 

•  OOOOOOOOi  03 

•lOOOOOOOi 

04 

•20000000. 

04 

•loaooooot 

04 

•49000000 t 

04 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi  00 

•  OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•  OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi  00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•00000000. 

00 

•00000000 « 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi  00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

tOOOOOOOO. 

00 

• 00000000 « 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi  00 

•00000000# 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi  00 

i OOOOOOOOi 

00 

•  OOOOOOOOi 

00 

ARRAY  KAA 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

oo 

•OOOOOOOOi 

00 

Il9929aaii>i04 

•49940727 i-Rl 

•72949794 i-Ol 

•OOOOOOOOi 

00 

•43299970 i-OA 

ARRAY  b 

•  MOOOOCCt 

00 

•107900001 

01 

•39990000a 

01 

•23130000# 

01 

•992000001  01 

•42290000. 

01 

••90900006. 

00 

iaoaoooooi 

01 

•990900001 

01 

•99940000* 

01 

••94000000i^01 

•i99ooeeoi  01 

•42990000. 

01 

•99920000* 

01 

i42990000i 

01 

•400SOOOO • 

01 

•44100000> 

01 

•31900000# 

01 

•92v20000i  01 

•92290000I 

01 

•19440000I 

01 

•40490000i 

01 

•  137100001 

01 

•99000000* 

00 

•27000000* 

00 

•i 17900000*  00 

•lOOOOOOOi-Ol 

•4iooooeo*-oi 

iiaooooooi 

00 

>• 929900001 

01 

••29990000* 

01 

-iiiaiooooi 

01 

•OOOOOOOOi  DO 

••12940000* 

01 

••12290000I 

01 

-.yaooooooi 

00 

•• 10770000 i 

01 

-•21960000* 

01 

••491000001 

00 

•iB3700000i  00 

-•zooiecooi 

01 

••290000001 

01 

•OOOOOOOOi 

00_ 

•oocooooo* 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

iOQOOOOOOi  00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

ARRAY  M 

•39200000« 

M 

•  Y9000000i>C)l 

•43100000# 

00 

•lllOOOOOl 

00 

•449000001  00 

•19400000* 

00 

•790e0000i>01 

•34900000. 

00 

•0*00000  « 

00 

•919000001 

00 

•99000000i-01 

•AOlOOOOOi  00 

•  20000000* 

00 

•90900000* 

00 

•90900000. 

90 

• 99900000 i 

00 

.«3130Q000i 

00 

•31900000i 

00 

illlOOOOOi  00 

•40000000. 

00 

•49700000. 

00 

•laoopooo* 

00 

»•  12*00000 • 

00 

#14400000# 

00 

•99000000 i«01 

•  22000000f01 

•19 000000 ••01 

•14100000. 

00 

•19100000. 

00 

••27900000. 

00 

•<»!.'7200000# 

00 

-•409000001 

00 

•  OOOOOOOOi  00 

-•T4000000.-01 

-•22000000*-01 

>.ITOOOOOO.«01 

-•20000000t« 

•01 

••912000001 

00 

-•140000001' 

•oi 

•29700000i  00 

-•40400000. 

00 

-•9170C000i 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi  00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

ARRAY  PHI 

•OOOOOOOO* 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•OOOOOOOOi 

00 

•172097991  01 

•19199041i 

01 

01 

•OOOOOOOOi 

00 

•eooooooot 

00 

ARRAY  X 


••A»AZ1»T9(  09  t3ea9919e«  OA  •4090199**  09  *I»2*A999»  06  tZaiTZlIOt  09  •Il99llf2t  09  tSOTOaOTO*  09 

*«9l06t91«  09  .2799 1091 t  07 
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Table 

III 

LIST  OF 

SYMBOLS 

Notation  of 

Notation  due 

Notation  of 

Notation  due 

computer 

primarily  to 

computer 

primarily  to 

program 

Miller 

program 

Miller 

A 

a 

ELLC 

Distance  from  center  to 

focus  of  ellipse  of  chosen 

ABSQ 

a^/b^ 

a/h 

contour  in  upwind  cloud 
region 

AH 

P 

F=1  if  W  >  28  ET 

ALP  (  ) 

ag,  a^,  o^, 

F=0  if  W  <:  28  KT 
(used  to  differentiate 

ALF23 

®2,3 

scaling  functions) 

ALP5 

ac 

G 

G=1  if  W  <  9  KH? 

p 

G=0  if  W  ^  9  KT 

ALP5P 

a' 

p 

H 

h 

AB 

a 

B 

s 

b 

I  (  ) 

Ipj  \>  lyj  ^8’  ^9 

ic  (  ) 

Intensities  of  chosen 

MfS,  BYS 

* 

(0,  I^),  (D,  1) 

contours 

MY8,  BY8 

* 

(0,  12^3),  (Y^,  1) 

TJ 

Intensity  of  chosen 
contour  l9  , 

M12,  B12 

M3i^,  B34 

* 

* 

O'!-  V’ 

(X^,  12^3),  Ii^) 

1023 

123 

2,3 

I 

M56,  B56 

* 

(X5.  15)^  (X6>  V 

J 

2,3 

Dummy  variable  used  as 

m67,  b67 

* 

(Xg,  Ig),  (X^,  I^) 

"i"  in  etc. 

M789,  B789»^ 

(X^,  I^),  (X^,  Ig) 

KA 

CA 

CA(I) 

KAA  (  ) 

Ki(l)A^,  K'(1)S„ 

CB 

X 

K12 

K  ,  2 

c 

1 

Cl 

LI7  X4U 

L  (  ) 

Array  containing  constants 
for  scaling  functions 

D 

D 

LA 

In  a 

DPR 

D' 

LARS 

In  a/R 

EDC 

Number  of  contours 

s 

studied 

LAS 

In  a 

54^ 


Talsle  III  (Continued) 
LIST  OP  SYMBOLS 


Notation  of 

Notation  due 

Notation  of 

Notation  due 

computer 

primarily  to 

computer 

primarily  to 

program 

Miller 

Miller 

lAO 

In  a 

XG6 

7,  9  ** 

0 

LI  (  ) 

In  1° 

XPRC 

X' 

LI? 

InCl?/!^) 

XPRS 

X' 

LW 

In  W 

XX 

X 

M  (  ) 

Array  containing 

XXC 

X 

constants  for 

•  T  /t. 

P 

scaling  functions 

Used  in  equations 

XXXY 

^  (  T  ) 
in  {W^9) 

involving  x' 

X44 

X98^  -  X87^ 

PHI 

a 

X7 

X?  in  feet 

QQ 

"2,3  -h 

X8 

,|gn  feet 

RS 

R 

s 

X87 

Xg  -X? 

SA 

SA(I) 

X9 

X^  in  feet 

SVW 

S(v  ) 

'  w 

X98 

X^  -  Xg  in  feet 

V 

VW 

V  /15 

w' 

YICC 

Maximum  cloud  width  at 
chosen  contour 

V 

VWN 

w 

YICS 

Maximum  stem  width  at 

n 

chosen  contour 

WSN 

n 

V 

w 

YS 

Y 

B 

W 

W 

YY 

Y 

X  (  ) 

x^. 

YOS 

Y° 

Xg,  X?,  Xg,  X^ 

Y8 

s 

XCl 

1,  2  ** 

Y815 

Yet  15) 

XC3 

3,  4  ** 

ZS 

z 

s 

XC4 

5,  7 

ZO 

Z 

0 
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Talale  III  (continued) 


LIST  OF  SYMBOLS 
Continued 


Notation  of 
computer 
■program 


Notation  due 
primarily  to 
Miller 


Notation  of 
computer 
program 


Z99WINDVELOCITirMPH  Z99YIELDKr 


Notation  due 
primarily  to 
Miller 


W 


The  slope  and  intercept^  respectively ^  for  a  semi -logarithmic  line 
through  the  two  points  whose  coordinates  are  given. 

**  The  downwind  distance  (for  the  chosen  intensity)  to  the  fallout  pattern 
between  the  points  Indicated. 


Data  Table  for  I' 


for 

N 

...  .  1(N) . . 

..  m(n) . 

uhf- 

1 

0.486 

0/262 

In  k  /'ir.^ 

2 

1.070 

0.098 

in  A 

3 

.3.389 

0.431 

E 

Q 

4 

2-319 

0.333 

ti. 

5 

3.820 

0.445 

6 

4.226 

0.164 

“23 

7 

-0.509 

0.076 

Ir* 

8 

2.880 

0-348 

>'l 

9 

3.308 

0.496 

10 

3.564 

0.319 

a,! 

b 

11 

-0.054 

0.095 

12 

3.850 

0.481 

13 

4.255 

0.200 

X,^ 

14 

3.862 

0.586 

15 

4.265 

0.305 

Y 

l6 

4.005 

0 . 596 

17 

4.4i0 

0.315 

V 

18 

5.190 

0.319 

^9 

19 

5.202 

0.311 

..  .O 

20 

3.223 

0.400 

Y 

21 

3.644 

0.46? 

A— 

y 

22 

4.049 

0.186 

rtf’ 

I  Rate  Contour  Program 


^9 

23 

1.371 

-0,124 

24 

0.980 

0.146 

25 

0.270 

0.089 

“5 

26 

-0.176 

0.022 

«6 

27 

0.030 

0.036 

“7 

28 

0.043 

0.l4l 

«8 

29 

0.185 

0.151 

V  fi 

30 

-3.286 

-0.298 

^5^a 

31 

-2.889 

-0 . 572 

32 

33 

-3.185 

-o.4o6 

% 

34 

-1.134 

-0.074 

35 

-1.225 

-0.022 

'<  s 

36 

-0.989 

-0.037 

Va 

37 

-1.079 

-0.020 

38 

-2.166 

-0.552 

a/h 

39 

-0.431 

-0.0l4 

40 

-0.837 

0.267 

'vf 

41 

-2.503 

-0.4o4 

^12 

42 

-2 . 600 

-.337 
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variatle,  for  example  H  (height  of  cloud),  the  program  determines  which 
formula  should  he  used.  The  equation  for  H  appears  on  line  0613  of  the 
contour  program. 

The  following  chain  is  then  started 


H  =  PF  (5  +  F,  W) 


line  0613 


PP(N,  W)  =  EXP(LFF'(N,  W)) 


line  0U54 


IPF(N,  W)  «  2.3025851  L(I!I)  +  LW  •  M(N)  line  046t 

The  program  then  finds  for  W  <  28  KT  that  F  =  0  since 

W  greater  28  is  false  and  F  is  a  Boolean  variable.  Therefore  for  II  •=  5, 
one  finds  in  the  Table  IV  that  L(5)  =  3 '820  and  M(5)  =  0.445.  The 
complete  equation  of  H  for  W  <  28  KT  is 


H  *  exp[(2.3O2585l)(3.820)  +  log^W  •  0.445]  ,  (59) 


log^H  =  (2. 3025851) (3 -820)  *  log^W  •  0.445  , 


I  59a) 


or 


log^^H  =  3.820  +  0.445  logj_QW. 


(59b) 


The  variation  of  and  with  wind  speed  depends  on  the  wind  modified 
parameter  a  (lines  0?13,  0717)*  The  equations  for  X..  and  X,  are  then 

<-  3  J  X  H 


_  g-303  log^Ig^3 


K 


1,2 


Ww  -  ^^’36) 

(1  +  9.273  X  10"\v^) 


(60) 


f6l) 
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or 

-  1900) 

^4  “  +  0.020 


(6la) 


for  15  mph  wind  velocity. 

The  equations  for  the  distances  Xg  through  (line  0995) 

given  hy 


For  i 
X. 


X. 

1 


6,  7^  8  9 

6.6o  X  ioV’^^5  a^,  W 

1.68  X  10 a. ,  w 

i’ 


=  1  to  28  KT  (62) 

=  28  to  10^  KT  {62a'l 


The  above  eq-oaticns  also  hold  for  X^  when  is  greater  than  a/h.  When 

a  is  less  than  or  equal  to  a/h,  then  X_  is  determined  by 
J  5 

X  =  6.6o  X  lO^^’^^^a'  -1.4o  X  10^W°'^°°  V  3.06w°’^^  +  a'^  j  W  »  1  to  28  KT  (63) 

5  5  5  ■ 

cr 

X^  =  1,68  X  -  1.40  X  10^° J  3.06W°’^^  +  a'^  ;  W=28  to  10^  KT  (63a) 

^  y  y 


w 


L^i'ext  the  intensities  1, ,  I,, ,  and  aie  set  equal  to  - r-  . 

1'  4’  9  ^  15  mph 

The  equations  foi'  Ij,,  and  are  on  lines  0951  through  0961.  Ig  is 


7 


determined  from  a  semi -logarithmic  line  between  (X  ,  1  )  and  (X  ,  I  ) 

17  9  9 

with  the  value  ot  Xg  as  previously  calculated. 

After  X  through  X  and  their  corresponding  intensities  have  been 

i  9 

found,  a  graph  can  then  be  drawn  of  the  intensity  profiles  or  elevations. 

A  plot  of  Log,^l  (1)  against  distance  X  for  a  10  MT  weapon  is  given  in 
f'igure  4.  For  any  set  of  intensities  the  intersection  (downwind  distance  X) 
of  the  contour  ellipses  with  y  =  0  can  be  read  off  the  graph. 


-f??- 

In  the  computer  program  semilog  equations  vere  found  for  the 
Intensity  vs.  Distance  curve  so  that  any  pre-selected  intensity  contour 
intervals  were  included  in  the  computer  printout. 

Nextj  to  find  the  lateral  locations  or  width  of  the  contour 
patterns  for  the  stem,  the  half  width  Yg  for  the  stem  fallout  is  needed. 

Y„  is  defined  as  the  lateral  distance  from  the  center-line  of  the  stem 

O 

pattern  to  the  Ir/hr  at  1  hr  (for  15  mph  wind)  contour  (lines  1216-1220). 

A  semilogarithmic  line  from  (l„  ,  Y  =  O)  to  (l,  Y„)  is  then  calculated 

c,  J  D 

on  line  1315*  Using  this  equation  the  lateral  displacement  of  any  intensity 
contour  of  the  stem  can  he  found. 

The  maximum  pattern  half  width  for  the  cloud,  Yg,  is  the  maximum 
distance  from  Y  =  0  to  the  Ir/hr  at  1  hr  contour  when  the  wind  velocity 
is  15  mph.  The  calciilation  of  Yg  with  its  variation  due  to  wind  speed 
occurs  on  lines  105^  to  1092 . 

To  find  the  lateral  displacements  of  the  contour  locations 
between  (x  ,  O)  aind  (Xg,  Yg)  a  semilogarithic  line  from  {ij,  Y  =  O)  to 
(l,  D  =  V  Yg^  +  (Xg  -  X.^)^)  is  calculated.  In  lines  1236  through  13^7 
the  X  +  Y  and  foci  are  then  found  for  the  ellipses.  Using  the  information 
in  the  computer  printout  an  ionization  rate  contour  map  can  he  drawn. 

(Figure  6) , 

B.  Fractionation  Numbers 
1.  First  period  of  condensation 

For  a  ground  surface  burst,  the  ratio  of  the  total  moles  of 
liquid  carrier  to  the  molar  volume  n(i)/V  is  computed  hy  Equations  (32)  and 
(33)  with  W  =  10  ^E^  or  10,000  KT,  and  T  =  l673°K.  They  yield 
t  =  53  sec 

[n(-«)/v]Kr  =  3.25  X  10”^  atm. 


Figure  6.  Selected  Ionization  Bate  Contours 


. . 

r 
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For  a  typical  element  Sr-89,  the  procedure  of  finding  its  fractionation 
number  of  first  period  of  condensation  r^(89)  is  given  below; 

Decay  chain  of  Sr-89; 


Here  a  neutron-emitter  is  involved.  Eut  the  Eolles  and  Ballou 
data  have  already  included  this  effect,  hence  no  alteration  needs  to  be  put 
on  the  number  of  atoms  N(A,  t,)  in  their  work.  Moreover,  the  values  eslc’uloted 
by  them  according  to  Glendenin  theory  are  generally  considered  to  be  more 
reliable  and  also  used  in  the  following  calculations . 


Se 

Br 

Kr 

Rb 

Sr 

k.  =  p_ 

1.01  X  10^ 

l.Olr 

0.605 

4.61  X  10‘^2 

J  J 

...0  .1 

j  "  TwITJ/vTot 

3. 108  X  lo'^ 

32.00 

18.62 

1.4l8  X  lo""" 

tj 

0.21 

309 

9t 

3 

r^'(A,  t) 

- 

0.21 

309.21 

(i03.21 

406.21 

:'ivA,  t) 

1  +  K. 

- 

6.363  X  10’^ 

0 

^ .  791 

3 

0 

-  N(A,  t) 

6 . 363  X  10  ^ 

(2)  6363 

I+.79? 

7-797 

j 

(89) 

- 

C 

0 

d 

0 

0.0119 

0 . 01919 

?;ote:  data  unlisted,  usually  negligible 

The  r^\S9)  for  Br,  Kr,  Fb  are  the  by-products. 

.  Second  period  of  condensatjon 

The  fractionation  number  of  second  period  of  condensation  r’  '.A,  t’; 

o 

depends  very  much  on  the  time  and  the  elements  involved  in  the  decay  chain. 
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When  the  time  Is  small,  the  sublimation  pressures  of  most 

P.V 

elements  are  still  quite  significant,  and  the  values  of  are  usually 
greater  than  those  of  0.23  BWY?.  In  the  second  period  of  condensation 

J 

a  portion  of  each  such  element  remains  in  the  vapor  phase.  Therefore, 
in  computing,  the  ratio  of  these  two  values  should  be  set  equal  to  unity. 
As  time  increases,  the  sublimation  pressures  of  all  elements  except  the 
rare  gases  and  As  and  Se,  decrease  sharply  and  soon  become  negligible. 
Therefore,  the  second  period  fractionation  number  for  a  mass  chain  not 
involving  the  rare  gases  and  As  and  Se  will  be  greatly  simplified: 


r'  (A,  t) 
o  '  ^ 


1  -  r^(A,  t). 


(6k) 


Four  typical  examples  are  given  below: 

(a)  10  MT  ground  surface  burst,  60  sec  after  burst,  Kr  (83)* 

fiecay  chain: 

,Se, 


Ge  -►“As! 


‘Se, 


Br->“Kr 


This  involves  a  split 

chain,  but  the 

Bolles  and 

Ballou 

data  also 

have 

taken  this 

effect  into 

account . 

Ge 

As 

®«1 

Br 

Kr 

ro(83) 

— 

0 

0 

0 

0.0066 

0.0066 

1-r 

o 

— 

1 

1 

1 

0.993»^ 

0.9934 

n( 83.60) 

— 

0.13 

21.80 

16.80 

12.30 

0 

IN 

— 

0.13 

P-1.93 

16.93 

51.03 

51.03 

N°-(l-r^)K 

— 

0.13 

21.80 

16.80 

12.20 

0 

f 
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m° 

— 

0.13 

21.93 

16.93 

50.93 

50.93 

"j 

(2)8216  (9)3578 

34.53 

34.53 

0.2086 

V 

0.23BWRT  ^ 

10^'  7A1 

7.41 

7.4i 

7.4i 

7.4i 

7.41 

SN" 

IIT 

— 

1 

1 

1 

0.9980 

0.9980 

r'(83) 

0 

— 

0 

0 

0 

o^<- 

0 

*  Note :  r '  i 
0 

(a,  t)  can  not  be  < 

a  negative 

value , 

hence  the 

smallest  is 

zero. 

(1d)  10  MT  grotand  surface  burst,  1A7 

min  after  burst, 

Sr(Y)(90). 

Decay  chain: 

TJ^  _ 

T/-^  -  TVU 

ijP  ^ 

■  ivr  ■'  x\D 

»j>r  1 

Br 

Kr 

Rb 

Sr 

V 

r^(90) 

— 

0 

0.036 

0.1378 

0.1378 

1-r 

0 

— 

1 

0.964 

0.8622 

0.8622 

N(90,88) 

— 

67 

338 

106 

0 

ZN 

— 

67 

405 

511 

511 

N"=(l-r  )N 

— 

67 

325.8 

91.4 

0 

ZN" 

— 

67 

392.8 

392.8 

392.8 

(^)3050 

00 

( 

6)4338 

0 

0 

V 

0.23BWRr  ^ 

10®  1.92 

1.92 

1.92 

1.92 

1.92 

ZN” 

ZN 

— 

1 

0 . 9698 

0.7686 

0.7686 

r^(90) 

— 

0 

c 

0.0936 

0.0936 
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(c.)  10  MI  ground  surface  lurstj  1.47  min  after  burst, Ge ( 75 ) . 

Decay  chain; 

Cu - fc-  Zn - -  Ga - ^  Ge 


Cu 

Zn 

Ga 

Ge 

r„(75) 

— 

1 

1 

0.6716 

l-r 

0 

— 

0 

0 

0.3284 

n(75,88) 

— 

(3)12 

0.0459 

0.0536 

— 

(3)12 

0.0460 

0.0996 

N°=(l-r  )N 

0 

— 

0 

0 

0.0176 

— 

0 

0 

(11)2288 

V  ,^8 

0.23BWTR  ^ 

1.92 

1.92 

1.92 

1.92 

,,  ^ 

-3 

— 

(3)42a 

0.23BWTT. 

J 

— 

0 

0 

0.0018 

r'(T5) 

0 

— 

0 

0 

0.3266 

*  Note:  In  computing  Y.,  the  individual  mass  number  dependent  conversion 

factor  o:(A)  for  to  should  be  multiplied  to  each  radionuclide. 

1.  =  S.a:(A)N(A,  t).  o;(a)  can  be  calculated  from  Table  I 

.1  " 


(d. )  10  MT 

ground  surface 

burst ,  6.77  niin 

after  burst, 

Sr(89).  Decay 

chain 

was  shown  in 

Section  (a.) 

Se 

Br 

Kr 

Rb 

Sr 

r  '89) 

0  • 

— 

0.0303 

0 

0.0119 

0.0192 

l-r 

n 

— 

0.9697 

1 

0 , 9881 

0.9808 
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N(8y,4o6) 

__ 

— 

84 

267 

56 

— 

— 

84 

351 

407 

W"*(l-r  )N 

0 

— 

— 

84 

ZN" 

84 

84 

84 

CD 

0 

0 

ZN"  ZN 

1 

0.2393 

0,2064 

r  (89) 

— 

— 

0 

0 . 7^86 

0.7804 

3.  Gross  Fractionation  Nianber 

The  gross  fractionation  niaober,  r^,  or  r  .  is  defined  ty  'E-^uation  !47) 
The  two  ionization  rates,  i^p('t)  and  i^^  (t^  in  this  form'^a  are  defined  hy 
Equations  (4-5)  and  (46)  respective.'Ly.  The  computations  are  stiaightfoiward, 
though  tedious,  once  the  values  of  the  fractionation  numbers,  disintegration 
multipliers  and  activities  are  known.  Ihe  variations  in  the  ionization  rates 
of  normal  fissxon  products  with  time,  or  simply  the  decay  of  normal  fission 
products,  for  U-235j  0-238  and  Pu-239  have  been  calculated  and  are  tabulated 
in  Table  V.  From  this  table,  the  normal  ionization  rate  of  ..r-238  at  H  +  1 
hour  can  be  interpolated  and  obtained  as 

i^p(l)  =  6.973  X  10  'r/hr  at  1  hr)  /  i^iission/sq.  ft-.  '  (65 

For  a  10  MT  surf'ace  burst,  the  condensed  ionization  rater,  .for  four  di.f'terent 
firebal],  rising  times  have  been  computed  and.  are  listed  with  the  \Q.lues  of 
gross  fractionation  nmbers  computed  from  them: 


t 

53  sec 


2.408  X  19' 


-  ft  ■ 
"X  ' 

0 . 368 
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6o  sec 

3.255  X  10"^^ 

0.467 

1.4-7  rain 

3.873  X  10'^^ 

0.554 

2,15  min 

5.029  X  10“^2 

0.721 

Since  the  fireball  rising  time,  t,  is  related  to  the  particle  size  para¬ 
meter,  a,  by  Equation  (12),  is  also  a  function  of  CC,  sometimes 

denoted  as 

C .  Activity-Concentration  Intensity  Ratio 

With  Equations  (^T)  "the  activity  concentration  intensity 

ratio  for  any  isotope  at  any  downwind  location  can  be  readily  computed.  The 
values  of  the  ratio  for  six  biologically  important  isotopes  in  a  10  MT 
surface  burst  are  calcxolated  and  shown  in  Table  VI. 


TABLE  1 


Decay  of  ffomal  Fission  Products  from  U-235,  tJ-238  and  Pu-239 
(r/hr  at'  3  ft  a!bove'  an  Infinite  plane  for  10^  fissions  per  sq.  ft.) 


Age 

U-235 

U-238 

Pu-239 

Years 

Days 

Hours 

Fission 

(8  Mev) 

Thermal 

Fission 

0.763 

/ 

(8)9977^ 

(8)9960 

(8)9329 

(8)8907 

(8)8750 

1.12 

(8)6648 

(8)6632 

(8)6172 

(8)5866 

(8 

‘5761 

1.64 

(8)4149 

(8)4153 

(8)3827 

(8)3592 

(s! 

>3537 

2.40 

(8)2453 

(8)2484 

(8)2256 

(8)2071 

(8)2065 

3.52 

(8)l4lO 

(8)1450 

(8)1303 

(8)1166 

(8 

11196 

5.16 

(9)8079 

(9)84i8 

(9)7582 

(9)6642 

(9i 

17098 

7.56 

(9)4786 

(9)5014 

(9)4587 

(9)3986 

(9)4398 

11.1 

(9)2964 

(9)3094 

(9)2897 

(9)2555 

(9)2821 

16.2 

(9)1804 

(9; 

11869 

(9)1792 

(9)1626 

(9)1761 

23.8 

(10)9716 

(9 

11094 

(9)1073 

(9)1010 

(9)1063 

34.8 

(10)6305 

(10; 

16428 

(10)6393 

(10)6235 

(10)6360 

2,13 

51.1 

(10)3730 

(10)3786 

(10)3817 

(10)3869 

(10)3811 

3.12 

74.9 

(10)2276 

(10)2319 

(10)2365 

(10)2470 

(10)2362 

4.57 

(10)1483 

(10)1524 

(10)1556 

(10)1645 

(10)1546 

6.70 

(11)9986 

(10)1031 

(10)1039 

(10)1099 

(10)1021 

9.82 

(11)6774 

(11)6972 

(11)6899 

(11)7244 

(11)6655 

l4.i<- 

(11)4490 

(11)4587 

(11)4462 

(11)4650 

(11)4226 

21.1 

(11)2910 

(11)2940 

(11)2845 

(11)2953 

(11)2660 

30.9 

(11)1813 

(11)1807 

(11)1762 

(11)1337 

(11)1645 

45.3 

(11)1061 

(11)1039 

(11)1034 

(11)1092 

(12)9777 

66.h 

(12)6055 

(12)5807 

(12)5910 

(12)6360 

(12)5728 

97.3 

(12)3676 

(12)3497 

(12)3559 

(12)3896 

(12)3543 

143 

(12)2170 

(12)2090 

(12)2079 

(12)2320 

(12)2200 

208 

(12)1194 

(12)1164 

(12)1133 

(12)1287 

(12)1180 

301 

(13)4874 

(13)4790 

(13)4733 

(13)5707 

(13)5170 

1.2 

438 

(13)1399 

(13)1373 

(13)1525 

(13)2135 

(13)1864 

1.78 

650 

(14)3884 

(14)3758 

(14)5517 

(14)9083 

(14)7690 

2. 60 

(14)2031 

(14)1975 

(i4)316o 

(14)4964 

(14)4352 

3.80 

(l4)l444 

(14)1432 

( 14)2213 

(14)2692 

(14)2594 

5-58 

(14)1154 

(14)1158 

(14)1603 

(14) 1442 

(14)1611 

8.18 

(14)1026 

( 14)1021 

(14)1291 

(15)9971 

(14)1225 

12.0 

(15)9432 

(15)9293 

(14)1094 

(15)8452 

(14)1057 

17.6 

(15)8310 

(15)8211 

(15)9164 

(15)7377 

(15)91&3 

25.7 

(15)7183 

(15)6987 

(15)7431 

(15)6219 

(15)7668 

Kumters  In  parentheses  indicate  the  numher  of  zeros  hetween  tlie  decimal  point 
and  the  first  significant  figure. 


Concentration  Intensity  Ratio  [N*(A)/i(1)] 
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D.  Particle  Size  Distribution 


There  is  almost  no  data  available  on  particle  size  distribution, 
nor  has  any  rigorous  theoretical  approach  been  made.  However,  using  the  scaling 
functions  of  the  previously  discussed  simplified  fallout  model,  estimates 
can  be  made  of  the  particle  size  groups  that  fall  at  any  downwind  location. 

Given  the  fall  velocity  and  the  initial  altitude,  the  size  of 
spherical  particles  falling  through  a  standard  atmosphere  can  be  found 
in  curves  prepared  by  D.  E.  Clark. At  any  downwind  distance  X  the 
first  information  needed  is  the  maximum  and  mlraimum  cc  size  for  the  particle 
groups  arriving  at  that  location.  Because  the  formulae  for  o:  and  ct  . 

depend  on  whether  the  a  group  descends  from  the  stem  or  the  cloud,  it  was 
decided  that  the  X-coordinate  intersection  of  the  logarithmic  line  connecting 
(X^,  I^)  "to  (X^^,  I^)  (point  3  to  point  ^4-  in  Figure  4)  with  the  logarithmic 
line  from  (X^,  I^)  to  (Xg,  Ig)  -  extended  if  necessary  -  would  be  the 
division  mark.  For  any  downwind  distance  less  than  this  intersection  X 
coordinate,  the  source  of  the  particle  group  is  assimied  to  be  the  stem. 

Then  o:  and  a  .  for  the  stem  may  be  solved  by  graphical  methods  or 

by  the  following  computer  approach.  Essentially  a  Newton-Raphson  iterative 
procedure  is  performed  on  the  following  function  to  find  separately 


log  [X  - 


V  P 

a(Z  a  - 


V 


k  is  the  exponential  constant  for  the  fireball  major  axis  expansion. 
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is  an  empirical  inverse  time  constant,  0.011  sec"^, 

V  is  the  wind  velocity,  usually  assigned  to  he  10  mph, 
w 

p  is  the  constant  in  the  empirical  formula  of  V  /V„,  0.95  when 

Z  X 

Z  =  5,000  to  50,000  ft.,  d  «  200  to  1,200  microns, 

q  is  the  constant  in  the  empirical  formula  of  V^/V^,  1.02  x  lO"^ 
when  Z  =  5^000  to  50,000  ft.,  d  -  200  to  1,200  microns,  and 
is  the  major  semi-axis  of  the  fireball  at  ground  zero. 

The  iterative  procedure  is 


a  sa  initial  estimate 

0 

(67) 

> 

I 

II 

.H 

+ 

(68) 

If  the  sequence  of  a's  approach  a  limit,  then  this  a  is  such  that  f(a)  *  0; 
or  a  is  approximately  equal  to  or  a^_.^  depending  on  the  sign  of  the 

number  enclosed  in  the  absolute  value  symbol  of  Equation  (56). 

When  the  downwind  distance  is  greater  than  or  equal  to  the  X 
Intersection  coordinate,  it  is  assumed  that  the  fallout  originates  from 
the  cloud.  For  particle  groups  falling  within  the  cloud  area  at  any  (X,  Y) 
location,  the  maximum  and  minimum  a  values  can  be  calculated  by: 

a  =  hX  t*/  X^b^d  -  yVa^)  +  (a^  -  y^)  h^  -  b^(l  -  y^/a^) 

“  h^  -  b^  (1  -  y2/a^) 


where  Y  is  translated  to  y  by  Equation  (57). 

For  any  a  value,  the  fall-rate  is  given  ^7  ~  • 

for  the  stem  fallout  is 


z 


(z  a  k  -  V  p) 
o  z  w 

a  k  +  V  q 
2 


The  height  of  fall 


(70) 
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For  the  particle  groups  that  originate  from  the  cloud  the  maximun  and 
minimum  heigjat  of  fall  is  given  by 

a(X  -  ah)b^  ±  ah  '/  (a^  +  Qph^)(l  y^/a^)  -~(X  - 

^  “  2  2^2 

a  +a  b 

Using  this  Information  calculated  in  the  particle  size  parameter  program, 
Figure  T ,  and  a  set  of  curves  prepared  by  D.  E.  Clark,  see  Figure  8,  the 
approximate  particle  size  distribution  at  (X,  Y)  can  be  estimated. 

For  a  10  MT  detonation,  and  for  various  dovmwind  distances  the 
maximum  and  minimum  d;  the  falling  velocity,  V^;  the  initial  altitude  Z 
and  particle  diameter  p.  are  given  in  Table  VII. 

( IP ) 

Recent  work  by  Clark  and  Cobbin-  is  closely  related  to  the 
subject  of  this  report.  These  authors  used  the  same  fallout  model  as 
described  herein  for  the  calcination  of  particle  size  and  radiation 
intensity  at  varying  distances  from  ground  zero.  The  results  obtained 
are  in  good  agreement  with  those  herewith  reported  for  weapons  of  similar 


yield. 
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Figure  7-  Computer  Program  for  Pn-H-Jr-iF.  FSizR  Parameters 
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IV.  INDUCED  RADIOACnVITI  IN  SOII^ 


Elements  can  te  inade  radioactive  ty  the  action  of  thermal  neutrons 
which  e;:cite  the  stable  atom  making  it  unstable  and  radioactive .  As  thermal 
neutrons  are  liberated  In  a  nuclear  explosion,  there  will  be  some  Induced 
radioactivity  in  the  soil  around  the  detonation  site .  A  small  portion  of 
this  soil  -will  be  taken  up  into  the  fireball  of  the  bomb  (for  a  ground  burst 
or  near  ground  burst)  and  become  part  of  the  total  fallout. 

(13) 

A..W.  Klement'''^  '  in  his  report  concerning  the  potential  radionuclides 

produced  in  weapons  detonations,  lists  the  activities  of  the  various  nuclides 

2  6 

that  are  induced  in  the  soil,  A  total  of  10  neutrons  per  KT  was  assumed  to 

27 

be  liberated,  which  agrees  within  a  factor  of  15  with  that  of  1,5  x  10 

(14) 

neutrons  per  KT  calculated  from  Glasstone  .  The  composition  of  the  soil 
that  was  used  for  these  calculations  is  given  in  Table  VIII. 

Figure  9  gives  the  decay  of  the  induced  radioactivity  in  soil.  The 
data  was  obtained  by  taking  the  contributions  of  the  various  nuclides  given 
in  Table  IX  and  summing  them  for  the  various  times  after  detonation  up  to  one 


year. 

Figure  10  ,  "integrated  Neutron  Flux  as  a  Function  of  the  Slant 

Range  in  Air  of  O.9  Sea-Level  Density  for  a  l-Kiloton  Explosion"  shows  the 

neutron  flux  versus  distance  for  a  1  KT  bomb.  These  curv'es,  as  given  by 

(Ik) 

Glasstone'  ,  were  integrated  over  an  area  with  a  2,500  yard  radius  and 
were  used  to  correlate  Klement's  va-kues . 

Senftle  and  Champion^ give  a  detailed  discussion  concerning 
induced  radioactivity.  A  method  for  calculating  the  induced  radioactivity 
involving  tables  and  formulas  is  also  g'ven.  This  method  was  used  to  prepare 

'’13) 

the  data  in  Table  X  in  order  to  correlate  the  data  given  by  Klement ' 


and  shown  in  Table  IX. 


Figure  9-  Total  Activity  of  Soil 


-78- 


TABLE  IX 

Radioactivity  in  Soils  (Mc/MT) 


Nuclide 


0 

1  hr 

1  day 

1  month 

1  year 

Na^^ 

1620 

1430 

542 

a 

a 

584 

7.05 

a 

A1^8 

83^870 

15.9 

a 

si3^ 

991 

773 

1.78 

a 

32 

1.21 

1.21 

1.15 

0.276 

a 

K 

a 

204.5 

194 

50.0 

a 

Ca^5 

0.300 

0.300 

0.300 

0.263 

0.057 

1100 

135 

a 

Tl^^ 

156 

0.0015 

a 

Mn^^ 

2122 

1655 

3-82 

a 

Fe55 

0.110 

0.110 

0.110 

0.104 

0.087 

Fe^^ 

0.014 

0.014 

0.014 

0.0095 

0.00006 

a  =  less  than  one  curie  per  megaton 


The  form’ula  for  the  radioactivity  Induced  in  any  element  as  given 
( 15 ) 

by  Sent  Me  and  Champion^  '  is; 


.  (c  rax)(l  -  e/®) 


(72) 


where 


activity,  in  disintegrations  per  second,  after  the  nuclide 
has  "been  removed  from  the  flux  for  a  period  6 
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0  SOO  1,000  1,500  2,000  2,500 


SLANT  RANGE  FROM  EXPLOSION  (YARDS) 


Figure  10.  Integrated  Neutron  Flux  as  a  Function  of  the  Slant  Range  in  Air 
of  0.9  Sea-Iievel  Density  for  a  1-Klloton  Explosion 
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Talile  X 

Activities  of  Various  Nuclides  la  Soy. 


Nuclide 

Induced  Activityj 

Na 

2li-,800 

27 

Mg 

11,100 

ro 

CD 

1,250,000 

12,300 

18 

CM 

3,350 

p  ^5 

Ca 

4.67 

l6,400 

4,730 

MD 

ITS 

i 

36,200 

Fe55 

0.580 

Fe59 

0.665 

Assumes  HTS  soil 

At  one  minute  after  detonation;  decay  time  equals  zero 
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cr 

f 

N 

k 

t 

Q 

T 


activation  cross  section  in  square  centimeters  for 
2200  m/sec.  neutrons j 

thermal  flux  in  neutrons  per  square  centimeter  per 
second; 

the  total  number  of  atoms  of  the  element  in  the  target; 
relative  abundance  of  the  isotope  from  which  the 
radionuclide  is  fonned; 
time  of  irradiation; 
time  of  decay; 

half-life  of  the  radionuclide  formed. 


The  product  (o  fNk)  in  Equation  (72)  has  been  calculated  for  a  flux 

of  10"^  neutrons/sq.cm- sec  and  for  the  number  of  atoms  in  1  g  of  the  target 

element.  The  values,  designated  A  ,  are  listed  for  each  element  by  Sesnftle 

s 

( 1^) 

and  Champion'  .  The  values  for  can  be  corrected  to  any  neutron  flux 
and  weight  of  target  element  by  multiplying  by  a  suitable  factor  of  propor¬ 
tionality.  The  corrected  value  of  A  is  given  by: 

s 


A’- 

s 


10 


A 

12^  s 


(73) 


where  W  is  the  weight  of  the  element  in  the  sample,  and 
f  is  the  neutron  flux. 

The  formula  used  in  the  calculations  was: 


A 


t 


A'  (1  -  e'°-^93t/T  ^  (g-0.693  e/T^ 


taking  t  =  1  min  and  9  =  0. 
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As  a  special  case,  where  t  is  known  to  be  less  than  15^  of  T, 
a  further  simplified  foimiula  was  used. 


0.693)  (e-^'693  s/T) 


Values  used  for  A^  obtained  by  Senftle  and  Champion^  are  listed 

below : 


TVT  2h 

Na 

TO 

1.39  X  10 

Ca**** 

1.78  X  10® 

1.3  X  10® 

2.5!;  X  10'^ 

Al^S  . 

h.79  X  10^ 

Ti^^ 

9.01  X  lo'^ 

Si3l  _ 

e.-jh  X 

Mn5® 

-  1.39  X  IC^^ 

- 

3.69  X  10^ 

Fe55 

-  h .  56  X  10® 

1.02  X  10^ 

Pe59 

-  2.26  X  lo”^ 

In  calculating  the  data  of  Table  X  ,  a  neutron  flux  of  O.9I  x  10^° 
neutrons/sq  cm  per  megaton  was  assumed.  The  total  weight  of  the  soil  involved 
was  taken  as  5-92  x  10^^  grams.  This  value  was  calculated  using  a  £500  yd 
radius,  a  1-ft  depth,  an  average  density  of  I.18  gm/cu  cm,  and  with  a  sci]. 
composition  as  that  given  in  Table  7131.  The  irradiation  time  (t)  was  taken 
os  one  minute . 

The  data  presented  in  Table  X  agree  reasonablv  well  with  the 
data  presented  by  Kleinenh  in  Table  IX.  Tnis  agreement  is  within  an  expected 
factor  o:  15  s'?  the  neutron  flux  used  was  a  factor  of  15  higher  bhan  that 
assumed  by  Klement, 

Manieville  '  states  that  thr  radioactivity  of  the  fission  products 
from  a  nuclear  explosion  overshadows  the  induced  radioactivity  by  a  factor  of 
approxiiiiately  10^ .  E'/en  though  this  tends  to  make  the  contribution  of  the 
induced  radioactivity  almost  negligible  there  are  conditions  under  which  the 
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induced  radioactivity  might  play  a  detectalale  part.  Tiie  induced  nuclides 
that  form  soluble  salts  will  be  readily  dissolved  in  both  surface  and 
ground  water  and  therefore  will  add  to  the  contamination  of  the  water 
supply . 

It  has  been  generally  recognized  that  the  activating  neutrons 
will  penetrate  the  soil  before  the  blast  wave  arrives.  Some  of  the 
activated  soil  will  then  be  talcen  up  into  the  cloud  due  to  the  cratering 
effect  and  thereby  becomes  part  of  the  total  fallout.  According  to 
Dr.  C.  F.  ^'•iller  ^  the  induced  activity  In  the  soil  will  be  about 
0.019  (l"9  pel'  cent)  of  the  ior,‘;.a'Mon  rate  due  to  normal  fission  [l^.,  (t)]. 

ip 

This  fraction  will  be  a  small  percentage  of  the  total  ionization  rate . 

This  esrimate  of  the  induced  activity  is  applicable  for  a  100^  fission 
yield,  assuming  0.8  neutron  capture  per  fission  from  U-2^8,  regardless 
of  weapon  size . 

It  has  also  become  established  that  the  induced  activity  will 
piedominate  for  a  t^tal  fusion  bomb.  For  a  50^  fusion  and  50^  fission 
bomb  ,  the  induced  activity  will  be  re.lativelv  jmal.i ,  except  for  a  short 
peri'-.i  Immediately  after  detonation  However,  activs.ted  sodium  and 

“ly  contribute  materially  at  a  later  tirae^^'''. 

It  is  therefore  concluded  that  the  contribution  of  induced 
radioac..ivlty  to  the  contamination  of  water  supplies  will  not  be  signi¬ 
ficant  . 
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V.  TRANSPORT  OP  FALLOUT  PARTICLES  BY  SURFACE  WATER 

Transport  of  the  insolutle  portion  of  fallout  may  be  divided 
into  four  general  phases-  transport  in  air,  transport  in  overland  flow, 
transport  in  stream  flow  and  transport  in  reser  vr.-xrs .  Each  of  the  four 
phases  may  be  treated  as  a  separate  tcpic.  Of  the  latter  three,  only 
transport  by  stream  flow  has  been  considered  in  any  detail;  overland 
transport  and  reservoir  transport  are  mentioned  only  briefly.  Transport 
in  air  is  discussed  in  Section  II  of  this  report. 

There  are  three  major  limitations  on  any  stream-flow  transport 
analysis.  These  are:  (l)  the  complex  nature  of  sediment  behavior  in 
streams;  (2)  the  paucity  of  data  concerning  fallout,  especially  data 
concerning  particle  size  distribution  or  the  relation  of  radioactivity 
to  particle  size;  and  (3)  the  obvious  impossibility  of  obtaining  stream- 
flow  data  for  every  reach  of  sti-eam.  These  limitations  precluded  an 
analj'-sis  of  the  quantity  of  sediment  resulting  from  the  fallout  of  a 
given  burst,  although  work  on  the  quantity  of  fallout  is  now  being  done. 

It  could  not  be  determined  that  a  certain  weight  of  sediment  would  fall 

on  a  water  surface  or  what  the  particle  size  distribution  would  be. 

(l4) 

Computations  are  available,  however,  from  Glasstone  and 
the  OCD  fallout  model  which  gives  certain  relationships  between  particle 
sizes,  downwind  distance  carried,  and  total  activity. 

Glasstone  shows  that  99  per  cent  of  the  total  activity  of 

fallout  is  carried  by  particles  of  kOO  microns  or  less  in  diameter. 

Tnis  relation  of  particle  size  to  total  activity  was  i.terlved  from 

( l4) 

Figure  9.18?  of  Glasstone 's  book'  '  and  is  shown  as  Figure  11  of  this 
report.  The  method  contained  in  this  section  is  based  on  this  relation¬ 
ship. 

fRECSQUD  PAOt  BUMUNOT  nUOl 
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E’articles  of  kCO  micron  diameter  are  on  the  order  of  magnitude 
of  normal  stream  sediment.  This  means  that  the  particles  carrying 
essentially  all  the  activity  will  remain  in  suspension,  or  at  least 
will  not  settle  out  immediately.  It  follow.s  that  most  of  the  fallout 
particles  which  land  on  the  water  surface  or  find  their  way  into  the 
stream  from  overland  transport  will  he  carried  hy  the  stream  and  will 
present  a  hazard  at  some  distance  downstream. 

A  relationship  of  particle  size  versus  downwind  distance 
carried  was  computed  in  Section  II  of  this  report.  (Table  VIl).  Data 
from  that  table  has  been  plotted  in  this  section  as  Figure  12. 

ihe  following  general  assumptions  were  made: 

(1)  Upt8.ke  hy  biological  organisms  is  negligible. 

(2)  The  stream  chemistry  is  not  such  as  to  dissolve  ordinarily 
insoluble  particles. 

(3)  Effect  of  flocculation  on  settling  rates  is  negligible. 

(4-)  Weapon  size  wind  speed  (l5mph)j  and  fallout  pattern 

are  u.sed  to  correspond  with  other  sections  of  this  report. 

(5)  Radioactivity  of  particles  over  500  microns  may  be  ignored. 

Further  assumptions  are  made  and  discussed  In  the  text.  Some 
of  there  assumptions  depend  on  the  natiur-e  of  a  specific  stream,  some 
require  further  investigation,  and  some  which  are  considered  negligible 
in  an  approximate  analysis  may  assume  a  proportionally  greater  importance 
in  cj,  more  detailed  computation. 

An  estimate  must  be  made  of  stream  width,  vol’jme,  and  surface 
area .  This  tiisy  be  done  by  equating  the  cross-sectional  area  of  the 
rtreara  1  at  a  gaging  rta+lon)  to  a  rect^angle  of  equivalent  area  whose  length 


CUMl/LATIVE  PERCENTAGE  OF  TOTAL  ACTIVITY 
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is  eqml  to  the  stream  width.  This  is  shown  in  Figure  13.  The  length 
of  the  stream  between  any  two  sections  is  considered  to  be  a  straight 
line  equal  to  the  distance  along  the  stream  channel  between  the  two 
sections ,  Each  section  of  the  stream  may  then  be  considered  to  be  a 
tnmcated  pyramid  lying  on  its  side  with  its  bases  the  equivalent  rec¬ 
tangles.  (Or  for  the  section  containing  the  source,  a  complete  pyramid, 
since  the  source  is  considered  to  be  a  point.)  Concentrations  may  be 
found  using  the  volume  of  the  pyramids. 

Intensity  (roentgen/hour)  contours  for  a  given  burst  are  super¬ 
imposed  on  the  watershed.  The  intensities  may  then  be  converted  to  activities 
(atoms/sq.  ft.)  for  insoluble  (and  soluble)  portions  of  the  various  isotopes 
under  consideration  by  use  of  Table  VI.  Details  of  this  method  are 
discussed  on  page  95  through  page  98  of  this  report. 

If  fallout  landing  on  the  water  surface  is  assumed  to  mix  uniformly 
within  the  reach  of  channel  enclosed  by  a  given  contour,  an  activity  con¬ 
centration  may  be  computed  using  the  idealized  channel  shape.  This  concen¬ 
tration  will  move  downstream  as  a  unit,  and  the  time  of  arrival  and  time  of 
lapse  of  this  concentration  may  be  found  at  a  point  downstream  by  considering 
the  particles  to  move  at  mean  stream  velocity. 

By  plotting  the  various  concentrations  from  units  under  given 
intensity  contours  as  ordinate  against  time  as  abcissa,  a  series  of  horizontal 
bars  is  obtained.  These  may  be  averaged  graphically  to  obtain  a  time- 
concentration  curve  for  the  stream.  Note  that  concentration  is  expressed 
as  activity  rather  than  particle  concentration. 

By  use  of  Figure  12  a  maximum  particle  size  at  a  given  point 
downwind  from  the  blast  may  be  found.  The  maximum  size  particle  landing 


a  --  channel  ien9th  berween  sections 


W  *  stream  width  «  length  of  eqwivolcnt  rectongle 
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Vs  “  woter  surface  between  sections 


Figure  13-  Idealized  Stream  Section 


on  the  stream  reach  under  consideration  is  compared  to  the  maximum  particle 
size  which  the  stream  will  carxy  in  suspension.  If  the  maximum  size  jx'.rticle 
is  carried  then  it  is  to  be  assumed  that  smaller  particles  will  also  he 
carried,  and  the  method  remains  valid. 

If  the  stream  will  not  carry  the  maximum  size  particle,  then 

some  percentage  of  the  activity  must  he  regarded  as  going  into  hed  load, 

which  travels  at  a  much  slower  rate.  Ilhe  only  way  now  available  to  find 

the  proportion  of  activity  which  would  go  into  bed  load  is  hy  assuming 

that  if  particles  of  a  given  size  are  found  to  settle,  that  all  particles 

of  this  size  or  larger  will  settle  out,  and  all  smaller  particles  remain 

in  suspension.  The  relationship  of  particle  size  to  activity  given  hy 
fll^) 

Glasstone^  ^  may  then  be  used  to  find  what  percentage  of  the  total  activity 
would  settle  out  with  the  larger  particles. 

The  ability  of  a  stream  to  carry  particles  of  a  civen  size  or 
smaller  is  best  defined  '  y  what  may  be  termed  the  effective  carrying 
velocity  of  the  stream,  i.e.  that  velocity  below  which  settling  of  particles 
of  a  given  size  is  most  likely  to  occur.  The  effective  carrying  velocity 
is  obviously  not  related  to  the  Idealized  stream  channel  used  to  compute 
coni-entrations,  hut  is  rather  a  function  of  both  channel  characteristics 
(actual  channel  width  and  depth)  and  flow,  as  well  as  particle  size  and 
sediment  distribution.  Stream  flow  data,  no  matter  how  accurate  or  com¬ 
prehensive,  cannot  cover  every  mile  of  channel.  Conditions  which  will 
produce  the  greatest  possibility  of  settling  must  therefore  be  generaliza¬ 
tions  of  what  is  known  of  the  nature  of  an  indivld'unl  stream. 

Particle  size  distribution  in  air  is,  unfortunately,  not  known, 
except  for  maximum  and  minimum  distance  which  a  particle  of  given  size  will 
travel.  There  is  no  way  at  present  to  compute  the  number  of  particles 


per  unit  area  landing  on  the  watershed.  Sediment  concentrations  cannot, 
tlierefore,  be  computed.  This  is  the  reason  for  computing  concentiatton 
in  terms  of  activity.  Adiitioral  fallo;,t  dat?,,  parti c^alarly  on  particle 
di.- ti'ibution,  and  additional  study  on  sediment  concentrations  are  needed 
to  support  the  assumptions  that  the  fallout  particles  mix  uniformly, 
that  particles  will  go  out  of  suspension  according  to  size,  and  that 
particles  move  downstream  at  mean  stream  velocity. 

In  order  to  approximate  actual  stream  conditions,  close  scrutiny 
of  maps  and  existing  stream-flow  records  will  be  required.  Large-scale 
topographical  maps  should  yield  much  irf orma.tion  on  channel  conditi.ons 
and  ctreara-bed  slopes.  The  L.  S.  Geological  Suiweys  and  Corps  of  Engineers 
publish  or  have  on  file  not  only  gaging-station  records  but  much  special 
measurement  data. 

:'\o  analysis  has  been  made  for  transport  res-olting  from  overland 
flow.  Fallout  landing  on  the  land  portion  of  the  watershed  is  considered 
less  of  an  immediate  problem  than  that  landing  on  the  water  surface, 

and  will  probably  present  no  immediate  problem  unless  It  is  raining 

at  time  of  falioiit  arrival,  (in  which  case,  incidentally,  the  fallout 
may  be  intensified,;  Subsequent  rain  may  produce  a  second  activity 
peak  st  a  water  intake  by  washing  particles  into  the  stream.  It  is  of 
interest  to  note  that  at  Coweeta  Hydrologic  Laboratory  it  was  found  that 
much  of  the  r-.noi'f  which  was  i.’orraerly  considered  to  flow’  overland  actually 
gcf.f  ur.der  the  surface,  or  altametes  between  sub-surface  and  above-surface 
flow.  I.a  woi'k  with  ividloactlve  tracers  at  Coweeta,  great,  difficulty  was 
eeperienced  in  t lading  a  tracer  which  would  follow  the  water  at  the  same 
speed  as  the  w.ater.  Mos*;  of  the  isotopes  used  at  Coweeta  tended  to  stay 

nf.ur  the  poinl  wheie  '.hey  were  applied.  These  studies  were  carried  out 
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on  steep  slopes  on  which  there  was  a  higher  percentage  of  surface  runoff 

than  on  the  same  soil  at  a  lesser  slope;  however,  the  isotopes  used  were 

( 19^ 

not  ah sorted  within  particles.^ 

Particles  which  go  into  hed  load  may  present  a  long  term  hazard 
"because  of  the  slower  and  more  irregu-lar  movement  of  bed  load.  A  bed 
load  analysis  will  be  required  if  it  is  found  that  large  particles,  in 
spite  of  their  possession  of  only  a  small  percentage  of  the  total  activity, 
are  found  to  actuially  carry  a  dangerous  amount  of  activity.  The  smaller 
the  particle  size  that  is  found  to  settle,  the  more  significant  becomes 
the  bed  load  analysis.  Note  that  such  a  bed-load  analysis  complements 
assumption  5  on  page  86, 

Transport  in  reservoirs  is  more  likely  to  produce  settling 
because  reservoirs  are  more  likely  to  approach  quiescent  conditions. 
Calculations  will  probably  be  simplified  because  the  conditions  of  reser¬ 
voir  transport  more  closely  approaches  the  condition  of  settling  in  a 
quiescent  basin  than  do  the  other  phases  of  the  transport  problem. 

This  method  was  presented  not  only  because  it  is  the  one  most 
likely  at  present  to  produce  usable  results,  but  because  it  provides  an 

outline  of  the  problem  of  sti'eam  transport.  Assumptions,  questions  as 
to  their  validity,  and  areas  requiring  further  study  have  been  pointed 
out.  Computations  have  not  been  included  because  it  is  wished  to  ascer¬ 
tain  the  validity  of  some  of  the  assumptions,  and  because  sufficient 
study  has  not  been  made  on  specific  streams.  This  treatment  therefore 
constitutes  only  a  good  start;  however,  further  study  in  the  same  manner 
should  produce  results  of  known  validity. 
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VI .  EVALUATION  OF  MUNICIPAL  WATER  SUPPLY  COITOAMINATIOW 

A.  Introduction 

An  approach  to  the  problem  of  water  contamination  by  radioactive 

fallout  was  taken  which  would  give  a  more  reliable  set  of  concentration 

values  than  those  calculated  previously.  The  values  presented  in  Quarterly 

* 

Technical  Report  No.  3  were  maximized j  in  that  they  represented  a  case  in 
which  the  highest  concentration  that  could  reasonably  be  expected  at  H  +  1 
hour  from  any  kind  of  attack  or  environment,  except  for  possibly  a  very  high 
megatonnage  attack.  Following  is  presented  a  more  detailed  and  realistic 
approach  to  the  problem.  The  values  herewith  reported  represent  specific 
cases  in  which  environmental  and  other  factors  have  been  considered,  and 
it  is  believed  they  are  at  least  of  the  same  order  of  magnitude  as  would 
be  expected  following  a  nuclear  attack. 

Only  the  cities  of  Houston,  Texas  and  New  York  were  chosen  for 
a  full  scale  evaluation.  Since  the  values  calculated  for  these  two  water¬ 
sheds  were  in  good  correlation,  it  was  felt  that  the  results  could  probably 
be  applied  without  serious  error  to  the  other  cities  previously  considered. 
Excellent  watershed  data  was  available  on  both  cities,  thus  facilitating 
the  calculations  considerably. 

A  physical  integration  was  performed  over  the  considered  watersheds 
and  reservoirs  to  give  realistic  specific  isotope  concentrations  at  H  +  1 
hour  for  variously  directed  15  raph  winds.  As  before,  the  case  of  fallout 
contamination  from  that  falling  in  the  reservoir  alone  was  considered  as 
well  as  the  case  of  runoff  contEuninatlon  from  the  entire  watershed. 

B.  Calculations 

Since  large  scale  maps  of  both  the  Houston  and  New  York  watersheds 
were  available,  transparent  overlays  were  superimposed  over  the  watershed 


See  Appendix  A 
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ai-ea-E  end  appropriate  grids  scaled  in.  These  grids  consisted  of  five  'mile 
squares  for  watershed  areas,  and  two  mile  squares  for  the  reservoir  areas 
when  large  enough  to  appl/.  Figure  It  shows  the  grid  on  the  watershed 
serving  Iiouston.  The  groiuid  cero  of  a  5  MT  weapon  detonation  was  located 
so  as  to  give  a  fallou.t  pattem  which  fell  -on  the  watershed  for  the  15  mph 
wind  model.  Various  wind  directions  were  also  chosen  so  as  to  give  some 
hasis  for  comparison  and  to  account  for  possible  seasonal  variation.  It 
is  otviou.s  from  the  geography  of  li'ev  York  and  Iiouston  that  the  wind 
directions  from  the  east  and  north  wo'old  give  rise  to  a  negligible  amount 
of  fallout  over  their  respective  watersheds,  so  that  these  wind  directions 
were  omitted. 

Ihe  attacks  that  were  evaluated  were  as  follows: 

Case  la  Houston  -  Cro’jnd  zero  is  downtown  Houston  and  the  wind  direction 

is  from  the  south. 

Case  lb  Houston  -  Ground  zero  is  downtown  Hous^ton  and  the  wind  direction 

Is  from  the  west. 

Case  Ila  New  York  -  Ground  zero  is  Central  Park  and  the  wind  direction 

is  from  the  south. 

Case  lib  New  York  -  Ground  zero  la  Elngh-amton.,  Kew  York,  and  the  wind 

direction  is  from  the  west. 

Separate  eval.uations  vere  raa-de  for  each  of  these  cases  with  respect 
to  contamination  from  v.he  reservoir  and  that  from,  runoff  of  the  entire  water- 
f’.hed.  In  th'r  case  of  Houston,  the-  concentration  of  activity  was  calculated 
for  one  of  the  ‘"ee.1  strp-M,ras  from  dlre-ct  fallout  contamination  for  coTnp.j,rison 
purposes . 

For  the  r'eoervoir  suudy,  tne  urea  of  the  respective  reservoirs 
wiihir.  each  two  mile  grid  was  found  by  planimeter.  ]h  the  Houston  watershed 


there  is  only  one  reservoir  and  in  the  New  York  watershed  all  six  reservoirs 
empty  into  one  so,  that  the  problem  is  somewhat  simplified  if  comp^^te  mixing 
is  assTjimed.  The  centroid  of  each  selected  grid  area  was  estimated,  and  its 
coordinates  measured  with  respect  to  groruad  zero.  An  intensity  was  then 
assigned  to  each  set  of  coordinate%  and  it  was  assumed  to  be  constant  over 
that  specific  area.  The  intensity  values  were  chosen  accoi'ding  to  whether 
the  specified  area  was  either  covered  by  the  upwind  or  downwind  cloud. 

These  values  related  to  the  upwind  cloud  were  aridved  at  by  interpolation 
between  previously  estimated  intensity  contours,  and  those  of  the  downwind 
cloud  were  arrived  at  by  interpolation  of  computer  calculated  contours , 

To  convert  from  the  Intensity  over  each  area  to  activity  in 
atoms/sq  ft  in  each  area,  the  value  of  from  the  fallout  model, 

for  each  isotope  considered,  was  taken  as  its  value  corresponding  to  the 
downwind  distance  X.  By  multiplying  the  appropriate  N  value  by 

the  intensity  over  each  square,  the  concentration  in  atoms/sq  ft  is  found 
for  each  area.  By  then  multiplying  each  atom  concentration  by  its  cor¬ 
responding  area,  the  total  number  of  atoms  in  each  square  of  the  grid  is 
found.  The  total  number  of  atoms  of  each  isotope  at  H  +  1  hour  in  the 
reservoir  was  then  found  by  adding  up  the  contributions  from  each  square 
over  the  reservoir.  By  assuming  that  complete  mixing  has  occurred,  the 
concentration  of  activity  in  the  reservoir  in  atoms/liter  may  be  obtained 
by  dividing  by  the  total  volume  of  water  in  the  reservoir.  The  concen¬ 
tration  of  each  isotope  considered  at  H  +  1  hour,  for  the  Houston  and 
New  York  reservoirs  is  shovjn  in  Tables  ll  and  Xj  along  with  the  eoulvalent 
activity  in  nc/ml. 
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Table  XI 


Activity  Concentrations  for  Direct  Contamination 


of  the  Hoiiston  Reservoir 


South  Wind 


West  Wind 


Isotope 


Sr-89 

Sr-90 

Ru-106 

1-131 

Cs-li? 

Ba-l40 


Atoms/ liter 


4.6  X  10^ 

10.9  X  10^ 
5.1  X  10^ 
9.8  X  10^ 
5-3  X  10^ 

13. T  X  10^ 


(Ac/ml 


1.9  X  10 

2.3  X  lO’’^ 
3.0  X  10“ 

— ii. 

2.7  X  10 

1.4  X  10"'^ 


.3  X  10" 


Atoms/ liter 


4.6  X  10^ 

10.9  X  10^ 
5.1  X  10^ 

9.8  X  10^ 

5.3  X  10^ 

13.7  X  10^ 


hc/ml 


1.9  X  10"^ 

2.3  X  lO"”^ 

3.0  X  10“^ 

2.7  X  10"^ 

1.4  X  lO""^ 


2.3  X  10‘ 


Table  XXi 


Activity  Concentrations  for  Direct  Contamination 


of  the  Few  York  City  Reservoirs 


Isotope 


Sr-89 

Sr-90 

Ru-106 

1-131 

Cb-137 


South  Wind 


Atoms/liter 


2.9  X  10^° 
5.6  X  10^° 

2.5  X  10^° 

5.6  X  10^° 


3.7  X  10-“ 


1.2  X  10 


1.2  X  10 


1.5  X  10 

1.5  X  10‘ 


7.3  X  10 


West  Wind 


Atoms/liter 


3.9  X  10^° 

6.9  X  10^° 

3.0  X  10^° 

6.8  X  10^^ 


5.1  X  10 


1.6  X  10 

1.5  X  10" 

1.8  X  10" 

.1.9  X  10" 


1.0  X  10 


Ba-l40 


7.0  X  10 


1.2  X  10 


8.4  X  10 


1.4  X  10 
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The  activity  concentrations  for  the  Houston  reservoir  (Table  Xl)  are 
the  same  for  south  wind  and  west  wind  because  the  same  intensity  contour  falls 
across  the  center  of  the  reservoir  in  both  cases.  When  the  wind  is  from  the 
souths  the  ri^t  hand  portion  of  the  contour  crosses  the  reservoir.  With  wind 
from  the  west,  the  left  hand  portion  of  the  contour  is  involved. 

A  separate  evaluation  of  stream  contamination  was  made  using  a 
randomly  selected  creek  feeding  the  Houston  reservoir  for  the  purpose  of  com¬ 
parison  only.  This  selected  feed  stream  (Cypress  Creek)  ran  perpendicular  to 
the  downwind  fallout  axis  and  emptied  into  the  Houston  reservoir.  The  evalua¬ 
tion  was  carried  out  in  the  same  manner  as  that  done  on  the  reser'/oir.  Wie 
stream  was  assumed  to  approximate  an  isoceles  triangle  with  the  base  repre¬ 
senting  the  width  of  the  stream  mouth.  Two  mile  lengths  were  laid  off  and  the 
trapezoidal  areas  calculated.  Intensity  values  and  activity  conversion  values 
were  chosen  as  before  to  obtain  the  total  activity  falling  directly  on  this 
stream.  To  present  an  estimation  of  the  concentration  of  activity  in  the 
stream,  complete  mixing  was  assumed.  However,  it  is  evident  that  the  activity 
will  really  move  in  a  slug  type  of  formation,  completely  mixing,  if  at  all, 
only  vrhen  the  radioactivity  reaches  the  reservoir.  Although  some  mixing  during 
stream  flow  will  occur  due  to  turbulence  and  diffusion,  the  presented  case 
will  never  be  actually  attained,  since  countercurrent  mixing  is  unlikely. 

Prom  the  presented  values  in  Table  XIII  it  may  be  seen  that  some  build-up  of 
activity  may  occur  in  areas  immediately  surrounding  the  stream  feed-in  points . 
However,  this  build-up  will  not  represent  a  serious  problem  in  the  studied 
case  because  of  its  distance  from  the  municipal  Intake,  and  the  relative  time 
elapsed  from  the  time  of  detonation  to  the  intake  time. 

The  second  case  studied  was  that  of  the  contamination  to  be  expected 
assuming  runoff  from  the  entire  watershed.  This  study  was  made  for  both  the 
New  York  and  Houston  watersheds,  the  method  being  simiJ.ar  to  that  of  the  first 
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Table  XIII 

Activity  Concentrations  for  Cypress  Creek  (South  Wind) 

Isotope  Atoms/liter  4c/ml 

Sr-89 

Sr- 90 

Ru-106 
1-131 
Cs-137 

Ba-li)-0 

case.  Five  mile  grids  were  used,  however,  instead  of  two  mile  grids.  Wind 
conditions  were  selected  Identical  to  those  of  Case  I.  To  calculate  the  actual 
amount  of  radioactivity  (soluble  fraction  only)  that  reaches  the  reservoir, 
the  runoff  coefficients  supplied  by  the* respective  municipal  water  works  were 
assumed  valid  in  that  they  represent  a  maximum  value  to  be  expected.  The 
calculated  activities  are  therefore  considered  to  be  maximum  values.  The 
actual  radioactive  runoff  coefficient  will  actually  be  less  than  the  aqueous 
runoff  coefficient  depending  on  (i )  the  instanteous  moisture  content  of  idle 
soil,  (2)  the  duration  of  time  from  detonation  to  rainfall,  (3)  ion-exchange 
and  absorption  in  and  on  the  soil,  and  (i)-)  plant  uptake.  The  calculated 
values  of  activity  concentrations  of  the  selected  isotopes  at  H  +  1  hour, 
assuming  complete  dilution  by  the  composite  lakes  and  streams  in  the  watershed 
for  both  Houston  and  Hew  York  is  presented  in  Tables  XIV  and  XV.  The  factors 
Involved  obviously  vary  widely  from  one  environment  to  another,  so  that  general 
assumptions  made  from  any  particular  reservoir  should  be  applied  with  caution 
to  other  situations . 


2.2  X  10 


12 


5.^4-  X  10 

2.6  X  10 

4.4  X  10 

2.5  X  10 

6.8  X  10 


12 


12 


12 


12 


12 


9.4  X  10"^ 

-4 

1.1  X  10  ^ 

1.6  X  io"2 

1.2  X  10"^ 
4.9  X  10"^ 
1.2  X  10"^ 
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Table  XIV 

Activity  Concentrations  for  Contamination  from  Runoff 
for  the  Houston  Watershed 


! 

South  Wind 

West 

Wind 

1 

Isotope 

Atoms/liter 

4c/ml 

Atoms/liter 

uc/ml 

Sr-89 

4.5  X  10^^ 

1.9  X  10"^ 

3.6  X  10^° 

/■  "4 

1.6  X  10 

Sr-90 

8.8  X  10^^ 

1.9  X  10"*^ 

8.5  X  10^^ 

1.8  X  10"^ 

Ru-106 

4.1  X  10^ 

2.5  X  10“3 

4.0  X  10^^ 

2.4  X  10"5 

1-131 

12 

9.2  X  10 

2.5  X  10“^ 

i  8.4  X  10^° 

2.3  X  10*^ 

Cs-137 

12 

5.3  X  10 

-4 

1.0  X  10 

1,  4.3  X  10^° 

8.5  X  lO"'^ 

Ba-l40 

12 

11.3  X  10-^ 

1.9  X  10"^ 

11.0  X  .10^^ 

1 

1 _ 

1.9  X  10"2 

Table  XV 

Activity  Concentrations  for  Contamination  from  Runoff 
for  the  New  York  City  Watershed 


Isotope 

South  Wind 

1 - - 

West  Wind 

Atoms/liter  M-c/ml 

Atoms/liter  pc/ml 

Sr-89 

0.36  X  10^^  1.5  X  10’^ 

1.6  X  10^  6.8  X  10’3 

Sr~90 

0.65  X  10^^  1.4  X  lO"^ 

2.8  X  10^  5.9  X  10‘^ 

Ru-106 

12  -4 

C.29  X  10-^  1.7  X  10 

12  -4 

1.2  X  10  7.1  X  10 

1-131 

0.65  X  10^  1.8  X  10"^ 

2.8  X  10^  7.6  X  10'^ 

Cs-137 

0,48  X  10^^  9.5  X  10~^ 

2.2  X  10^  4.4  X  10‘5 

3a-l40 

0.8  X  10^^  1.4  X  10'^ 

3.4  X  10^^  5.8  X  lO'^ 
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C .  Oonclusloas ; 

It  is  not  the  purpose  here  to  draw  any  conclusions  in  regard  to 
the  'biological  hazard  resulting  from  the  concentrations  of  activity  presented 
in  ITahles  XI  -  XV.  Tiiis  aspect  is  being  considered  elsewhere  in  this  report 
and  also  by  other  investigators. 

It  may  be  prudent,  however,  to  consider  the  relative  value  of 
the  previous  calculations.  It  must  he  remembered  that  only  a  5  MT  weapon 
was  considered  and  that  all  concentrations  are  for  H  +  1  hour.  Thus 
for  any  other  size  weapon,  appropriate  revisions  in  the  fallout  model  will 
have  to  be  made  before  it  may  be  applied.  As  for  a  multi-bomb  attack, 
one  cannot  but  accept  some  method  based  on  additivity.  Complications 
arrive  here,  as  far  as  forming  a  model  is  concerned,  because  of  ground 
zero  locations.  Since  H  +  1  hour  values  of  concentration  are  given,  one 
may  apply  suitable  decay  curves  to  adapt  them  to  einy  specific  stream  flow, 
reservoir  flow-through,  or  intake  time,  etc.  for  a  particular  watershed 
system. 

It  is  not  to  be  implied  here  that  the  Houston  and  New  York  values 
are  representative  of  all  municipal  reser\'‘oirs .  However,  they  do  present 
a  reasonable  value  to  be  expected  in  a  realistic  case,  and  also  demonstrate 
the  ease  with  which  an  evaluation  may  be  made .  By  making  such  an  integrated 
evaluation  for  any  specific  watershed  and  applying  local  conditions,  the 
relative  radiological  hazard  may  be  e.stimated. 

The  Houston  and  New  York  values  will  be  seen  to  vary  from  the 
potential  values  if  unusual  ccnditions  prevailed  such  as  would  occur 
during  a  di^^r  spell.  Ilov^ever,  the  effects  will  be  somewhat  counterbalancing 


in  that  high  reservoir  concentration  due  to  low  volume  will  be  somewhat 
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offset  "bj  the  lack  of  runoff  contamination.  It  is  douhtecl  that  unnatural 
conditions,  vniless  very  severe,  will  change  the  concentration  of  activity 
hy  more  than  an  order  of  magnitude. 
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AHALYSIS  OF  RADIOACTIVITY  IN  WATER 

A.  Radiochemical  Methods 
1 ■  Introduction 

In  most  radiochemical  analyses  the  steps  are  concentration, 
separation,  purification,  and  counting.  Evaporation  and  ion  exchange 
ere  the  two  principal  methods  of  concentration,  hut  may  not  he  necessary 
if  the  level  of  activity  is  sufficiently  high.  The  individual  radio¬ 
nuclides  are  usually  separated  by  means  of  precipitation,  solvent 
extraction,  or  ion  exchange  elution.  Purification  of  the  individual 
radionuclides  usually  follows  standard  chemical  procedures  which  depend 
mainly  on  the  solubility  chareeteristics  of  compounds  of  the  radionuclides . 
Coryell  and  Sugaman  ^ ^  have  described  purification  of  many  of  the  more 
important  radionuclides.  Radioisotopes  are  counted  using  appropriate 
instruments,  depending  on  the  specific  isotope  and  the  desired  accuracy. 

"Carriers"  are  of  prime  importance  in  wet  radiochemical  procedures. 

The  mechanisms  and  usage  of  radioisotopic  carriers  have  been  well  defined 

( Cl) 

by  Overman  and  Clark  ^  A  carrier,  which  is  normally  the  stable  isotope 

of  the  element  being  determined  or  a  stable  element  with  very  similar 
properties,  is  added  to  increase  the  total  concentration  of  the  element, 
because  the  radioactivity  will  normally  be  present  in  micro  quantities  only. 
Hence,  if  the  carrier  is  not  added,  part  of  the  radioactive  material  will 
usually  be  lost,  even  though  employing  the  most  scrupulous  analytical 
techniques.  By  weighing  the  final  solid  to  he  counted,  and  knowing  the 
original  amount  of  carrier  added,  the  percentage  chemical  and  physical  loss 
may  be  calculated.  By  assuming  complete  interchange  of  carrier  and  radio¬ 
nuclide,  the  same  percentage  of  each  will  be  lost,  and  the  original  concen¬ 
tration  of  activity  may  readily  be  calculated.  The  carrier  itself  is 
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nomally  added  as  a  solution  of  a  soluble  salt  of  the  carrier  element 
prior  to  any  chemicai.  procedure.  Carrier  addition  is  used  extensively 
throughout  the  descrihed  analytical  procedures,  although  some  carrier-free 
determinations  are  now  coming  into  use . 

There  are  several  compilations  available  which  describe  detailed 

(20  22  23  2k-) 

procedures  for  the  quantitative  determination  of  many  radionuclides  ' 

Most  of  these  methods  have  been  proven  and  are  in  general  use. 

2.  Procedures  for  the  Determination  of  Specific  Radionuclides 
a.  Radio strontium 

Strontium-90  has  been  generally  recognized  as  being  a  very  bio¬ 
logically  hazardous  radionuclide  as  it  accumulates  in  bone .  Due  to  its 
relatively  dangerous  nature,  wide  interest  has  been  aroused  to  find  a  simple 
and  rapid  method  of  analysis  for  strontl\mi-90.  Techniques  that  have  been 
in  use  for  several  years  are  still  finding  widest  use  because  of  their 
accuracy,  the  time  consideration  being  less  important  in  peacetime . 

Strontium-89,  perhaps  somewhat  less  of  a  hazard  than  strontlum-90, 
because  of  its  shorter  half-life,  will  also  be  present  in  fallout  contaminated 
water.  Separate  analyses  of  strontlum-89  and  strontium-90  are  thus  desirable. 
Strontium-89  activity  is  normally  found  by  subtraction  of  strontlum-90  from 
total  radiostrontlum  concentration. 

(25) 

(1)  Radiostrontlum  and  Radiobarium  by  Ritrate  Separation^ 

Strontium  and  barivim  carriers  are  added  and  the  Group  II  cations 
ere  precipitated  as  the  carbonates.  Partial  separation  from  calcium  is 
accomplished  by  nitrate  precipitation  in  fiuning  nitric  acid.  The  remaining 
calcium  nitrate  is  then  extracted  with  acetone.  Rare  earths  and  other 
trlvalent  cations  are  removed  by  two  hydroxide  scavengings.  Barium  is 
finally  separated  from  strontium  by  precipitation  as  barium  chromate.  Stron¬ 
tium  is  then  collected  as  the  oxalate  and  counted.  The  barium  is  converted 
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from  the  chromate  to  the  chloride  for  counting.  Ihe  lowest  level  of  detection 

“8 

for  radiostrontium  hy  this  method  is  4  x  lO"  |j.c/ml  and  for  radioharlum, 

10  |jc/ml.  A  precision  of  about  10^  may  be  obtained. 

The  calculation  of  the  total  amount  of  radioactivity  due  to  either 
radiostrontium  or  radiobarium,  using  an  internal  proportional  counter,  is 
made  by  use  of  the  following  equation: 


net  cpm  , 

(A)(B)(c)tD)(2.22) 


where  A 
B 
C 
D 


efficiency  factor, 
chemical  yield, 
self-absorption  factor,  and 
volume  in  liters. 


This  formula  also  holds  xrue  for  determination  of  the  other  radionuclides 
A  variation  of  the  preceding  technique  was  Introduced  by  Kboi^^^^ 

After  the  separation  of  strontiiun  by  carbonate  and  nitrate  precipitations, 
barium  removal  is  effected  by  barium  chloride  precipitation  in  a  hydrochloric 
acid-ether  system.  Axij  lanthamam-l40  which  may  be  present  from  barium-l40 
decay  is  removed  by  a  ferric  hydroxide  scavenge.  Sensitivity  in  the  10  jic/ml 
range  is  obtainable  by  this  procedure . 

(2)  Strontium-90  by  Solvent  Extraction  of  Yttrium-gO^  ' 

Strontl\am  carrier  is  added  and  carbonate  precipitation  is  performed. 
The  carbonate  precipitate  is  dissolved  in  hydrochloric  acid  and  the  solution 
scavenged  with  hydroxide.  The  strontium  is  again  precipitated  as  the  carbonate, 
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which  is  then  allowed  to  stand,  thus  peimitting  yttrium-90  ingrowth.  The 
yttrium  is  extracted  with  tributyl  phosphate  3q.uilibrated  with  liiTJ  nitric 
acid.  Following  extraction,  the  yttrium  is  washed  from  the  organic  phase 
with  0.1  N  HNO^i  dried  in  a  planchet,  and  counted  in  an  anti-coincidence 
beta  counter. 

4 

A  decontamination  factor  of  approximately  10  is  obtained.  The 
chemical  yield  is  greater  than  70^  for  strontium;  the  recovery  of  yttrixm 
by  the  solvent  extraction  is  85^. 

Strontium- 90  determinations,  such  as  that  shown  above,  are  normally 
made  by  measuring  the  activity  of  its  beta  decay  daughter,  yttrium-90.  After 
a  suitable  period  has  elapsed,  the  yttrium-90  ingrowth  is  separated  from  the 
bulk  of  the  strontium  by  some  sultab].e  means.  The  yttrium-90  is  then  pujrified 
and  counted.  The  strontium  may  then  be  determined  by  calculation  using  the 
following  formula: 


strontium-go,  -  Xa)(b)(c)(dKe!(f)(2  .32)’ 


where : 


A  =5  .i.^’f iciency  factor, 

B  =  per  cent  extraction  factor, 

C  =  per  cent  ingrowth  factor, 

D  «  chemical  yield, 

E  *  sample  volume  (liters),  and 

F  =  decay  factor;  calculated  from 


A  = 


A 

o 


e 


-\t 


vhere : 
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A  =  activity  remaining  after  a  time  interval  t^ 

A^  =  activity  of  sample  at  some  original  time, 

^  _ 0.693 _ 

half -life' of  yttrim-90  (  62.4  hrs}  ' 
t  =  time  (hours)  from  separation  to  counting. 

By  knowing  the  total  stronti-jm-90  activity  and  the  total  activity  due 
to  radiostrontlumj  the  strontium-89  concentration  may  be  obtained  from  the 
difference . 

(3)  Strontium-90  by  Direct  Precipitation- of  Yttrium-90^^^^ 
Strontium-90  activity  may  also  be  easily  measured  by  a  variation  of 

the  preceding  method.  After  a  suitable  time  has  elapsed  for  yttrium-90  in¬ 
growth,  yttrium  carrier  is  added  and  subsequently  precipitated  as  the  hydroxide 
by  the  addition  of  sodium  hydroxide.  The  radioyttrium  is  then  purified  by 
standard  radiochemical  techniques,  which  will  depend  on  the  concentration  and 
species  of  the  other  rare  earth  radionuclides  present.  Resolution  of  the 
decay  curve  will  be  necessary  if  yttrium-91  is  present. 

(4)  Radiostrontium  by  Ion -Exchange  Methods 

Kahn  and  Reynolds^ utilized  ion-exchange  resins  for  the  concentra¬ 
tion  of  radiostrontium.  Strontium  carrier  is  added  to  a  10-llter  water  sample 
and  the  solution  is  passed  upflow  through  a  cation  exchange  resin  in  hydrogen 
form.  The  strontium  is  then  eluted  with  l^J-K  nitric  acid.  The  radio  strontium 

is  then  separated  from  other  alkaline  earths  and  purified  by  previously 

(26)(27)(29) 

mentioned  methods,'  '  '  depending  on  whether  total  radiostrontium  con¬ 

centration  is  desired  or  strontium-90  alone. 

By  means  of  concentration  of  activity  with  ion-exchange  resins,  great 
sensitivity  can  be  obtained,  possibly  down  to  the  10  |ic/ml  range. 
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A  similar  ion-exchange  method  was  described  by  Kahn,  Eastwood,  and 
Lacy ^31 ) ,  Using  appropriate  radiochemical  techniques,  other  long-lived  radio¬ 
nuclides  could  also  be  determined.  In  eveiy  case,  the  recovery  of  activity 
was  greater  than  99^- 

Bryant,  Sattlzahn  and  Warren  utilized  an  ion-exchange  procedure 
to  separate  yttrium  from  strontium  after  a  ten  day  ingrowth  period.  Tne 
yttrium  was  then  selectively  eluted  from  a  cation  exchange  resin,  followed  by 
a  readsorption  on  another  resin  and  finally  counted.  The  radiochemical  re¬ 
covery  was  greater  than  97^*  Gravimetric  measurement  is  not  required  when 
using  this  procedure. 


According  to  Libby  ’  strontium-90  may  readily  be  determined  in  low 
concentrations  by  concentrating  the  activity  by  evaporation. 

The  sample  (l  liter)  is  evaporated  to  dryness,  and  dry  strontium 
nitrate  carrier  is  added.  Following  addition  of  calcium  chloride,  to  the  dry 
residue,  the  sample  is  dissolved  in  phosphoric  acid.  Yttrium  is  then  removed 
by  the  addition  of  lanthanum  carrier,  followed  by  hydroxide  scavenging.  The 
precipitate  is  redlssolved  and  the  milking  repeated  twice  more.  The  third 
precipitate  is  saved,  ignited  and  counted.  The  analysis  itself  takes  two 
hours  to  perform.  It  is  necessary  to  utilize  decay  curves  to  determine  the 
yttrium-90  activity,  due  to  the  presence  of  yttrium-91  and  other  rare  earths. 

b .  Radiocesium 

Radiocesium  is  considered  to  be  one  of  the  more  important  of  the  bio¬ 
logically  hazardous  radionuclides.  Cesium  is  chemically  similar  to  the  other 
alkali  metals,  such  as  sodium  and  potassium,  which  are  commonly  found  in  the 
internal  organs.  Thus  ionic  interchange  within  the  body  is  easily  accomplished. 
Cs-137  and  Cs-134  with  half-lives  of  33  years  and  2.2  years  respectively 
constitute  the  greatest  radiocesium  danger. 
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Selective  alkali  precipitation  and  ion-exchange  elution  have  proven 
to  he  the  hest  methods  of  separation  for  the  analysis  of  radiocesium.  Deter¬ 
mination  of  the  various  isotopes  appears  to  he  most  feasible  hy  gamma  spectro¬ 
scopy  . 

(1)  Radioeeslum  hy  Phosphomolyhdate  Precipitation^ 

After  the  addition  of  cesium  carrier  to  the  sample,  phosphoric  acid 
and  ammonium  molybdate  are  added  and  the  cesium  is  precipitated  as  cesium 
ammonium  phosphomolyhdate.  After  dissolving  the  precipitate  in  sodium 
hydroxide,  the  cesium  is  reprecipitated  as  cesium  cohaltinitrite .  This 
precipitate  is  dried  at  100°C  on  an  aluminum  planchet  and  counted  in  an  in¬ 
ternal  proportional  counter.  A  93?^  removal  of  cesium  activity  was  attained 
hy  the  phosphomolyhdate  precipitation.  Concentrations  as  high  as  lO"^  nc/ml 
and  as  low  as  lO"^*^  ia.c/ml  were  measured  with  good  accuracy.  The  two  principal 
cesium  radioisotopes  measured  hy  this  method  were  cesium-137  and  cesiura-13^. 

(2)  Radiocesium  hy  Cohaltinitrite  Precipitation^ 

After  the  addition  of  cesium  carrier,  the  cesium  and  other  Group  I 
cations  were  precipitated  as  the  cohaltinitrite  after  the  addition  of  sodliun 
nitrite  and  cohaltous  chloride.  After  washing,  the  cohaltinitrite  is  dis¬ 
solved  in  hydrochloric  acid.  Cesium  silicotungstate  is  then  precipitated  hy 
the  addition  of  silicotungstic  acid.  The  silicotungstate  is  then  dissolved 
in  dilute  case  and  the  solution  scavenged  with  ferric  hydroxide.  The  excess 
tungsten  is  removed  as  the  insoluble  trioxide.  Cesium  and  sodium  perchlorates 
are  then  precipitated  hy  the  addition  of  perchloric  acid  and  absolute  alcohol 
to  the  solution.  The  sodium  perchlorate  is  then  removed  hy  washing  the  pre¬ 
cipitate  with  absolute  alcohol.  The  final  cesium  precipitate  is  then  washed, 
dried  and  coimted.  This  method  has  been  proven  to  he  quite  accurate  and  is 


now  in  wide  use. 
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Oamond,  et  al.^  co-precipitated  cesium  cobalt inltrite  on  potassiom 
cobaltinitrite  and  followed  with  a  similar  purification  procedure.  This 
method  has  found  little  popular  use. 

(3)  Radiocesj-um  hy  Silicotimgstate  Precipitation^ ^ 

This  procedure  is  based  on  a  method  originally  described  by  Yamagata 
and  ^amagata'"^  ,  The  cesl\am  is  first  precipitated  as  the  sillcotxmgstate . 
Dissolution  of  the  precipitate  is  then  followed  by  a  ferric  hydroxide 
scavenging.  The  cesium  is  then  reprecipitated  as  the  dlpicrylaminate .  This 
salt  is  dissolved  in  4 -methyl-2 -pentanone  and  the  cesium  is  extracted  by  means 
of  2M  hydrochloric  acid.  The  cesium  is  finally  pi^cipltated  as  the  perchlorate, 
in  which  form  it  is  dried,  weighed,  and  counted.  The  chemical  yield  is  80^,  and 
eight  analyses  may  be  performed  in  eight  hours. 

(23) 

(4)  Radiocesium  by  a  Co-crystallization  Procedure^ 

A  method  has  been  described  by  which  cesi'um  is  separated  from  the  bulk 
of  the  alkali  elements  and  the  mixed  fission  products  by  co-crystallization 
with  ammonium  aluminum  sulfate.  The  ammonium  salts  are  decomposed  by  heating, 
and  the  cesium  is  precipitated  for  counting  from  a  dilute  hydrochloric  acid 
solution  as  the  chloroplatinate .  This  method  has  been  in  use  as  standard 
procedure  in  sc*ue  laboratories. 

(5)  Radiocesium  by  Ion-Exchange  Methods 

(31) 

Kahn,  Eastwood,  and  Lacy  have  developed  a  separation  scheme  for 
the  analysis  of  the  more  hazardous  radionuclides  by  ion  exchange .  The  specific 
radioelements  studied  were  cesium,  cerium,  cobalt  and  strontium.  The  lower 
limit  of  detection  was  decreased  a  hundredfold  by  concentrating  the  radio¬ 
activity  of  a  large  sample  with  a  cation  exchange  resin. 

Cesium  was  selectively  eluted  by  means  of  04  hydrochloric  acid  from 
the  cation  exchange  resins  studied.  Standard  purification  procedures  were 
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then  used  to  prepare  the  sample  for  counting.  The  cesium  was  separated  from 
the  other  alkali  cation  contaminants  hy  precipitating  it  as  the  silicotung- 
state,  dissolving  the  precipitate  in  sodium  hydroxide,  and  then  reprecipitating 
it  as  the  perchlorate  for  counting.  The  activity  recovered  for  all  radio¬ 
nuclides  was  greater  than 

(37) 

Tsuhota  and  Kitano  found  that  an  ammonium  formate  -  formic  acid 
buffer  of  pH  3.2  selectively  eluted  the  alkali  metals  from  a  cation  exchange 
resin.  This  buffer  was  employed  particularly  for  the  determination  of  radio¬ 
cesium,  being  somewhat  superior  to  hydrochloric  acid  or  citrate  buffer  as 
an  eluant. 

c .  Radio iodine 

Radio iodine  may  be  either  beta  or  gamma  counted,  although  gamma 
scintillation  is  preferred  for  determination  of  the  individual  radioisotopes 
of  iodine. 


(1)  Radiolodtne  by  Chemical  Methods 

Glendenin  and  Metcalf  determined  radioiodine  activity  by  an 
extraction  purification  procedure.  Carrier  sodium  iodide  is  added  to  the  water 
sample,  and  interchange  is  accomplished  by  oxidation  to  the  iodate  with  sodium 
hypochlorite  in  basic  solution,  followed  by  reduction  to  the  iodide  by  sodium 
bisulfite  in  acid  solution.  Sodium  nitrite  is  then  added  to  oxidize  the  iodide 
to  elemental  iodine  which  is  then  extracted  into  carbon  tetrachloride .  The 
iodine  is  further  purified  and  concentrated  by  back-extraction  into  sodium 
bisulfite  solution  which  is  finally  gamma  counted  at  the  iodine-131  photopeak. 

g 

The  sensitivity  of  the  method  is  approximately  lO"  c/ml  for  a  liter  sample. 

The  problem  of  incomplete  carrier  interchange  is  overcome  by  this 
procedure  which  employs  oxidation-reduction. 
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The  same  authors  '  also  describe  a  procedure  in  which  the  iodate  is 
reduced  directly  to  the  iodide  by  hydroxylamine  hydrochloride.  The  iodide 
obtained  after  the  sodium  bisumfite  back-extraction  is  precipitated  with 
silver  nitrate  as  silver  iodide  which  is  then  dried,  weighed,  and  beta 
counted . 


A  similar  but  faster  method  _ ..  radioiodine  determination  was  de- 

( 39) 

veloped  by  Lewis  He  developed  a  continuous  extractor  employing  two 

centrifugal  pumps,  to  promote  mixing  and  extraction.  The  iodine  was  extracted 
into  carbon  tetrachloride  and  backextracted  into  bisulfite  solution  in  one 
operation.  The  gamma  emission  of  a  pipetted  sample  was  then  measured.  The 
time  of  analysis  was  cut  from  2  hours  to  30  minutes.  An  accuracy  of 
93*0  -  8.0^  was  achieved. 

d.  Total  Radio  Rare  Earths  Determination 

(1)  Total  Rare  Earth  Activity  by  Fluoride  Precipitation 

Hume  and  Martens  ^^^^ave  developed  a  rapid  method  for  the  determina¬ 
tion  of  rare  earth  activity  by  which  the  beta  and  gamma  activities  are  deter¬ 
mined  separately.  The  method  is  quite  rapid,  the  gamma  determination  taking 
only  a  half-hour. 

Two  rare  earth  fluoride  precipitations  serve  to  remove  zirconium  and 
niobium,  the  principle  gamma  emitters  of  fission  material.  Barium  and  stron¬ 
tium,  which  are  heavily  coprecipitated,  are  removed  by  hydroxide  precipitation 
after  the  addition  of  holdback  carrier.  The  sample  may  then  be  gamma  counted. 
To  estimate  the  beta  activity,  the  rare  earths  must  be  precipitated  as  the 
oxalates  before  counting. 

(2)  Total  Rare  Earth  Activity  by  an  Alternate  Method 

(4l) 

Boldridge  and  Hume  developed  a  method  similar  to  that  above 
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for  the  detemination  of  total  rare  earth  activity.  It  is  somewhat  lengthier, 
hut  more  accurate,  particularly  for  beta  counting. 

Stable  cerium  is  added  as  a  carrier  for  the  entire  group  of  rare 

earths.  The  rare  earths  are  then  precipitated  as  fluorides,  redissolved  in 

boric  and  nitric  acids,  and  reprecipitated  as  the  hydroxides.  Zirconium  is 

then  removed  by  precipitation  as  the  iodate  after  reduction  of  all  the  rare 

earths  to  the  trivalent  state .  Following  hydroxide  scavenging  for  alkaline 

earth  removal,  the  rare  earths  are  precipitated  as  the  oxalates  and  beta- 

counted  in  an  internal  proportional  counter. 

(42) 

Alstad  and  Pappas '  '  employed  a  similar  technique,  using  lanthan\aa 

as  a  carrier  for  all  the  rare  earths .  Following  separation  and  purification 
of  the  group  of  rare  earths,  individual,  separation  was  attained  by  elution 
from  a  cation  exchange  column  with  ammonium  lactate .  The  ion-exchange  separa¬ 
tion  step  was  performed  without  additional  carrier,  and  the  recovery  of 
activity  for  each  radionuclide  was  assumed  to  be  the  same  as  that  for  the 
lanthanum.  The  advantage  of  the  non-carrier  separation  is  that  a  "weightless" 
sample  is  obtained,  thus  eliminating  self-absorption  and  self-scattering 
factors . 

e .  Miscellaneous  Radionuclides 

(l)  Radiobarium 

Determination  of  radiobarium  may  easily  be  made  during  the  radio¬ 
strontium  analysis.  Both  strontium  and  barium  form  insoluble  carbonates,  and 
subsequent  separation  may  be  made  by  precipitation  of  barium  as  the  chromate. 

(25) 

This  method,  by  Hahn  and  Straub,  was  described  earlier  under  radio strontium. 

(45) 

Hunter  and  Perkins  followed  the  same  procedure  except  that  the  precipitate 
is  counted  as  the  carbonate  instead  of  the  chloride . 
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(  2k') 

Minltklnen  '  has  developed  a  procedure  for  the  analysis  of  radio- 
harium  in  which  harium  chloride  carrier  is  specifically  precipitated  as  the 
monohyd:'ate  from  a  c  .ncentrated  hydrochloric  acid-ethyl  ether  mixture.  This 
precipitate  is  dissolved  in  water  to  which  ferric  hydroxide  is  su'.sequently 
added  as  a  scavenger.  The  harium  is  then  isolated  as  the  chromate.  The 
chemical  yield  of  harium  chromate  is  about  Tiie  harium  chromate  precipi¬ 

tate  is  then  set  aside  for  134  hours  to  allow  the  harium- l4o  and  its  daughter 
lan.thanum-l40  to  come  to  equilibrium. 

(2)  Radiocerium 

( 

B'urgus  and  Engelkemelr  ‘  have  developed  t-,  method  of  analysis  for 
radiocerltim  ( cerium- iWt- ) ,  employing  a  series  of  precipitations.  After  the 
addition  of  cerium  carrier  and  alkaline  earth  h<.ildhack  carriers,  ceri-um  was 
separated  along  with  the  rr.re  earths  hy  ■/luoriuc  p-^ecipitation.  Alternate 
hydroxide  and  fluoride  precipitations  are  then  employed  to  remove  zirconium 
and  .'5l.k.iline  earths  completely.  The  cerium  is  then  oxidized  to  the  tetravalent 
state  with  perchloric  acid  and  precipitated  as  the  iodate.  Thorium  activity  is 
then  removed  as  the  iodate  following  the  reduction  of  cerium  with  sulfur 
dioxide,  lifter  scavenging  with  thorium  and  f'urther  reprecipitation,  the  cerium 
is  finally  precipitated  as  cerous  oxalate,  which  is  subsequently  ignited  to 
the  oxide.  After  weighing,  the  cerium  activity  is  counted  on  an  internal 
proportional  counter. 

The  percentage  of  chemical  recovery  is  only  about  30^,  i>ut  the  radio¬ 
active  purity  of  the  final  precipitate  is  very  high. 

It  is  likely  that  this  process  could  he  shortened  because  of  the 

probable  absence  of  thorium  from  fallout  contaminated  water. 

(4s) 

Ames  '  employed  a  similar  procedure  for  the  analysis  of  total 

radiocerium.  The  concentration  of  cerium-l44  is  obtained  by  counting  im- 

2 

mediately,  using  a  217  mg/cm  aluminum  absorber.  This  shields  all  the 
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cerium  lietas,  cotonting  only  those  from  praseodymi\jm-l44,  the  daughter  product 
of  cerium-l44.  The  difference  hetween  the  total  radiocerium  concentration  and 
cerium-l44  concentration  is  assumed  due  to  cerium-l4l,  a  neutron  induction 
product . 

(  C’V 

Barnes  describes  a  method  hy  which  radiocerium  is  isolated  hy  an 

extraction  nKthod  and  counted  for  cerium-l44.  Carrier  interchange  is  accom¬ 
plished  hy  means  of  an  oxidation-reduction  cycle.  Two  dibutyl  phosphate  ex¬ 
tractions  are  then  carried  out  to  remove  the  heavy  elements.  The  cerium  is 
then  precipitated  as  the  fluoride,  oxidized  to  the  tetravalent  state,  and  ex¬ 
tracted  into  hexone.  Finally  it  is  converted  to  the  dioxide  and  counted  for 
cerium- l44  hy  means  of  a  beta  coianter  using  a  217  mg/cm  aluminum  absorber. 

The  chemical  yield  varies  from  50-6o^. 

Kahn  and  Reynolds^  have  developed  another  method  for  the  determina¬ 
tion  of  radiocerium  activity  employing  the  use  of  ion-exchange  resins.  Con¬ 
centration  of  radiocerium  and  added  carrier  was  attained  by  passing  a  one-liter 
water  sample  through  either  Dowex-50  on  the  hydrogen  cycTe  or  IR-220  on  the 
sodium  cycle.  Ceriun  was  selectively  eluted  from  the  resin  with  25  ml  of 
3M  nitric  acid.  Purification  of  the  radiocerium  is  obtained  by  successive 
precipitations  as  the  fluoride,  the  hydroxide  and  finally  as  the  oxalate.  The 

cerium  oxalate  is  weired  to  determine  carrier  loss  and  then  counted. 

-8 

Sensitivity  of  detection  was  increased  to  approximately  lO”  uc/ml. 

(3)  Radloruthenlum 
(  4t  j 

Glendenin'  ‘  has  developed  an  analysis  scheme  for  the  determination 
of  radioruthenium  which  has  been  utilized  as  a  standard  method  of  procedure. 
Separation  of  the  radioruthenium  with  carrier  ruthenium  is  accomplished  by 
oxidation  with  perchloric  acid  to  the  volatile  tetraoxlde  which  is  subsequent¬ 
ly  distilled  off.  Sodium  bismuthate  is  added  prior  to  distillation  to  prevent 


-118- 


volatllization  of  the  halides  by  oxidizing  them  to  the  oryacids .  !Eie  ruthenium 
is  absorbed  into  a  sodi\am  hydroxide  solution  and  precipitated  in  the  form  of  its 
lower  oxides  by  reduction  with  ethanol.  The  ruthenium  oxides  are  tnen  dis¬ 
solved  in  hydrochloric  acid^  and  ruthenium  is  precipitated  in  the  metallic 
state  by  reduction  with  magnesium  metal.  The  ruthenium  metal  is  weighed  for 
carrier  loss  and  counted.  The  chemical  yield  is  about  65^. 

Melnick^^^^  ,  using  an  almost  identical  chemical  procedure,  describes 
a  process  for  the  determination  of  rathenium-106.  By  using  a  series  of  heavy 
aluminum  absorbers  with  a  beta  proportional  counter,  a  correction  factor  for 
the  presence  of  ruthenium- IO3  (half-life  =  h2d)  may  be  obtained  by  extrapola¬ 
tion.  Extrapolation  for  the  correction  is  necessary,  in  that  while  counting 
the  betas  from  ruthenium- IO6  throiagh  a  210  mg/cm  aluminum  absorber,  the  gamma 

rays  from  ruthenium-103  are  also  counted.  If  the  activity  due  to  ruthenium-103 

2 

alone  is  desired,  gamma  scintillation  may  be  used  with  a  2,000  mg/cm  aluminum 
absorber. 

Merritt^^^^  describes  a  method  of  radioruthenium  analysis  in  which  the 
ruthenium  is  determined  in  the  presence  of  strontium,  cesium  and  cerium  by 
fusion  and  extraction.  Following  the  additJ.on  of  carrier  ruthenium  to  the 
water  sample,  the  sample  is  fusee,  with  potassium  hydroxide,  sodium  nitrite, 
and  sodium  carbonate  at  550°C  for  two  hours .  The  melt  is  leached  twice 
with  water  to  dissolve  the  ruthenate  and  cesium,  leaving  strontium  and  cerium 
in  the  residue .  Strontium,  cerium  and  cesium  are  then  analyzed  by  previously 
mentioned  procedures.  Tlie  ruthenium  is  extracted  from  the  leach,  8,fter  the 
addition  of  periodate,  with  carbon  tetrachloride.  The  ruthenium  is  then 
stripped  from  the  carbon  tetrachloride  with  6N  hydrochloric  acid.  Magnesium 
is  then  used  to  reduce  the  ruthenium  to  its  elemental  metallic  state  in  which 
form  i  is  subsequently  counted.  This  .method  has  proven  satisfactory  for 
soi.'’  leaching  analysis,  etc.,  but  could  well  be  used  on  water  samples. 
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(4)  Radlozlrconium-Radloniotalum 

Steinberg has  developed  a  method  for  the  determination  of  zirconium- 
niobium  activity  which  has  been  tested  and  proven.  Following  the  addition  of 
both  zirconium  and  niobitm  carriers,  the  two  elements  are  complexed  by  the  ad¬ 
dition  of  oxalic  acid.  Thorium  oxalate  scavenging  is  carried  out  by  the  ad¬ 
dition  of  potassium  chlorate  to  the  acidified  solution.  The  niobium  oxide 
is  redissolved  and  reprecipitated  as  the  oxide  h;9'<it’ate .  The  water  is  then 
driven  off  by  ignition  and  the  niobi'um  oxide  is  weighed  and  mounted  for  count¬ 
ing.  The  zirconium  is  precipitated  from  the  original  solutj  ---i  as  the  phen- 
ylarsonate  by  the  addition  of  phenylarsonic  acid.  This  pres  pj  ite  is  dissolved 
in  oxalic  acid,  forming  the  soluble  oxalate.  The  zirconium  is  reprecipitated 
as  the  phenylarsonate ,  which  is  subsequently  ignited  to  zirconium  dioxide. 

This  precipitate  is  weighed  and  counted.  Both  precipitates  (niobium  and 
zirconium)  are  beta  counted  on  an  internal  proportional  counter. 

Brady  and  Engelkeraeir have  developed  a  phosphate  method  for  the 
determination  of  activity  due  to  zirconium-niobium  which  is  somewhat  lengthier 
than  the  preceding  method.  Separation  of  zirconium  is  based  upon  the  precipi¬ 
tation  of  zirconium  phosphate .  The  phosphate  is  then  redissolved  as  the 
fluoride  and  precipitated  as  the  phenylarsonate  which  is  subsequently  ignited 
for  counting.  Niobium  is  precipitated  in  acid  solution,  dissolved  as  the 
fluoride  complex  and  finally  weighed  and  counted  as  the  oxide.  About  10  hours 
is  required  for  the  entire  analysis .  The  chemical  recovery  of  zirconium  is 
about  and  that  of  niobium  is  about 

( 52) 

A  method  is  also  described  by  Stanley  for  determination  of 

zlrconiuHi-95 '  Carrier  interchange  is  effected  by  the  formation  of  the  fluoro- 

-2 

zirconate  complex,  ZrFg  •  Rare  earth  activities  are  removed  by  lanthanum 
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fluorlde  scavenging.  The  zirconium  is  then  separated  hy  means  of  three  'barium  ^ 

fluorozirconate  precipitations.  The  zirconium  is  finally  precipitated  with 
mandelic  acid  from  hydrochloric  acid  solution  and  ignited  to  the  dioxide,  in 
which  form  it  is  weighed  and  coixnted.  A  chemical  yield  of  about  751^  is 
obtained.  It  is  essential  to  not  begin  the  analysis  until  zirconium-97 

°  hours)  has  had  time  to  decay  to  a  negligible  amount.  After  the 
final  precipitate  has  been  obtained,  counting  should  begin  immediately  so 
that  no  appreciable  nlobium-95  has  time  to  grow  in.  A  beta  proportional 
counter  should  be  employed  for  the  counting. 

3.  Summary  of  Radionuclide  Analyses 

For  the  convenience  of  the  reader,  additional  references,  as  well 
as  those  previously  cited,  are  summarized  in  Table  XVI. 
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Table  XVI 

SuTnmary  of  Refereaces  of  Radionuclide  Analyses 


Element 

References  Cited 

Additional 

References 

1. 

Radiost rent ium 

25, 

26, 

27,  29, 

53, 

54, 

55, 

30, 

31, 

32,  33 

56, 

57, 

61, 

62 

2. 

Radioceslum 

23, 

31, 

34,  35, 

29, 

^9, 

54, 

55, 

36, 

37 

58, 

64, 

65 

3. 

Radioiodine 

28, 

38, 

39 

55, 

59 

k. 

Total  Radio  Rare 
Earths 

ko, 

^1, 

42 

55, 

6o 

5. 

Radiobarlum 

2k, 

25, 

43 

29 

6. 

Radiocerium 

30, 

lA, 

45,  46 

29, 

31, 

^9, 

63 

7- 

Radloruthenlum 

^7, 

48, 

49 

30 

8. 

Radiozirconium  - 
Radionloblum 

50, 

51, 

52 

30 

-122- 


B .  Instrumentation 

Quantitative  analysis  of  radioactive  elements  in  trace  concen¬ 
trations  is  done  most  conveniently  by  meeuas  of  radiation-detecting  instru¬ 
ments.  Due  to  the  veiy  small  amounts  of  radionuclides  present,  such 
instrments  are  usually  superior  in  speed  and  sensitivity  to  conventional 
wet  chemical  analyses.  However,  when  complete  chemical  analyses  are 
required,  extensive  use  is  made  of  wet  chemical  procedures  for  the  separa¬ 
tion  of  various  isotopes,  as  described  on  pages  105  ~  121.  Such  separations 
are  not  required  when  only  gross  fission  product  activity  is  to  he  mea¬ 
sured  or  when  only  one  radioisotope  is  present. 

If  the  activity  of  a  sample  is  sufficiently  high,  this  activity 
may  he  determined  directly  using  a  simple  detector  and  ratemeter.  At 
veiy  low  concentrations  in  aqueous  samples  or  for  beta  emitters,  for  which 
self-absorption  in  the  sample  is  a  major  factor,  preconcentration  of  the 
sample  and  specially  designed  instruments  such  as  a  flow  counter  may  be 
required . 

A  survey  of  commercially  available  portable  radiation  instruments 
has  been  carried  out  and  the  information  obtained  is  summarized  in 
Tables  XvIH  and  XIX.  These  are  all  instruments  of  the  ratemeter  type 
capable  of  recording  count  rates  down  to  a  few  co’,mts  per  minute.  To 
establish  their  usefulness  for  the  deteimination  of  fallout  product 
concentrations  in  water  it  is  necessary  to  correlate  the  sensitivity  of 
the  equipment  with  different  types  of  detectors  and  reasonable  sample 
volumes  with  the  level  of  contamination  concentration  expected  in  emergencies. 

For  raw,  untreated  aqueous  solutions,  Figure  15  shows  the  relation 
of  nuclide  concentration  to  measured  dose  rate  for  various  isotopes  of 
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inteie&t  for  a.  representatiYe  commercial  instr.^eat .  It  has  heen  shown  in 

Tshles  A..*'  and  Xf  that  fallotit  product  concentrations  of  the  order  of 

-  ^  -L  , 

10  -  10  may  he  expected  in  the  water  supplies  of  typical  v^atersneds 

during  the  early  post -attack  period.  Table  X' .  L  taken  from  Lacy  and 
Kahn'  '■  shows  the  relative  detectability  of  selected  radioisotopes  with 
survey  meters  w'lth  a  lovest-ecale  sensitivity  down  to  0.01  mr/hr  ass'aming 
ideal  counter  geometry. 

If  it  is  merely  desired  to  determine  the  presence  or  absence  of 
these  high-level  lallout  concentrations,  it  is  evident  that  available 
survey  instiuraents  sho'old  be  capable  of  suppjyiiig  this  information,  pro¬ 
vided  thau  they  are  fitted  with  a  detector  head  suitable  ei.the.r  for  Immersion 
into  the  liquid  sample  or  for  flow-through  measurements. 

Movever,  in  most  cases  the  concentrations  measured  6.re  well 

(  67  ^ 

above  the  :1CLP  tnaxim’ora  permissible  levels  of  concentration  in  water’  '  or 
close  tc  the  limits  of  detectability  of  the  instr-jments  listed  in  I'ahle  XV aj..!. 
Por  this  reason  it  seems  necessary  to  develop  and  make  available  procedures 
and  facilities  for  the  more  accurate  measurement  of  low-Level  fission  product 
concent ra.’t ions  to  ensure  the  safety  of  available  sources  of  drinking  water. 
This  development  can  follow  twc  paths: 

1.  Pxeconcentration  of  the  sample  to  increase  the  number  of 

radioactive  atoms  to  a  level  that  is  readily  detected  by  conventional 

>• 

survey  instruments;  or 

2.  Increasing  the  sensitivity  of  the  detectors  by  an  increase 
in  sample  ''olume,  background  teduction  by  added  shielding  and  coincidence 
circuit  arrangements,  or  more  intimate  contact  between  sample  and  detector 
volumes,  aF  rn  proportional  flow  co-mters  or  in  liquid  scir.tiiiation 


counting. 
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Tai;i  e  XvTI 
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preconcentration  Methods 

Many  methods  of  preconcentration  have  disadvantages  vhich  limit 

their  use  in  an  emergency.  Concentratior  hy  evaporation,  for  example, 

requires  a  great  deal  of  care  and  skill  to  obtain  reasonably  accurate 

results.  The  use  of  ion  exchange  resins  to  remove  the  activity  from  the 

sample  can  give  very  good  reproduceability  and  in  some  applications  can 

(68) 

be  almost  completely  automated.  Emmons  sind  Lauderdale  have  developed 
a  continuous  water  monitor  which  employs  a  long  G.M.  tube  at  the  center  of 
a  column  of  cation  exchange  resin.  By  employing  two  columns,  one  may  be 
used  while  the  other  is  oeing  regenerated  with  strong  acid.  With  this 
system  sensitivities  of  around  lO”'^  pc/ml  have  been  obtained. 

The  mixed  resin  bed  is  potentially  capable  of  sensitivities  in 
the  lO”  '  -  lO”  |i'c/ml  range  when  used  with  low  background  counters. 

Although  the  measurement  of  radioactivity  on  resin  presents  problems  of 
self  absorption  in  the  sample,  reasonable  efficiencies  and  reproduceabilitles 
can  be  obtained  if  sufficient  care  is  taken  in  sample  preparation. 

One  method  which  may  be  employed  to  overcome  self-absorption 
in  the  resin  is  to  detect  the  beta  particles  with  a  scintillation  resin. 
Little  work  has  been  done  on  scintillation  resins,  but  when  techniques 
are  developed  this  should  become  a  very  valuable  tool  in  contamination 
evaluation  work. 

A  highly  accurate  way  in  which  the  radioisotope  content  of  water 
can  be  determined  is  by  radiochemical  analysis .  By  chemically  separating 
the  various  elements  of  interest  and  then  measuring  the  activity  of  each 
constituent,  the  activity  concentration  as  well  as  the  identity  of  the 
radio  contaminants  may  be  obtained.  This  type  of  analysis  requires  a 
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skllled  radiochemist  and  such  an  individual  may  be  difficult  to  find  in 
emergencies.  Solvent  extraction  methods  have  also  been  ^lsed,  but  at 
present  a  method  involving  a  mixed -bed  ion-exchange  resin  seems  to  offer 
the  simplest  practical  approach^ . 

Detector  Development 

Improvement  in  detector  sensitivity  can  be  achieved  by  a  suitable 
choice  of  detector  characteristics,  detector  geometry,  and  electronic 
circuit  parameters.  In  order  to  measure  the  contamination  level  of  water 
supplies  the  radiation  emitted  from  the  water  sample  must  be  determined 
as  efficiently  and  as  accurately  as  possible.  Because  of  the  statistical 
nature  of  the  radiation  emission,  higher  accuracies  can  be  obtained  only 
by  increasing  the  total  nmber  of  counts  registered  by  the  detector, 
while  at  the  same  time  keeping  background  counts  to  a  minimum.  For  this 
reason,  at  low  levels  of  activity  individual  counting  circuits  or  scalers 
are  superior  to  integrating  counting  circuits  such  as  ratemeters. 

If  the  radioactivity  is  composed  largely  of  beta  particles,  a 
major  problem  in  measuring  beta  activity  in  water  is  the  short  distance  which 
a  beta  particle  can  travel  in  this  dense  medium.  It  is,  therefore,  desirable 
to  bring  as  large  a  volume  of  water  as  possible  into  the  immediate  proximity 
of  the  detector. 

One  of  the  most  convenient  monitors  for  this  purpose  is  the  flow- 
throtigh  type  which  allows  the  water  to  come  in  contact  with  the  detector  for 
a  brief  period  and  then  be  flushed  out .  The  most  common  detector  used  in 
flow-through  monitors  is  the  G.M.  tube^'^^\  Various  configurations  of  the 
G.M.  monitor  type  have  been  used  to  check  cooling  water  effluent  at  infeuay 
of  the  water-cooled  reactors  in  the  co\mtry. 
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Lately,  the  geometrical  flexlhility  of  plastic  scintillators 

has  increased  the  volume  of  vater  which  can  be  brought  close  to  the  active 

volume  of  the  detector.  There  are  presently  on  the  market,  water  monitors 

which  use  large,  spherical  scintillators  capable  of  measuring  activities  in 

the  range  of  10  |ic/ml^'^^\  Another  method  which  displays  much  of  the 

geometrical  advantage  of  liquid  scintillation  but  has  none  of  the  chemical 

problems  is  the  use  of  scintillation  fibers'  .  Although  this  method 

lij- 

has  been  used  for  the  most  part  to  measure  C  and  tritium  in  water,  it 
should  yield  good  results  with  fission  products  and  extend  direct  monitoring 
capabilities  below  the  10  |ic/ml  range. 

Liquid  scintillation  techniques  have  yielded  very  accirrate  results, 
hi^  efficiencies  and  good  sensitivity  in  the  determination  of  radioactivity 
in  water.  This  method  may  prove  quite  valxiable  in  the  hands  of  an  experienced 
radiochemist . 

The  limiting  factor  in  the  sensitivity  of  a  detection  device 
is  the  background  or  "zero  activity"  response  of  the  device.  Both  electrical 
and  physical  approaches  have  been  pursued  in  the  reduction  of  this  background. 
Physical  methods  include  the  shielding  of  detectors  and  the  reduction  of 
detector  sensitivity  to  the  gamma  radiation  which  can  penetrate  the  shielding. 
Since  the  cross  section  of  plastic  scintillators  is  quite  high  for  beta 
radiation  while  being  fairly  low  for  gamma  radiation,  very  thin  scintillation 
discs  have  been  used  to  produce  good  beta  detection  efficiences  with  a 
background  of  less  than  one  half  count  per  minute  for  a  one  inch  diameter 
sample 

The  electrical  method  of  reducing  background  involves  a  shielding 


detector  and  anti-coincidence  circuits.  The  detection  device  is  thus  made 
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insensitive  to  most  of  the  hackgroiind  radiation  which  enters  the  detector. 
Most  of  the  anti-coincidence  detectors  presently  on  the  market  reduce  the 
apparent  hackground  to  a  range  of  0.5  to  2  counts  per  minute. 

The  only  rellahle  method  presently  applicable  to  civil  defense 
use  for  measuring  the  amount  of  specific  isotopes  in  water  is  gamma-energy 
analysis.  Most  gamma  analyzers  are  large  and  quite  complicated  to  opera±e, 
although  some  that  are  portable  and  easy  to  use  are  now  appearing  on  the 

(74) 

market.  With  a  small  amount  of  training  civil  defense  personnel  should 

he  able  to  make  good  use  of  these  instruments. 

Summary 

Preparation  for  the  post  attack  monitoring  of  drinking  water 
supplies  for  radioactive  contaminants,  may  now  be  divided  into  four  areas 
of  required  action.  These  include  the  installation  of  continuous  monitoring 
equiiaient  at  treatment  plants,  dissemination  of  information  on  the  use  of 
available  instruments  for  water  activity,  insuring  the  availability  of 
battery-powered  equipment  designed  for  water  or  general  liquid  radioactivity 
assay,  and  continuing  the  research  on  equipment  capable  of  meeting  emergency 
requirements . 

A  large  number  of  water  treatment  plants  would  probably  still  be 
operational  after  a  nuclear  attack.  The  danger  of  contamination  from 
fallout  would  require  that  the  water  be  continuously  monitored  for  contamina¬ 
tion  even  at  levels  well  below  the  maximum  permissible  levels.  In  this  way 
alternative  water  sources  could  be  located  or  decontamination  procedures 
could  be  brought  iru.j  operation  as  soon  as  the  possibility  arises  that 
existing  water  supplies  become  unfit  for  consumption.  It  would  be  desirable 
for  this  permanent  equipment  to  be  able  to  determine  the  relative  concentra¬ 
tions  of  the  various  radiocontaminants. 
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Many  commercial  instruments  which  detect  nuclear  radiations  can 
he  set  up  to  measure  the  amount  of  radioactivity  in  water  to  a  degree  of 
accuracy  suitable  for  emergency  use.  This  measurement  would  require  a 
given  procedure  and  a  set  of  conversion  tables  for  each  type  of  instrument 
or  class  of  instruments. 

Of  the  widely  available  instruments,  the  most  common  type,  the 
G.M.  survey  meter,  can  be  adapted  for  emergency  use  by  either  dipping  the 
probe  into  the  water  or  holding  the  probe  close  to  the  water  surface,  to 
measure  radio-contamination  levels  down  to  lO”^  jic/ml  of  some  of  the 
isotopes.  Some  scintillation  meters  can  reduce  this  lower  limit  below 
10  ^  |j,c/ml  with  the  absolute  limit  depending  on  background  radiation. 

Most  ion  chamber  survey  meters  are  not  as  sensitive  as  the  G.M.  or  scintil¬ 
lation  types,  although  some  of  the  more  sensitive  ones  might  be  used  when 
no  other  equipment  is  available. 

Companies  marketing  nuclear  detection  equipment  were  asked  to 
supply  detailed  specifications  on  their  survey-type  instruments.  A 
summary  of  information  received  by  January  30,  19^3  is  given  in  Table  XV" II. 
The  instruments  are  listed  in  order  of  decreasing  sensitivity  (minimum 
full  scale  readings ) , 

Pertinent  specifications  in  the  table  include  the  maximum  full 
scale  range,  battery  life,  type  of  radiation  detected,  type  of  detector 
used,  weight,  the  manufacturer's  model  number,  and  unit  cost  of  the 
instrument .  The  ability  of  an  instrument  to  discriminate  between  various 
types  of  radiation  is  also  included,  as  vrell  as  the  type  of  radiation  to 
which  a  given  instrument  is  sensitive,  is  indicated  by  an  "X"  in  the 
appropriate  column.  If  an  "o"  appears  under  the  indicated  type  of  radiation 
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XABI£  XXII  Capabilities  of  Portable  Survey  lustnitoents 
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it  simply  means  that  it  can  be  shielded  out  by  a  movable  absorber.  In  the 
event  an  instrument  may  be  powered  by  standard  "d"  cells,  the  number  of 
batteries  required  is  Indicated.  Instruments  which  are  specified  to  be 
moistiire  proof  or  immersion  proof  are  indicated  to  be  "weather  proof"  by 
an  "X"  in  the  appropriate  colxomn.  It  is  believed  that  the  above  cited 
characteristics  are  the  most  important  to  be  considered  in  selecting  an 
instrument  for  surveying  and  measuring  water  supplies  and  for  other  general 
uses.  The  manufacturer  of  each  insti^ument  may  be  determined  by  referring 
to  the  "Commercial  Source"  column  in  the  table  and  then  to  the  company 
number  listed  in  Table  XIX. 


Table  XIX 

Equipment  Manufacturers 


1.  Atomic  Accessories,  lac. 

8ll  W.  Merrick  Road 
Valley  Stream,  New  York 

2.  Baird  Atomic,  Inc. 

33  University  Road 
Cambridge  38,  Mass. 

3.  Eberline  Instrument  Corp. 

P.  0.  Box  279 

Santa  Fe,  New  Mexico 

4.  Victoreen  Instrument  Co. 

5806  Hough  Avenue 
Cleveland,  Ohio 

5.  Lionel/Aaton  Electronic  Laboratories 
1226  Flushing  Avenue 

Brooklyn  37,  New  York 

6.  Nucleonic  Corporation  of  America 
196  Degraw  Street 

Brooklyn  31,  New  York 

7.  Radiation  Counter  Laboratories,  Inc. 
5121  W.  Grove  Street 

Skokie,  Illinois 


8.  Radiation  Equipment  and 

Accessories  Corp. 

665  Merrick  Road 
lynbrook,  New  York 

9.  Technical  Associates 

l40  W.  Provldencia  Avenue 
Burbank,  California 

10.  Tracerlab 

1601  Trapelo  Road 
Waltham  5^,  Massachusetts 

11.  Nuclear-Chicago  Corporation 
333  East  Howard  Avenue 

Des  Plaines,  Illinois 

12 .  Franklin  Systems 
2734  Hillsboro  Road 
West  Palm  Beach,  Florida 

13.  Nuclear  Corporation  of  America 
Instrument  and  Control  Division 
2  Rlchwood  Place 

Denville,  New  Jersey 
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Using  the  information  in  Figure  15  fco  convert  the  instrument 
sensitivities  in  Table  XVIII  into  equivalent  concentration  values^  Table  XVII 
has  been  compiled  to  shov  the  concentrations  of  typical  isotopes  which  can 
just  be  detected  with  instruments  of  the  sensitivities  stated^ It  must 
be  borne  in  mind  that  the  dose  calibration  in  mr/hr  for  most  instruments 
is  done  with  a  cobalt  or  radium  groxand  source,  and  that  the  dose  extrapola¬ 
tion  to  other  isotopes  and  energy  ranges  will  not  be  the  same  for  different 
detectors . 

It  is  evident  from  the  foregoing  discussion  that  there  appear  to 
be  adequate  commercial  instruments  on  the  market  which  can  be  adapted  to 
give  a  reliable  indication  of  high-level  contamination  of  water  supplies. 

The  assay  of  isotopic  contamination  and  the  accurate  determination  of  safe 
concentrations  of  radioisotopes  in  water  by  field  instruments  under  emergency 
conditions,  however,  does  require  some  further  research  and  development 
work.  Such  equipment  will  have  to  be  made  available  in  sufficient  quantity 
for  extensive  wat«r  analyses  as  soon  after  an  attack  as  possible. 
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Figure  15. 


Relation  of  Dose  Rate  to  Nuclide  Concentration 
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VIII.  DECONTAMINATION  OP  WATER  SUPPLIES 

As  discussed  in  Quarterly  Technical  Report  #3  (see  Appendix  A) 
radiologically  contaminated  waters  have  been  subjected  to  decontamination 
procedures  of  both  conventional  and  non-convent ional  nature .  Of  the  six 
radioisotopes  of  biological  Importance,  all  have  been  found  to  be  amenable 
to  removal  to  some  degree.  The  degree  of  treatment  required  has  been 
predicted  in  previous  reports,  predicated  by  the  Maximum  Permissible 
Concentration  for  peacetime  consumption  of  water  or  exposure  over  a  con¬ 
siderable  period  of  time.  These  requirements  are  very  unrealistic  when 
considered  in  the  light  of  the  occurence  of  a  nuclear  incident.  The 
removal  efficiencies  for  the  most  part  have  been  derived  from  laboratory 
data  on  small  volumes  of  water  treated  under  rigidly  controlled  conditions. 
Of  practical  Importance,  only  those  data  obtained  by  actual  treatment  plant 
operation  are  to  be  considered  reliable  on  a  major  scale.  Even  in  this 
respect,  it  must  be  pointed  out  that  these  data  were  derived  from  analysis 
of  water  supplies  which  had  been  contaminated  by  long-range  weapon  debris. 

It  has  been  discussed  elsewhere  in  this  report  that  the  long-range  fallout 
is  more  soluble  than  short-range  or  early  fallout.  It  is  for  this  reason 
that  we  may  expect  full-scale  operations  to  be  more  efficient  in  radio¬ 
isotope  removal  of  early  fallout,  than  the  presently  reported  efficiencies. 

It  appears  that  water  decontamination  in  the  early  phases  of  a 
nuclear  incident  aftermath  are  not  beyond  jthe  realm  of  practical  application 
of  the  basic  sanitary  engineering  principles  of  water  treatment .  The  major 
problem  involved  in  water  decontamination  will  be  the  long-term  exposure  of 
the  surviving  population  to  concentrations  of  the  long-halflived  isotopes 
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of  biological  importance.  The  reduction  of  these  contaminants  to  deslreable 
levels  will  constitute  the  primaiy  long-term  problem  to  the  water  treatment 
field,  but  is  by  no  means  an  insurmountable  task. 
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IX.  BODY  BUBDEN  YRm  CONTAMINATED  WATER 

In  the  analysis  of  Body  Burden  from  continued  ingestion  of  water 
contaminated  with  radioactivity  resulting  from  nuclear  attack,  it  has  Been 
found  that  the  development  of  a  rational  formula  is  desirable.  A  formula 
has  Been  developed  in  this  report  which  gives  values  of  Body  Burden  com¬ 
paring  closely  with  those  obtained  By  empirical  means. 

Many  of  the  empirical  values  used  in  this  analysis  were  derived 
from  animal  studies,  and  may  or  may  not  Be  applicable  to  man  per  se; 
however,  the  values  obtained  By  formulas  developed  in  this  report  are  in 
good  general  agreement  with  data  collected  from  experiments  with  man. 

Data  from  man  has  Been  used  wherever  possible;  unfortunately  experiments 
with  man  are  not  sufficiently  comprehensive  to  serve  as  the  entire  Basis 
of  an  analysis. 

The  relation  of  Body  Burden  to  many  of  the  radioisotopes  found 
in  fallout  is  not  easily  understood  for  short  periods  of  time,  that  is, 
when  the  study  period  is  a  very  small  fraction  of  the  physical  half-life 
of  the  radioisotope.  In  this  report,  in  addition  to  study  of  the  Immediate 
effects,  several  cases  are  analyzed  for  periods  of  time  in  excess  of  fifty 
years.  (Figures  l8  and  20) 

Nijmerous  studies  have  Been  made  on  the  affinity  of  certain 
types  of  tissue  for  particular  radioisotopes.  These  selectivity  studies 
are  of  special  interest  to  anyone  who  wishes  to  set  up  theoretical  values 
which  may  not  Be  exceeded  without  permanent  damage,  that  is,  such  values 
as  "Maximtun  Permissible  Concentration,"  "emergenc''-  limits,"  etc.  Although 
it  is  not  within  the  scope  of  this  report  to  determine  such  values,  it  is 
of  interest  to  note  that  allowable  Body  Burden  is  the  total  concentration 
of  a  given  Isotope  allowed  in  the  most  critical  organ.  The  most  critical 
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organ  is  that  organ  which  receives  its  damaging  dose  first  during  the  chain 
of  hiological  processes. 

The  phenomenon  involved  in  the  rate  of  turnover  of  various  isotopes 
hy  a  particular  organ  is  of  interest  since  this  rate  may  he  equal  to, 

(75) 

greater  than,  or  less  that  the  rate  for  the  entire  hody.'’'-^'  A  sample  plot 
of  the  critical  organ  versus  physical  half-life  is  shown  in  Figure  I6. 

Both  hiological  and  physical  decay  rates  are  considered  in  this 
report.  In  order  to  demonstrate  the  difference  in  body  burden  levels  ob¬ 
tained  by  (1)  considering  only  physical  decay  and  (2)  considering  both 
physical  and  biological  decay,  i.e.  effective  decay,  A^  values  for  both 
cases  afe  plotted  in  Figures  16  through  21  as  curves  B  and  A  respectively. 

values  plotted  as  curve  A  are  parameters  of  the  rate  at  which 
the  radioactivity  level  in  the  body  is  reduced  by  a  combination  of  bodily 
processes  and  physical  decay.  These  A^  values  are  therefore  parameters  of 
the  effective  decay  rates  of  the  various  radioisotopes.  The  amount  of  any 
radioisotope  which  passes  through  the  body  must  further  be  related  to  the 
amount  retained  by  the  critical  organ.  This  relation  is  represented  by  an 
f  value  which  gives  the  amount  of  a  radioisotope  taken  up  by  the  organ 

"W 

critical  for  that  isotope. 

In  the  interest  of  generality  we  have  developed  a  method  of  compu¬ 
tation  (using  A^  values)  rather  than  a  list  of  specific  figures  for  a  limited 
number  of  cases. 

No  attempt  has  been  made  to  evaluate  cell  damage  resulting  from 
various  levels  of  body  burdens.  It  is  obvious  that  the  greater  the  body 
burden  the  more  extensive  is  the  radiation  damage.  Examples  are  cited  in 
this  report  of  some  of  the  biological  blocking  processes  which  may  be 
applied  to  effectively  reduce  the  body  burden.  These  A^  values  will  give 
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the  highest  burden  resulting  only  from  drinking  contaminated  water  while 
eating  a  normal  uncontaminated  diet  and  breathing  uncontaminated  air. 

With  an  easily  derived  scaling  factor,  however,  these  same  graphs  may  be 
applied  to  ingested  solids.  values  may  be  obtained  either  from  fomulas 
to  be  derived  later  in  this  section  or  from  the  graphs  in  Figures  l6 
through  21.  The  graphs  are  easier  to  use  and  are  included  as  an  aid  to 
computation . 

No  limitations  of  physical  death  or  early  genetic  death  of  the 
individual  under  study  have  been  placed  on  the  data  shown  here,  since 
both  these  limitations  require  an  analysis  of  damage.  We  cannot  overlook 
the  possibility  that  at  certain  levels  of  contamination  or  consumption 
physical  or  genetic  death  will  limit  the  continuity  of  our  data,  as  will 
any  departure  from  the  size  and  age  of  a  staindard  man. 

Establishment  of  a  Criteria 

(a)  The  data  to  be  presented  should  be  in  the  most  useful  and 
general  form  applicable  to  the  subject  matter. 

(b)  Analysis  of  damage  that  may  occur  at  the  various  levels  of 
body  burden  is  not  to  be  attempted;  however,  comment  regarding  damage  in 
general  will  be  incorporated  to  give  depth  to  the  conclusions. 

(c)  The  development  of  a  rational  method  for  evaluating  the  body 
burden  is  to  be  outlined  and  explained.  The  correlation  of  the  values 
obtained  by  the  use  of  this  method  and  those  obtained  by  empirical  means 
will  be  discussed. 

(d)  Assumptions  should  be  as  few  and  as  reasonable  as  possible 
in  order  to  maintain  the  Integrity  of  the  data  to  be  presented. 

(e)  All  values  used  for  the  various  calculations  shall  be 
those  derived  from  the  most  authoritative  Information  available. 
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As  sumptions 

(a)  The  hiological  decay  rate  used  is  that  of  a  standard  man. 
Since  the  agp  and  overall  metaholism  of  the  Individual  will  affect  the 
turnover  and  replacement  rate  of  radioactive  materials  in  the  hody^  and 
therefore  the  hiological  decay  rate^  it  is  necessary  to  standardize  these 
parameters  for  the  present  study. 

(h)  The  contamination  of  an  individual's  water  supply  takes 
place  all  at  one  time,  and  no  further  contamination  occurs.  The  analysis 
of  a  fluctuating  contamination  level  would  require  that  we  assume  rates 
of  fallout,  wind  direction  and  numerous  other  meteorological  conditions. 

In  order  to  retain  the  generality  of  our  investigation  we  have  taken 
a  case  in  which  weather  and  attack  data  were  not  applicable .  curves 
may  be  added  to  each  other  if  new  contamination  occurs  (possibly  from 
fallout  from  a  different  burst  or  fallout  return  from  wind  reversal). 

(c)  The  radioactive  contamination  of  the  water  supply  is  reduced 
only  by  physical  decay.  This  assumption  is  conservative  in  nat’jre,  since 
it  disregards  entirely  the  amount  of  activity  which  might  settle  out  or 

be  consumed  by  life  forms  found  in  the  water  supply,  e.g.  plankton,  algae, 
fish. 

(d)  Effective  decay  rates,  even  if  extrapolated  from  animal 
data,  are  applicable  to  man.  In  the  absence  of  experimental  data  on  man, 
the  animal  data  used  ’,ras  the  best  available.  It  may  be  of  interest  to  note 
that  many  contemporary  investigators  are  not  satisfied  with  many  of  the 
extrapolations  of  animal  body  burden  data  that  have  heretofore  been  accepted 
as  correct,  or  at  least  useable,  for  the  case  of  humans. 

(e)  The  dally  intake  of  water  is  constant  and  takes  place 
Instantai-ieously.  Virtually  all  of  an  individual's  water  intake  occurs 


in  sixteen  hours,  and  over  half  occurs  in  ei^t  hours;  however,  the  errors 
associated  with  this  simplification  tend  to  he  compensating  if  the  individual's 
drinking  hahits  are  consistent.  For  reasons  of  generaliny  we  are  forced  to 
consider  that  the  man  under  study  is  a  creature  of  hahit.  By  drawing  a 
smooth  curve  through  the  daily  points  we  have  essentially  integrated 
over  the  time  periods  between  the  daily  values  and  have  described  a  contin¬ 
uous  function  as  if  1die  individual  had  spread  his  intake  over  the  entire 
day.  The  combination  of  the  two  effects  results  in  a  close  approximation 
of  the  actual  build-up  of  body  burden. 

Development  of  a  Rational  Method  of  Analysis  of  Body  Burden  for  Continuous 
Consumption  of  Water  Contaminated  with  Radioactive  Elements 
Based  on  the  Geiger-Kutall  Law  of  radioactive  decay  we  may  say  that  the 
amount  of  activity  present  at  any  time  is  equal  to 

— \  t 

A  =  A  e  ^ 
o 

when  only  physical  decay  is  considered.  is  the  physical  decay  constant 
for  a  particular  radioactive  element,  and  A^  is  the  amount  of  activity 
initially  present. 

Furthermore,  if  only  biological  decay  is  considered  we  may  say 

that 


where  is  the  biological  rate  of  elimination  constant. 

If  we  define  "Effective  Half  Life"  as  the  time  required  for  a 
radioactive  element  fixed  in  the  tissues  of  an  animal  body  to  be  diminished 
by  fifty  per  cent  as  a  result  of  the  combined  action  of  radioactive  decay 
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and  "biological  elimination,,  then  we  may  say  that  the  effective  half-life 


is 


(76) 


Biological  half-life  x  Radioactive  half-life  _  ^ 

Biological  half-life  +  Radioactive  half-life  '  "^eff 

An  effective  decay  constant  (^g^)  defined 

V  -  Q‘693 

eff  T 

"^eff 

Let  us  assume  iritially  that  for  any  activity  taken  into  the 
body  the  entire  amount  will  be  fixed  in  tissue.  It  is  then  apparent  that 
immediately  before  initiation  of  consumption  of  contaminated  water^  the  body 
burden  is  equal  to  whatever  radioactivity  is  present  as  a  result  of  pre¬ 
vious  activity.  For  the  development  of  our  argument  we  shall  take  this 
initial  body  burden  to  be  zero  and  correct  it  in  later  developments  for 
the  case  of  initial  body  burden  not  equal  zero. 

The  following  reasoning  may  then  by  applied: 

(1)  Just  before  initiation  of  consumption  activity  in  the  body 
equals  zero. 

(2)  Immediately  following  the  start  of  consumption  the  activity 
in  the  body  equals  A^,  that  is,  the  activity  per  unit  volume  of  the  water 
(pc/ml,  atoms/ gal,  etc.)  times  the  volume  of  water  consumed.  Since  we 
have  assumed  that  the  amount  of  liquid  is  constant  we  may  drop  the  volume 
term  and  reincorporate  it  later  in  the  development.  Furthermore,  by 
considering  A^  equal  to  unity  we  will  develop  an  equation  which  when 
multiplied  by  the  actual  initial  activity  will  give  activity  present  at 
any  given  day  n.  We  wish  to  develop  a  plot  which  is  independent  of  both 
amount  consumed  and  initial  activity,  and  is  dependent  only  on  the  radio¬ 


isotope  in  question. 
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(3)  Immediately  prior  to  consumption  on  the  second  day  the  "body 
burden  would  be 


A,  =  A  e 
1  o 


-w 


and  just  ai’ter  consumption  the  body  burden  is 

-4  (1) 

A'  =  Ae  '+Ae  ^ 

1  o  ■  o 


(4)  Just  before  consumption  on  the  third  day  the  body  burden 


would  be 


A,  =  A|e  =  (A^e 


and  immediately  following  ingestion,  the  expression  for  body  burden  may 
be  written 

1  1  1  1  1 

Ai  *  (A  e  ^  )e  ^  +  (A  e  +  A  e  ^  )e 

d  '  o  '  '  o  o 


(5)  Continuing  this  reasoning,  a  recursion  formula  is  developed 
of  the  form 


A  =  A  [e  +  e 

n  o 


eff 


+  (a+l)>>.  1]  1 


-  -  -  +  e 


eff 


_  -nX  1 

^  +e  ^]-Ae  ^ 


However,  rather  than  deal  with  the  (n  +  l)  terms  necessary  to  evaluate 
this  equation  at  n  days  we  simplify  the  equation  as  follows 

n  E 


A 


A  b  ) 

o  '  b  -  p  ' 


where  b  =  A  e 
o 


and 


p  =  A  e 
o 
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(5)  'j?his  seneral  expression  for  is  the  value  plotted  as 
curve  A  in  Flijurcs  i6  bhrough  21. 

(7)  We  may  now  incorporate  into  our  formula  the  portions  of  total 
body  '‘lurden  which  have  thus  far  been  neglected. 

(a)  Since  we  have  considered  to  be  unity,  we  may  find 
the  total  body  burden  by  multiplying  A^  by  the  true  original  activity. 


A  ,  ,  and  the  volume  consumed  daily,  D. 

o  true  '  ’ 

atoms \  .  / atoms  > 


Th: 


jy  Ti  ^  -n  ■  /  auoms  \  .  /  aijoins  \  liters  \ 

rafore  Body  Buraanj  {^)  .  x  x 


A  /days  of  consumption \ 
n  body 


(b)  Tlie  assumption  that  ail  the  activity  taken  into  the 
body  is  immediately  incorporated  into  the  intestine  is  highly  improbable. 
We  may  however,  consider  the  amount  incorporated  into  tissue  as  a  fraction 
of  the  ingested  activity  reaching  an  or^r.  (f  )  and  multiply  the  body 

V 

burden  of  (a)  above  by  f  to  obtain  a  more  reasonable  value. 

w 

Thus  Body  Burden.rT  =  f  xA^  xDxA 

II  wo  true  n 

Some  values  of  f^^  for  isotopes  under  consideration  are  given  in  Table 

(c)  While  reducing  the  body  burden  by  the  factor  f  we  have 

w 

neglected  the  activity,  A„_,  which  is  present  in  the  intestinal  tract 

ui 

on  day  n.  This  activity  is  an  addition  to  the  body  burden,  therefore; 


Body  Burdera-TT  =  f  xA.  xDxA  + 

111  w  o  true  n  GI 


where  by  the  reasoning  and  definitions  previously  expressed 

■1.25k, 


"GI 


-k  (n-0.25) 

(1  -  f  )(A  ,  e  P  )(e 

w  o  true 
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Since  the  time  required  for  ingested  material  to  pass  through  the  tody 

is  31  hours  or  about  1.25  days,  the  term  is  strictly  transitory. 

It  may  or  may  not  he  negligible  depending  on  the  period  in  which  it  takes 

place.  On  the  first  day,  for  example,  it  would  be  the  largest  portion  of 

the  total,  but  by  the  time  the  fractional  half  life  (n/'Ii)  is  eqioal  to 

2 

one,  the  term  may  in  some  cases  be  considered  negligible,  (n  equals 

time  in  days  since  beginning  of  consumption  and  Tj^  eqtxals  radioactive 

2 

physical  half-life  in  days.) 

(d)  The  only  term  still  missing  from  the  mathematical  expression 
is  that  for  the  activity  present  in  the  body  at  the  time  drinking  of 
contaminated  water  started  (A^).  That  which  was  already  present  in  the 
tissue  at  time  of  Initial  Ingestion  will,  after  n  days  of  effective  decay, 
equal 


\  -  *1' 


That  present  in  the  gastro-lntestinal  tract  at  the  time  of  initial  ingestion 
will  give  an  additional  activity  equal  to 


where  A  is  the  pre-ingestion  activity  in  the  Intestinal  tract . 
u 

(8)  After  combining  the  various  components,  the  expression 
for  body  burden  from  any  particular  radioisotope  at  day  n  may  be  written: 

(n-1.25)  -1.25V 

P  D 


Body  Burden  = 


n 


o  true 


)(e 


)  + 


w'''o  true^^°^'\^  ^GI  \  ^ 
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The  full  explanation  of  this  development  has  been  presented  in 
order  to  allow  the  user  of  the  quations  to  decide  which,  if  any,  terms 
may  be  neglected  for  any  specific  analysis.  Should  the  individual  in 
question  be  an  average  U.  S.  citizen  w^ho  has  so'uglit  shelter  during  an  attack. 


the  term 


*1  *  'A  ^  \  \ 


will  for  practical  purposes  equal  zero.  As  previously  noted,  the  value 

of  will  generally  be  negligible  after  a  relatively  short  period  of 

consumption.  If  the  (A^  +  A^)  and  A^^  terms  are  dropped,  the  practical 

expression  for  body  burden  is  (f  )(a  .  )(D)(A  ).  The  A_^  term,  however, 

w  o  true  n  GI 

may  be  of  major  importance  when  the  physical  half-life  is  much  greater 
than  the  effective  half-life. 


It  must  also  be  stressed  that  to  evaluate  the  effect  of  Ar..T> 

Ay  and  A^  on  any  individual  that  individual's  history  of  exposure  to 
radiation  must  be  known.  It  is  only  by  thorough  knowledge  of  the  develop¬ 
ment  of  this  rational  method  that  one  may  make  a  reasonable  calculation  of 
body  burden. 

Reduction  of  Body  Burden  by  Blockin 


It  has  been  found  that  in  certain  cases  uptake  of  radioisotopes 

may  be  reduced  by  biological  blocking.  For  example,  by  administering  a 

soluble,  non-toxic,  inert  compound  of  iodine  one  can.  prevent  uptake  of 

radioactive  iodine  by  the  thyroid.  Ten  to  one  hundred  milligrams  of 

potassim  iodide  taken  at  one  time  will  prevent  the  uptake  of  any  more 

iodine,  radioactive  or  otherwise,  by  the  thyroid  for  about  one  week.  To 

continue  this  effect,  it  is  believed  that  an  additional  200  milligrams  eveiy 

( 77) 

two  days  would  be  required , '  Blocking  lowers  the  value  of  f  in  body 
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■burden  calculations. 

In  the  case  of  Strontium  89  and  Strontium  90,  which  "behave  lihe 
calcium,  it  is  "believed  that  "blocking  may  be  completely  or  partially 
accomplisiied  by  consuming  a  calcium  compound  far  in  excess  of  norml 
requirements.  This  is,  however,  a  conclusion  made  on  the  basis  of  extrapo¬ 
lated  dataj  the  only  case  in  which  it  has  definitely  been  proved  that  blocking 
is  effective  is  that  of  iodine. 

Should  the  body  be  initially  deficient  in  a  given  stable  isotope 
(possibly  from  deficiency  of  a  common  mineral  in  the  pre-ingestion  diet) 
the  value  of  f^  for  the  radioactive  isotope  of  that  mineral  or  a  chemically 
similar  mineral  may  be  increased  because  of  the  increased  affinity  of  the 
body  for  that  mineral. 

Blocking  should  not  be  considered  a  panacea  to  the  problem  of 
body  burden,  since  it  may  result  in  the  cure  being  worse  than  the  disease. 

In  the  iodine  example,  for  instance,  the  prescribe(’  '!  ’Ounts  of  potassium 
iodide  constitute  a  daily  inr^t  of  iodine  into  t>'  )od  »ne  thousaiid 
times  normal.  Although  it  is  believed  that  for  t  iority  of  the  U.  S. 
population  this  treatment  would  not  be  harmful  for  twL  three  weeks, 
the  secretion  rate  of  the  thyroid  is  decreased.  Extended  treatment  will 
produce  other  still  more  harmful  side  effects  which  have  not  yet  been 
fully  evaluated.  Blocking,  however,  may  play  an  important  role  in  com¬ 
puting  body  burden  values  and  should  not  be  overlooked  when  selecting 

f  values . 
w 

Discussion 

(1)  In  Figures  I6  through  21,  curve  A  gives  the  alue  resulting 
from  the  effective  decay  rate.  Curve  B,  included  for  purposes  of  comparison, 
represents  the  rate  at  which  physical  decay  alone  takes  place.  It  seems  to 


-1*1-8- 


point  out  that  because  of  the  body's  tendency  to  eliminate  certain  radio- 
isotopes j  or  not  to  ihcorporate  them  into  tissue  at  all.  body  burdens  are 
significantly  lowered. 

(2)  T-ible  XXIV  gives  a  few  selected  values  of  f^  for  six 
radioisotopes  and  various  organs.  151606  are  values  which  result  from  a 
normal  uncontaminated  diet.  Values  of  f^  may  vary  from  this  table  according 
to  the  excess  or  deficiency  of  certain  minerals  in  the  body.  Note  that  the 
concentration  of  specific  isotopes  must  be  known  before  body  burdens  can 
be  computed. 


Table  XX 

Selected  Values  of  f  for  Various  Organs 

w 


Isotope 

Organ 

f 

w 

Bone 

2.5  X  10"^ 

Bone 

2.5  X  10’^ 

Ru^°^ 

Kidney 

2.0  X  10"^ 

1I3I 

Thj’roid 

2.0  X  10"^ 

Cs137 

Whole  Body 

1.0  X  10° 

Cs137 

Liver 

*<-.8  X  10"^ 

Bone 

7.0  X  10"^ 

Notes : 

(1)  The  f^  values  shown  are  for  critical  organs.  Ihe  critical 
organ  will  determine  the  since  it  is  the  one  used  in  computing 

(2)  The  f  value  for  the  entire  body  is  the  sijmmation  of  f 
values  for  all  the  individual  organs  concerned . 


1  2 
FRACTIONAL  HALF-LIFE 

versus  Time  for  Cesium  137 


Figure  20, 
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FRACTIONAL  HALF-LIFE 


A  versus  Time  for  Barium  ikO 
n 


Figure  21. 
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A 

A 

^GI 


A' 

s 


a 


empirical  constant  for  sutilimation  or  vaporization  reaction 
radioactivity  after  some  period  of  decay 

radioactivity  already  present  in  the  intestinal  tract  at  time 
of  ingestion  of  contaminated  water 

radioactivity  present  in  the  intestinal  tract  at  some  time  after 
ingestion  of  contaminated  water 

radioactivity  already  incorporated  into  body  tissue  at  time  of 
ingestion  of  contaminated  water 

a  parameter  of  the  effective  decay  rate  in  the  "body  of  a  given 
isotope 


Initial  radioactivity  per  unit  volume  of  contaminated  water 
taken  into  the  body 


(a  f  Nk) 
f 


V 


10' 


12^ 


radioactivity  present  in  "body  tissue  resulting  from  decay  of 
that  activity  Initially  present  in  tissue  (A  )  after  some  time 
of  effective  decay 


activity  at  time  t 

activity,  in  disintegrations  per  second,  of  the  radionuclide 
in  a  target,  after  the  nuclide  has  been  removed  from  the  flux 
for  a  period  Oj 

radioactivity  incorporated  into  tissue  from  that  activity  originally 
present  in  the  intestinal  tract  (A_) 

cloud  horizontal  semi-axis 


a^  fireball  horizontal  semi-axis  at  ground  zero 

a„  fireball  horizontal  semi-axis  at  Z 

Li 

B  empirical  constant  for  sublimation  or  vaporization  reaction 

B  ratio  of  fission  to  total  yield 

Body  Burden  radioactivity  incorporated  into  body  tissues  as  a  result  of 
ingestion  of  radioactivity 


b 


cloud  vertical  semi-axis 
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C 

D 

D 

°X 

d(t,  1) 

FD^(t) 

FP  (t) 
r 

f 

g 

h 

W*) 


a 

"j 


GLOSSARY 

(Continued) 

fireball  vertical  semi-axis  at  ground  zero 
fireball  vertical  semi -axis  at  Z 

empirical  constant  for  sublimation  or  vaporization  reaction 
instrument  response  factor 
daily  water  Intake,  volume 
exposure  dose 

decay  correct  ion  factor  for  H  +  1  hour 
fraction  of  device  contour  ratio 
fission-product  contour  ratio 

thermal  flux  in  neutrons  per  square  centimeter  per  second; 

fugacity  of  the  element  in  the  liquid  phase 

a  coefficient  relating  ingested  radioactivity  to  Body  Burden 

acceleration  due  to  gravity 

cloud  center  height 

fallout  intensity  at  time  t 

air  ionization  rate  per  fission  at  3  ft  above  an  infinite, 
ideal  plane  for  a  uniform  distribution  of  the  normal  fission 
product  mixture 

air  ionization  rate  per  fission  at  3  ft  above  an  infinite, 
ideal  plane  for  a  vuiiform  distribution  of  the  condensed 
fission  product  mixture 

air  ionization  rate  per  fission  at  3  ft  above  an  Infinite, 
ideal  plane  for  a  uniform  distribution  of  the  neutron  induced 
activities 

relative  abundance  of  the  isotope  from  which  the  radionuclide 
is  formed 

arbitrary  constant 

fireball  horizontal  semi-axis  expansion  constant 
fireball  vertical  semi-axis  expansion  constant 
Henry's  Law  constant 
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(Continued) 


M^(t) 

m 

m 

N 


N 


fp 


N, 


H*(A) 


n 

n(i) 


n. 

J 

n'. 


J 


P 


k^/[n(i)A]RP 

inverse  time  constant 
mass  contour  ratio 

mass  of  fallout  per  unit  area  at  any  downwind  distance,  X 
mass 

disintegration  multiplier 

the  total  number  of  atoms  of  the  element  in  the  target; 
niomber  of  atoms,  or  moles,  of  fission  products  per  unit  area 
mole  fraction,  of  element  J  in  the  liquid  phase 
mole  fraction  of  element  j  in  the  vapor  phase 

number  of  atoms  (of  the  end  member  of  mass  chain  and  condensed  on 
the  outside  of  the  particle)  which  land  per  square  foot  of  ground 

number  of  moles  of  vapor 

total  moles  of  liquid  carrier 

abbreviation  of  n..(t) 

JA 

amount  of  element  j  condensed  on  the  surface  of  the  solid 
particles 

amount  of  element  j  in  the  vapor  phase 

abbreviation  of  n°  (t) 
jA 

number  of  moles  of  element  J  with  mass  number  A  dissolved 
in  the  n(i)  moles  of  liquid  carrier,  which  is  the  particle 
in  the  liquid  phase,  prior  to  solidification 

number  of  moles  of  element  j  with  mass  number  A,  which  is 
mixed  with  moles  from  other  mass  chains  to  form  n  moles  of 
vapor  and  not  condensed  in  the  liquid  carrier 

empirical  constant  in  particle  falling  velocities  relation 

sublimation  or  vaporization  pressure 

partial  pressure  of  the  liquid 
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s 


q 


R 

r^CA,  t) 
r;(A,  t) 
r^Ct) 

T 

T 

T 

eff 


T. 

t 

t 


V 

V, 


w 


W 

X 

X 


suTDlimation  pressure 

empirical  constant  in  particle  falling  velocities  relation 

terrain  shielding  factor 

molar  gas  constant 

initial  firehall  spherical  radius 

fractionation  number  of  the  first  period  of  condensation 

fractionation  number  of  the  second  period  of  condensation 

gross  fission  product  fractionation  nxmiber 

absolute  temperature 

half-life  of  the  radionuclide  formed 

time  reqxiired  for  a  radioactive  element  fixed  in  the  tissue 
of  an  animal  body  to  be  diminished  50  per  cent  as  a  result 
of  combined  action  of  physical  decay  and  biological  elimination 

radioactive  (physical)  half-life 

time  of  irradiation 

time 

time  of  arrival  for  particle  falling  out  of  the  cloud 
time  of  cessation  for  particle  falling  out  of  the  cloud 
particle  falling  time  in  the  fireball 
particle  rising  time  in  the  fireball 
molar  volume 

average  particle  falling  velocity 
wind  velocity 

instantaneous  particle  falling  velocity 
weapon  yield  in  klloton  TNT  units 
longitude  axis  in  ground  coordinates  system 
longitude  axis  in  cloud  coordinates  system 
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i^(t) 


latitude  axis  in  ground  coordinate  system 

total  amount  of  radionuclides  for  the  entire  chain  yield  of 
mass  number  A 

total  amount  of  radionuclides  of  element  j  present  at  time  t 

total  number  of  atoms  of  element  j  of  all  mass  numbers 

maximum  half -width  of  the  crosswind  distance  on  a  1  r/hr 
contour 

latitude  axis  in  cloud  coordinate  system 

amount  of  radionuclides  of  element  j  and  mass  number  A 
present  at  time  t  after  fission 

number  of  atoms  of  element  j  of  mass  number  A  not  condensed 
in  the  liquid  carrier 

altitude  axis  of  ground  coordinate  system 

yield  dependent  mxlLtiplier  in  fireball  altitude  equation 

altitude  axis  in  cloud  coordinate  system 

smallest  Intercept  of  particle  falling  slope  with  cloud 

largest  intercept  of  particle  falling  slope  with  cloud 

particle  size  parameter 

maximum  or  minimum 


a  constant  expressing  the  biological  rate  of  elimination  of 
radioactivity 

a  constant  expressing  the  rate  of  attenuation  of  radioactivity 
in  the  body  as  a  result  of  both  physical  decay  and  biological 
elimination 


disintegration  constant}  the  fraction  of  the  number  of  atoms 
of  a  radioactive  nuclide  which  decay  in  unit  time 

activation  cross  section  in  square  centimeters  for  2200  m/ sec 
neutrons 


e 


time  of  decay 
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AFPEHDIX  A 


Excerpts  from  Quarterly  Technical  Report  #3 


F.  Preliminary  gyaluation  of  Water  Supply  Contamination  for  Selected  Cities 
1 .  Int  roduct 1 on 

In  order  to  ha-ve  some  basis  for  evaluating  the  relative  hazard  of 
nuclear  fallout  to  the  vater  supplies  of  different  cities,  attacks  have  been 
assumed  on  each  city  theit  would  result  in  the  same  intensity  over  each 
watershed  involved. 

From  this,  the  contamination  of  the  water  due  to  the  suluble  portion 
of  the  fallout  has  been  calc\iLated  for  H  +  1  hour.  Simplifying  assumptions 
were  made  in  regard  to  an  average  intensity  over  the  specific  watersheds  and 
8.S  to  the  time  factors  involved.  Ibese  assumptions  were  made  to  provide  a 
basis  for  the  calculations  of  the  activity  concentrations  shown  in  Table  IV. 
These  values  are  shown  principally  to  demonstrate  the  approach  being  taken 
and  to  give  an  idea  of  the  maximum  contamination  that  could  be  expected 
immediately  following  a  nuclear  emergency. 

The  approach  to  be  taken  in  future  study  will  be.  more  specific. 
Factors  that  will  be  taken  Into  account  to  give  somewhat  more  reasonable 
values,  will  ber.  (l)  a  realistic  superimposition  of  intensity  contours  upon 
the  specific  watersheds  with  regard  to  seasonal  wind  directions  and  target 
sites,  (2)  consideration  of  decay  factors  to  give  expected  activity  at  times 
later  than  1  +  1  hour,  (3)  consideration  of  watershed  characteristics  such 
as  rorioff  coefficient,  feed-stream  velocities,  reservoir  draft,  transport 
phenomena,  and  reservoir  mixing  properties,  (4)  integration  of  intensity 
contours  to  give  a  more  precise  value  for  specific  isotopic  activities  on 
the  watershed  and  (5)  a  conversion  of  Intensity  to  a-ctivity  for  sites  off  the 
downwind  axis  employing  a  known  relationship  to  the  predetermined  values 


directly  downwind. 
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Table  XV 


:iT..y 


Activity  Concentrations  la  Water  of  8elected  leotopes 
Following  5  MX'  ^tuclear  Atts.eks 

Radloncciide  Cioncentration  in  Water  (lac/ml) 

Sr- 90  Ri-lOo  1-131 _  C^-X3'^ 


Sr-89 


i-lkO 


erA-er 

* 


2 .,7  x  lo'-^  2.4.  X  10  3-0  X  10  “■  3-0  x:  10  1-9  x  10  ""  2.2  x  10 


-4 


c;  ' 


i.o'"’  4-6  X  10"-^  5.6  X  10 


i.7 


3.5  X  10'-^  4.1 


1-r  A.;"*:  b 
Minn » 


7.3  X  10"^  7.2  X  10"^  6.2  X  10'^  8.0  X  lO"'^  4,9  x  lo"^  8.7  :x  lo"”^ 
lot  appl,i.cabl,6. 


Texa.-- 


4.5  X  lO*'^  4.1  X  10""  5v0  X  10'“^  p. 


0  X  10 


4.5  X  lO'-^'  4.1  X  10  '■  5.0  X  10"^  5.0 


3.1  X  iO  ^  3.6  X  10 
3,.l  X  10’-'  3.6 


-2 


"•ew  lork 

■I’J  •  i  . 


10-1  X  10'^  9-3  X  10'^  .L.l  X  iO'^  1.2  X  10"^^  7.2  x  lO"^  6.3  x  lO"'^ 
3.2  X  10"^  2.9  X.  lo"^  3.5  X  10' 3.5  X  lO"-"  2.2  X  lO"''  2,6  x  lO'" 


^pringfleld  1 


1.2  X  10*'^  1.5  X  10  ^  1.4  X  10"^ 


X 


10“"  8-5  X  lO"  9.9  X  10 


9.3  X  10  ^  8.5  X  !.-0  '  1, .  1  X  10  ^  1. .  i 


6.5  X.  10 


-4 


,6  X  iO" 


Ij;  'a-'e  1  applies  to  water  furiace  only 

OaEG  11  applle?  to  total  aiea  of  waterehed 
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By  employing  the  aforementioned  characteristics  and  variables,  it  is 
believed  that  activity  concentrations  will  be  reached  that  would  approach  a 
realistic  case. 

2.  Calculation  of  Activity  Concentration 

The  calculation  of  concentration  of  activity  used  in  the  preliminary 
study  is  based  upon  a  consideration  of  the  area  of  the  water  surface,  the  volume 
of  water,  and  the  total  watershed  area.  In  the  case  of  Duluth,  Minnesota,  where 
the  entire  municipal  water  supply  is  obtained  from  Lake  Superior,  only  the  depth 
of  the  lake  at  the  intake  point  was  considered. 

Current  information  on  the  characteristics  of  the  watersheds  was  obtained 
in  each  case  from  officials  of  the  municipal  water  supply  system  involved  or 
their  engineering  consultants.  Data  used  in  this  report  is  tabulated  in  Table  VI. 

TABLE  VI 

Characteristics  of  Watersheds  Serving  Selected  Cities 


Total  Watershed 

Total  Water 

Total  Water 

City 

Area  (ft.^) 

2 

Surface  (ft.  ) 

Volume  (liters) 

Denver,  Colo. 

8.29  X  10^° 

4.36  X  10® 

2.34  X  10^^ 

Houston,  Tex. 

7.8k  X  10^° 

7.85  X  10® 

2.50  X  10^^ 

New  York,  N.  Y. 

4.24  X  10^° 

1.37  X  10^ 

12 

1.90  X  10 

Springfield,  Mass. 

1.30  X  10^ 

1.68  X  10^ 

1.97  X  10^° 

As  an  example  of  the  procedme  used,  the  derivations  of  the  concentra¬ 
tions  In  the  Denver,  Colorado  water  supply  prior  to  any  municipal  treatment  are 
presented  here  in  detail.  The  other  cities  listed  in  Table  IV  were  analyzed  in 
the  same  manner  as  Denver,  except  for  Duluth  which  is  presented  separately. 
Denver  Case  I 

In  case  I,  only  that  radioactive  fallout  deposited  directly  on  water 
surfaces  was  considered  as  the  so'orce  of  the  water  contamination.  The  following 


assumptions  were  made  prior  to  the  calculation. 


-176- 


1.  5  M2?  10211)  attack. 

2.  Wind  speed  of  15  mph. 

3.  An  intensity  of  5^500  r/hr  over  the  entire  watershed. 

4.  Complete  homogeneous  solution  of  the  soluble  portion  of  the  fallout  in  the 
water  within  one  hour. 

5.  A  ratio  of  soluble  radioactive  atoms  per  square  foot  to  intensity  was 
chosen  as  that  value  for  each  isotope  at  51*5  miles  directly  downwind. 

6.  That  there  will  he  no  change  in  concentration  of  radioactivity  either  due 
to  reservoir  draft  or  feed-in. 

7.  Tt.at  only  Sr-89,  Sr-90,  Ru-106,  I-I3I,  Cs-137^  aud  Ba-l40  would  he  considered. 

0.  That  the  volume  and  area  of  feed  streams  will  not  affect  the  final  concentra¬ 
tion  appreciably. 


The  soluble  activity  density,  for  each  of  the  six  isotopes  being  coh- 
2 

sidered,  in  atoms/ft”^  was  calculated  as  follows: 

Given;  l(l)  »  5,500  r/hr, 

\(89) 


-TUT 

\(90) 

-lUT 

\(106) 

~iUT~ 

lUJ 

\(137) 

~iurr 

^(140) 

mr" 


6.17  X  10^*^  atoms/ft^^r/hr, 

11.1  X  10^^  atoms/ft^^/hr, 

4.87  X  10^^  atoms/ft^r/hr, 

10.7  X  10^^  atoms/ft^/:/hr, 

9*19  X  10^*^  atoms/ft^^/hr,  and 
12.4  X  10^*^  atoms/ft'^r/hr, 


where  the  subscript  4  is  the  value  of  the  particle  size  parameter,  a,  corresponding 
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to  the  downward  distance. 

By  multiplying  5,500  r/hr  by  IiJJ(a)/i(i)  for  each  isotope,  the  total 
number  of  soluble  atoms  of  each  isotope  per  square  foot,  n]|^(A),  will  be  obtained. 
This  value  for  each  isotope  was  found  to  be  as  follows: 


<(A)  .  1(1) 

'TOTT 

N*(89) 

=  3*^0  X  10“^  atoms/ft^ 

\(90) 

“  6.10  X  10^^  atoms/ft^ 

\il06) 

a  2.68  X  10^^  atoms/ft^ 

\il3l) 

=  5*89  X  10^^  atoms/ft^ 

nJ^(137) 

=  5*65  X  10^^  atoms/ft^ 

Nj(li^O) 

a  6.82  X  10^^  atoms/ft^ 

Multiplying  the  Nj^(a)  value  for  each  of  the  isotopes  by  the  total  area 
of  water  surface  will  give  the  total  number  of  atoms,  N(A),  of  each  isotope  in 
soluble  form,  deposited  in  the  water. 

N2^(A)  X  Area  =  N(A)  atoms 


N(89) 

= 

3. AO  X  10^^  X  A. 36  X  10^ 

=  lA  .8  X  10^^ 

N(90) 

= 

6.10  X  10^^  X  A. 36  X  10® 

=  26.6  X  10^^ 

n(io6) 

= 

2.68  X  10^^  X  A. 36  X  10® 

=  11.7  X  10^^ 

N(131) 

5.89  X  10^^  X  A. 36  X  10® 

PP 

=  25.7  X  10 

N(l37) 

s 

5.05  X  10^^  X  A. 36  X  10® 

PP 

=  22.0  X  10 

N(lko) 

s 

6.82  X  10^^  X  A. 36  X  10® 

22 

=  29.7  X  10 

atoms 

atoms 

atoms 

atoms 

atoms 

atoms 


Thus  knowing  the  total  number  of  soluble  atoms  of  each  isotope;  division  by 


-i  (5- 


the  total  water  volume,  V,  in  liters  will  give  the  oonoentration  of  each  isotope 
in  atoms  per  liter. 

* 


Ijl(A)  atoms 


liter 


f:.i) 


X  10^^  «  6.33  X  10^’  atoms/liter 

2.3k  X  10^ 

26.6  X  10^^  a  11.4  X  10^^  atoms/liter 
=  4.99  X  10^^  atoms/liter 


'89 


C* 

^90 


'106 


'131 


'137 


^140 


2.34  X  10' 


11 


,22 


11.7  X  10 

2.34  X  10^ 


.22 


^^•■7..^.  =  XI. 0  X  10^  atoms/liter 

2.34  X  10 


22.0  X  10 


22 


9.40  X  10^  atoms/liter 


2.34  X  10 


,11 


PP  1  1 

29.7  X  10*"*^  a  12.7  X  10-^  atoms/liter 


2.34  X  lo¬ 


ll 


(22) 


We  may  now  convert  the  activity  in  atoms/liter  to  microcuries/milliliter 
(|ic/ml)  -by  using  an  appropriate  conversion  factor  for  each  isotope. 

First,  to  change  from  atoms/liter  to  disintegrations  per  minute  per  liter, 

dpm  ,  the  foll.owing  relationshln  is  used, 
i 

C  ^  C*\ 
n  A 

Where:  •  Concentration  of  activity  in  dpm/i. 

=  Concentration  of  activity  Jji  atoms/i,  and 
^  =  Characterj,stic  time  constant  defined  by  the  equation 

-  4^ 

V 


Where :  T. 


1/2 


physical  half-life  of  the  specific  isotope  in  minutes, 
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p 

Since  there  are  2.22  x  lO"^  dpra/curie,  a  constant  for  all  isotopes, 
division  of  hy  this  value  will  give  the  concentration  In  curles/liter. 


n 


2.22  X  10' 


12 


2.22  X  10' 


12 


where  C  is  the  concentration  of  radioactivity  in  curies  per  liter, 
c 

Multiplication  of  hy  10^  will  convert  ciiries/liter  to  microcuries/ 
milliliter.  The  entire  conversion  factor  for  each  Isotope  is  then 


K 


•A 


2.22  X  10" 


The  conversion  factor,  K^,  for  each  of  the  six  isotopes  heing  studied 


are: 


K, 


■89 


K 

*^0 


=  1.25  X  10 


■15 


=  2.12  X  10 


-IT 


^106 


=  5-95  X  10 


-16 


K 


131 


=  2.72  X  10 


•ll 


K 


137 


K. 


140 


=  1.93  X  10 


=  1.70  X  1C 


-17 


•ll 


Bv  multiplying  each  conversion  factor  by  the  respective  concentration 
of  each  isotope  in  atoms  per  liter,  the  following  values  were  obtained. 

0(89)  =  2.7  X  10'^  ^lc/Inl 

C(90)  =  2.1  X  10"^  |ic/ml 

C(106)  =  3.0  X  lO"^  pc/mi 

C'(l3l)  =  3-0  X  10  ^  pc/ml 
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C(137)  =  1.9  X  10“^  nc/ml 
0(l40)  »  2.2  X  10"^  nc/ml 


Case  II 

In  this  case  the  soluhle  fraction  of  the  fallout  particles  deposited  over 
the  entire  watershed  was  assumed  to  have  been  transported  to  and  homogeneously 
dispersed  in  the  streams  and  reservoirs.  Other  assumptions  were  the  same  as 
for  Case  I.  To  continue  the  Denver  example,  in  Case  II,  the  watershed  area  of 
8.29  X  10^^  ft^  was  multiplied  by  the  Kj^  value  for  each  isotope  to  give  the 
total  number  of  atoms  of  each  Isotope.  The  remaining  calculations  were  carried 
out  in  the  same  fashion  as  for  Case  I. 

Duluth,  Minn, 

The  entire  water  supply  for  the  city  of  Diiluth  is  obtained  from  Lake  Superior, 
and  hence,  this  situation  is  different  from  that  at  the  other  cities  studied.  The 
Duluth  Intake  line  extends  out  II58  feet  into  the  lake  and  terminates  at  a  lake 
depth  of  72  feet.  The  fallout  deposited  over  this  Intake  was  considered  to  be 
evenly  dispersed  throughout  this  depth  of  water,  and  from  this  assuB^ition  the 
number  of  atoms  per  unit  volume  and  subsequently  the  activity  in  ^c/ml  was  calcu¬ 
lated. 

j.  Conclusions 

In  all  of  the  cases  thvis  far  studied,  the  time  factor  has  not  been  considered. 
EVen  at  an  assumed  wind  velocity  of  I5  miles  per  hour,  a  number  of  hours  would 
have  to  elapse  before  the  arrival  of  particulate  matter  at  the  far  reaches  of  a 
distant  watershed.  Dven  if  the  material  landed  directly  in  the  water,  much  more 
time  would  be  involved  in  transportation  in  streams  and  through  reservoirs  to  the 
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Intake  of  the  water  works.  There  would  he  further  delay  in  getting  into  streams 
those  particles  deposited  initially  on  dry  surfaces.  Another  delay  would  he  in 
the  occurence  of  rainfall  of  sufficient  intensity  to  dissolve  and/or  wash  this 
material  into  streams.  All  of  this  time  will  allow  radioactive  decay  to  proceed 
thus  ameliorating  the  hazardous  effects  of  the  fallout. 
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I .  Decontamination  of  Water  Supplies* 

1.  Removal  ty  Conventional  Water  -Treatment  Processes 

Conventional  municipal  water  purification  processes  generally  include 

aeration,  chemical  coagulation  with  sedimentation,  rapid  sand  filtration  and 
(79) 

chlorination'  .  To  a  lesser  extent,  lime-soda  ash  softening  and  slow  sand 
filtration  are  used.  Except  for  aeration,  the  processes  are  capable  of  removing" 
radioactive  contamination  to  some  degree,  either  singularly  or  in  combination. 

The  decontamination  capability  of  each  type  of  process  is  discussed  separately 
below,  with  particular  reference  to  the  six  elements,  barium,  cesiuDi,  iodine, 
lanthanum,  ruthenium  and  strontium, 
a.  Chemical  Coagulation 

The  most  common  coagulants  In  water  treatment,  aluminum,  and  iron  salts, 
form  alviminum  or  ferric  hydroxides  which  precipitate  as  floe.  The  chemical 
floe  acts  as  an  efficient  scavenger  by  adsorbing,  entrapping  or  otherwise 
bringing  together  suspended  ma,tter,  particularly  that  which  is  colloidal  in 
nature.  The  artificial  increase  of  the  alkalinity  in  water  may  also  form  the 
hydroxides  of  hea^’y  metals,  which  co-precipitate  with  the  aluminum  or  ferric 
hydroxide , 

Coagulation,  followed  b^,  sedimentation,  has  been  extensively  studied  by 
ORiTL^®*^^  with  reported  removals  of  only  36^  and  51^  for  Cs  and  Sr,  respectively. 
It  appears  that  the  method  is  effective  for  the  removal  of  suspended  or  colloidal 
material  and  for  most  cations  of  valence  3.  ^  or  5,  including  the  rare-earths 

group.  Removals  In  excess  of  98?^  have  been  reported  for  .P-32  as  the 

,  ,  ,  (8l)f82) 

orthophosphate' 

Removal  of  fa,  T  and  La  by  coagulation  was  investigated  by  Lacy'" 
with  a  maximum  of  found  for  both  Ba  and  La,  b"at  only  4$^  for  I. 


*  References  begin  on  page  155. 
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(84) 

Matsumura'  '  states  that  Ru  can  he  removed  up  to  925t>  hy  chemical 
coagulation. 

h.  Rapid  Sand  Filtration 

The  amount  of  radioactivity  removed  hy  filtration  -will  vary  depending  on 
the  nature  of  the  material.  Y  and  Zx,  prohahly  present  in  the  colloidal  state, 
■were  removed  from  up  to  93?^  hy  sand  filtration  alone  as  reported  in  a  com¬ 
prehensive  ORNL  report^®*^^,  whereas  materials  such  as  Sr  (4^)  and  Cs  (50?^) 
present  in  true  solution,  were  not  greatly  reduced  hy  passage  throng  sand 
filters.  They  further  reported  that  Ba  and  La  co'uld  he  removed  up  to  9?'t 
74^ ,  respectively . 

No  data  were  reported  on  the  removal  of  1  or  Ru  hy  rapid  sand 
filtration. 

c .  Chlorination 

Hannah,  et  studied  various  methods  for  the  removal  of  I-13I 

from  water,  and  found  that  small  dosages  of  chlorine  in  the  presence  of  100  ppm 
activated  carhon  produced  up  to  80^  removal.  The  authors^^^^  concluded  that 
the  only  effective  method  found  for  removing  I-I3I  with  materials  noimally 
available  in  water  treatment  plants  involves  chlorination  followed  hy  adsorption 
of  liberated  iodine  on  activated  carhon.  The  optimum  chlorine  dosages  were 
q-oite  small,  i.e.,  O.05  to  0.1  ppm.  The  removal  of  I-I3I  decreases  to  less 
than  2.Qfjo  when  the  dosage  of  chlorine  increases  to  Ippm.  Therefore,  normal 
pre -chlorination  employed  hy  the  water  treatment  plant  could  not  he  used  for 
iodine  removal  because  the  chlorine  residuals  would  generally  exceed  the  re¬ 
quired  dosage  for  activity  removal.  Stable  iodine,  in  dosages  greater  than 
0.01  ppm,  inhibit, ed  removal  of  I-I31  with  chlorine  and  Aqua  Nuchar  A 
(activated  charcoal).  Variation  of  pH,  achieved  hy  adding  s'alfuric  acid 
sodium  hydroxide,  was  found  to  have  i-'-'-tle  effect  on  the  removal  of  iodine 
hy  chlorine  and  activated  charcoal. 
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d.  Lime-Soda  Ash  Softening 

SofteEdng  plants  remove  excessive  amounts  of  scale-forming,  soap- 

consuming  compounds,  chiefly  comprised  of  the  cations  of  calcium  and  magnesium, 

hy  the  addition  of  lime  and  soda  ash  which  precipitates  calcium  as  a  carbonate 

and  magnesium  as  a  hydrate.  Because  of  the  similarity  in  the  chemical  properties 

between  strontium  and  calcium,  the  former  ion  is  co-precipitated.  For  most 

(or  ) 

satisfactory  removals  of  strontium,  Hoyt  '  fouM  that  excess  dosages  of  both 
lime  and  .soda  ash  were  required.  Removal  efficiencies  of  99 *7?^  ^  are 

possible  under  favorable  conditions.  Other  isotopes  which  can  be  effectively- 
removed  by  the  lime-soda  ash  are:  Ba,  Cd,  Y,  Sc,  Zr  and  Nb.  The  removal  of  Ba 
(99^t)  and  La  (905^)  was  reported  by  McCauley,  et  al  ORHL  ^  reported 

that  Cs  could  he  removed  by  lime-soda  softening,  but  not  efficiently  {23i>) . 

e.  Slow  Sand  Filtration 

Dovming,  et  al  '  ^ '  reported  that  the  slow  sand  filter  was  very 
effective  In  the  removal  of  radiostronti'um  in  the  first  feir  days  of  operation, 
but  that  the  efficiency  decreased  rapidly  and  essentially  drops  to  zero  at  ll<- 
days.  In  another  paper,  presented  by  Eden,  et  al  on  experiments  carried 

o^It  in  the  same  Water  Pollution  Research  Laboratory  in  Britain,  the  slow  sand 
filter  appeared  to  be  effective  for  only  a  few  hours.  The  activity  of  the 
effluent  rose  steadily  and  reached  30^6  of  the  initial  (imfiltered)  water  in 
one  day,  in  2  days,  and  95^^  after  7  days.  Eden,  et  aJ.  ^  also  reported 
that  iodine  was  not  efficiently  (50^)  removed  by  the  slow  sand  filter. 

2 .  Removal  by  Non-Convent ional  Treatment  Methods 

Below  are  listed  several  methods  of  decontamination  of  water  not 
commonly  enq?loyed  in  municipal  treatment  plants.  Economic  considerations  may 
limit  the  use  of  many  of  these  methods,  except  in  emergency  cases. 

a.  Ion  Exchftnpp 

A  few  water  treatment  plants  use  ion  exchange  as  a  softening  measure. 
Resins  have  been  found  to  provide  one  of  the  most  effective  methods  of 


decontamination.  Over  99 *9^  removal  of  Ba,  Cd,  Ce,  Cs,  I,  La  and  over  SSffo 

removal  of  P  and  Zr  was  q.uoted  by  OEKL  \  an-.^.  Amphlett  and  Sammon  ^ 

L  5 

reported  work  using  mixed  beds.  A  decontamination  factor  of  10  to  10  inay  be 
achieved  if  leakage  is  eliminated  ^ .  Che  natural  greensands ,  used  in 
some  municipal  softening  installations,  are  not  effective  in  the  removal  of 
anions . 

Although  the  ion  exchange  process  offers  one  of  the  most  promising 
methods  for  removing  radio-contaminants,  the  cost  may  preclude  widespread 
application  unless  suitable  regeneration  teehniq.ues  are  developed. 

The  possibility  of  using  home  water  softeners,  of  the  ion  exchange  type, 

(9^) 

has  been  discussed,  in  the  literature 


b.  Phosphate  Coagulation 

Because  the  conventional  method  employing  alum  or  iron  salt  as  the 

(9^  ) 

coagulant  does  not  effectively  remove  Sr  and  many  other  isotopes,  Lauderdale 
investigated  phosphate  coagulation  as  a  means  of  increasing  removals,  A  more 

) 

comprehensive  study  vb.s  ma,de  inter  by  Hesbltt,  et  al  .  Based  on  the 

theoretical  considerations,  as  well  as  experimental  observations,  Nesbitt  "  * 

concluded  that  the  mechaulsms  of  the  renoval  of  radioactivity  are  coprecipitation 

and  adsorption.  Nesbitt  reported  a  99^  strontium  removal.  The  three  variables 

which  were  foiuid  to  exeit  the  greatest  influence  on  coagulation  were  pH, 

calcium  concentmtion,  and  phosphate  concentration.  Careful  control  of  these 

vardables  was  necessary  to  achieve  optimum  strontium  removals. 

(^6 ) 

Tamuira  asLd  Sbruxness'''  '  reported  that  lew  cost  natural  phosphates 
could  be  used  in  place  of  sodium  phosphate.  The  authors  reported  that  over 


9C^  removals  were  achl.eved  with  these  minerals. 


Removal  of  87.8^  for  Ba,  98^  for  La  and  85^  for  Ru  have  been  reported 


by  Hatsumura 


'Jcr'’«'Eer  et  al 


report  only  351^'  removal  of  Cs  using 


phosphate  coagifLation. 
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Modified  Coagulation  and  Llme-Soda  Ash  Softening 

At  least  two  modified  coagulation  processes  have  teen  studied,  one  of 

which  is  the  continuous  addition  of  clay.  Cne  modification  of  the  process  is 

the  addition  of  small  amounts  of  silver  niti-ate  to  increase  the  removal  of 

iodide.  The  removal  of  95-3%  1-131  has  been  reported  hy  Eden,  et  al^®  \ 

(§9 ) 

Lacy  and  de  Laguna  ^  ^  report  removals  jf  99. 7^^  i'or  Ba  and  Os,  and 

;97 ) 

97. for  Ku  when  coag-jlation  is  supple iic....ced  b>  clay  addition.  Co-wser,  et  al’  '  ^ 
using  a  lime-soda  ash  softening  process  -with  the  addition  of  clay,  reported 
ruthenium  removals  of  l&f).  The  removal  of  cesium  averaged  Q6%  with  clay  doses 
of  200  ppm. 


d.  Metallic  Dusts 
flOO) 

Lacy ■  ^  conducted  laboratory  jar-test  studies  for  the  removal  of 

radioactive  contaminants  by  employing  various  concentrations  of  powdered 

aluminum,  copper,  iron  and  zinc.  In  general,  highest  removals  were  obtained 

(80) 


with  iron  powder.  Removals  of  37«2j&  of  I  were  reported 

(101) 


.  In  1962,  a 

patent  was  granted  to  Silier  '■‘‘"‘•'for  removing;?,  As,  Mn,  the  rare  earth  metals, 
and  actinides  from  aqueous  solutions  by  sorption  on  particles  of  aluminum  metal. 


e.  Clay  Materials 

Straub,  et  al  and  others  reported  the  use  of  various 

types  of  clays  (haolinite,  montmorillonlte  and  shales)  for  the  removal  of 
.■pecific  radioisotopes  from  watei'  and  waste  solutions.  The  clays  were  found 
to  be  especially  suitable  for  the  removal  of  (!s,  Zr  and  ITb.  With  the  addition 
of  5,000  ppm  of  clay,  removal  of  over  of  Cs  was  reported  Amphlett  and 

Sammon  reported  95%  removal  of  Sr.  Ttils  method  can  not  be  economically 

justified  because  large  vol'ames  of  clay  must  be  handled  both  initially  and  as 
a  contaminated  sludge.  This  method  has  its  greatest  potential  use  where 
geologic  and  hydrologic  conditions  are  such  that  the  ground  itself  may  be 
usea  for  disposal  of  the  sludge.  Studies  have  also  been  made  '  J 


use  of  clay  for  high  level  waste  disposal  by  fixing  specific  radioisotopes 
into  clay  by  means  of  high  temperat’ore  treatment. 
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f .  Oxidat ion»Reduct ion 

Kils  method  vas  especially  designed  to  inqprove  Ru  removal.  Ruthenium 

is  often  present  in  the  form  of  a  nitrosy  1-ruthenium  complex  which  is  difficult 

to  remove  by  coagulation  procedures.  An  approach  to  the  problem,  described 

(930 

by  Amphlett  and  Sammon  ,  Involves  pre-treatment  with  an  oxidizing  agent  to 
destroy  the  rathenium  complexes,  followed  by  a  calcium  phosphate-ferrous 
phosphate  floe  in  the  presence-  of  HaHSG^’  to  maintain  the  ruthenium  in  a 
reduced  form.  As  hi^  as  9Sii  of  ruthenim  was  removed  when  100  ppm  RaHSO^  was 
used. 

g.  Foaming 

A  foam  separation  technique,  based  on  surface  activity  phenomena,  was 

(lot) 

recent.ly  developed  by  Radiation  Applications,  Inc,,  Long  Island  City,  New  York'^ 

A  reduction  of  the  concentration  of  radioactive  strontim-90  and  cesium-137 

-<5  -12 

from  10  to  10  molar  in  low-level  nuclear  waste  was  claimed.  The  concen¬ 
tration  of  the  foaming  agent  is  one  of  the  most  critical  factors  in  detenninj.ng 
the  separation  factor  (surface  adsorption  factor).  A  detailed  discussion  of 
the  theoiy  as  well  as  experimental  observations  of  the  foam  separation  technique 

was  presented  bj-^  Schoen  and  Mazzella 
(tTci 

Lacy'*‘'^^  investigated  the  flotation  process  in  1957  and  reported 
r-'movals  of  up  to  'Jb’jh  for  Ba-l40-Ija-l40,  %jfo  for  Cs-137j  74^  for  Ru  and  68^  for 
Sr- 90,  These  studies  were  conducted  using  cetylpyridlnium  chloride  as  the  sur¬ 
face-active  agent  and  foaming  by  aeration. 

3.  Emergency  Methods  for  the  Decontami.nation  of  Radioactive  Water 

According  to  the  literature,  veiy  few  emergency  decontamination  units 
are  now  available  for  municipal  use,  although  the  U.S.  Amy  has  developed 
several  mobile  units  which  could  find  wide  application  for  public  use.  Several 
smaller  decontamination  units  are  commercially  available,  most  of  which  employ 
ion  exchange  techniques . 


-  j.Ou- 

Supplementary  treatment  to  the  conventional  rater  treatment  processes 
wo'jd-d  proha, bly  he  the  most  efficient  means  of  supplying  potahle  vater  to  a 
large  ntmiber  of  people,  folloving  radioactive  contamination  of  the  rater, 
a.  Mimicipal  Size  Decontamination  Unlt.e 

Woodrard  and  Roheck  report  on  an  ion  exchange  process  that  may  be 

used  to  supplement  the  normal  rater  treatment  proced'ore.  An  ion  exchange 
column  :in  the  form  of  a  cartridge  is  insert-ed  into  the  system,  and  is  disposed 
of  after  hr-eakthrougb  occurs,  ^fhe  type  of  resin  viH  vary  according  to  the 
time  after  detonation  that  the  rater  will  he  put  into  use,  A  mixed-hed  resin 
■'>"111  he  necessaiy  for  Immediate  use  of  the  rater,  vhereas  a  cation  resin  -will 
suffice  for  long^^term  usage,  due  to  the  decay  of  shcrt;-lived  iodine~131- 
ihe  life  of  each  Ion  cartridge  will  depend  -upon  the  total  solids  present, 
Altho'o^  large,  centrally  located,  decontamination  onits  ac'e  more 
economical,  smaller  decentralized  units  wo'^fLi  hs  preferable  as  the  problem 
of  ■'/■ater  distribution  vould  be  lessened, 
b •  Field  Decontamination  "nits 

Ihe  U.S.  Aimy  has  developed  several  mobile  decontamination 

units  primarily  for  use  by  troops  in  the  field,  bu-i;  these  co-^jld  also  be  used  to 
supply  rater  foi'  small  pop".ilation  groups  in  case  of  nuclear  attack.  One  such 
unit  is  comprised  of  a  floc-culator,  filters,  ■dual-bed  ion  exchange  column 
anc  a  chlorinator.  fhe  encire  unit  is  mounted  on  two  2-l/2  ton  true;  ,  A 
1,500  gph  output  is  obtainable.  Regeneration  of  the  ion  exchange  resins  is 
effected  by  hydrochloric  acid  and  soda  ash. 

e ,  Mlecellaneoas  recontam-infitlon  units 

hamotte  Uhemlcal  Company  has  developed  a  mixed  ion  exchange  resin 
bed  which  has  been  found  usefifL  for  the  decontamination  of  radioactive 


water.  The  -unit  consists  of  a  "bed  of  Amberlite  resin  MP-3  (lE-120  and  IRA-4lO) 
and  is  capable  of  producing  10  gallons  of  water  of  triple  distilled  quality. 

A  unit  developed  by  Lauderdale  and  Emmons  consists  of  twin 

columns  of  steel  vraol,  burnt  clay  and  activated  carbon  in  one  column  and  a 
mixed  bed  ion  exchange  resin  in  the  other-  Approximately  30  liters  of  water 
can  be  decontaminated  with  an  efficiency  ex.cee'i’ing  99-9^. 

Kease  and  Lacy  describe  a  method  of  radioactive  water  decon“ 

t  ami  nation  whereby  the  water  is  passed  through  a  column  of  natural  materials, 

.such  as  clay,  leaves,  humus,  gravel  and  sand.  Removals  of  over  90^  were 
reported.  The  problem  wo’uLd  be  to  obtain  uneontamlnated  materials  from  the 
environment  for  use  in  stich  a  column. 
k.  Diacusslon  and  Oonclusions 

ISie  more  important  radioisotopes  that  appear  initially  in  contaminated 
water  are  Ba,  Cs,  I,  La,  Ru  and  Sr.  Laboratory  data  indicate  that  the 
combination  of  conventional  water  treatment  methods,  namely,  chemical  coagu¬ 
lation,  lime-soda  ash  softening,  and  sand  filtration,  is  capable  of  removing 
99^  strontium.  The  removal  of  95?o  iodine  can  be  accomplished  by  the  addition 
of  silver  ions.  86^  cesium  can  be  removed  by  the  lime-soda  ash  softening  proced’ure 
when  proper  amounts  of  clay  are  added. 

The  efficiency  of  gross  activity  removal  lies  generally  between  50  and 

0*17' 

75^  ^.n  plant  scale  operation'  ' .  The  lower  plant  scale  efficiencies,  in 
comparison  to  the  labciratory  results,  are  believed  to  result  from  the  fact 
that:  (a)  only  one  or  two  decontamination  processes,  among  those  mentioned, 
are  employed  by  uhe  usual  water  treatment  plant,  Eind  the  operating  conditions 
may  not  be  optimized;  and  (b)  in  the  laboratory,  the  radioactive  materials  are 
generally  pre.sent  in  the  form  of  simple  salts.  In  contrast,  the  radioactive 
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materials  reaching  the  water  treatment  plant  are  generally  'belie-ved  to  he  In 

chemical  formn  which  are  more  difficult  to  remove  because  the  portion  which  is 

readily  removed  has  already  been  eliminated  by  natural  decontamination 

fLl8) 

mechanisms  before  reaching  the  water  plant.  Setter  and  Russell  ■  stated  that 
the  longer  the  activity  has  been  subjected  to  natural  purification,  the  more 
difficult  it  is  to  remove  the  remaining  soluble  activity  by  coagulation. 

It  should  be  pointed  out  here  that  the  decontamination  data  reported 
In  the  literature,  both  on  laboiatory  and  plant  scale,  are  based  on  one  of  the 
following  three  sotirces  of  radioactivity: 

(1)  added  radioactive  salts 

(2)  low-level  radioactivity  waste 

(3)  long-range  (world-wide)  fallout  from  atomic  bomb  tests 

The  chemical  and  physical  characteristics  of  radioactive  materials  from 
any  of  the  above  sources  is  different  from  that  present  in  local  fallout.  The 
latter  is  characterized  by  its  low  water  solubility  of  1^  to  2$  in  con^jarison 
to  the  30ia  or  higher  solubility  of  radioactive  materials  of  the  above  three 
categories.  Low  water  solubility  generally  will  Increase  the  removal  of  local 
fallout  by  conventional  water  treatment  methods  which  are  efficient  in  dealing 
with  particulate  or  colloid  materials. 

The  maximum  decontamination  factor  for  each  of  six  selected  elements 
above  are  listed  in  Table  X.  These  values  are  based  on  the  best  available  data 
in  the  literatiore.  It  is  en5)hasized  that  the  percentage  of  removal  is  greatly 
dependent  on  the  chemical  form  and  physical  state  of  the  radioactive  element  as 
well  as  the  concentration  of  the  treatment  additives  and  other  related  conditions 
such  as  pH  and  temperature  of  the  water.  As  an  example  of  the  decontamination 
factors  required  to  reduce  activity  to  the  MFC,  based  on  hypothetical 
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Table  X 

Maximum  Decontamination  Factors  Reported 


in  the  Literature  for  Selected 


opet; 


Decontami- 

naticn 


Decontamination  Factors  (D.F.)  and  References 


Method 

Ba-l40 

Gs-137 

1-131 

La-l40 

Eu-106 

Sr-89,90 

1.  Chemical 

3.^ 

1.5 

1.8 

3A 

12.5 

(84) 

2.0 

Coagulation 

(83) 

(80) 

(83) 

(8f) 

(80^ 

2 .  Lime-Soda 

10^ 

1.3 

NR 

10^ 

NR 

3-3  X  10' 

Softening 

(87) 

(88) 

(87) 

(80) 

3-  Rapid  Sand 

2  X  10^ 

2 

NR 

3.8 

NR 

1 .0- 

Filtration 

(80) 

(80) 

(80) 

^:8o) 

4.  Ion 

10^ 

10^ 

10^ 

10^ 

NR 

3.3  X  10 

Exchange 

(80) 

(80) 

(91) 

(80) 

(83) 

5.  Phosphate 

3.2 

1.5 

NR 

5  X  10^ 

6.7 

4,s  X  10' 

Coagulation 

(84) 

(97) 

(84) 

(84^ 

(94) 

6.  Clay 

NR 

5  X  10^ 

NR 

m 

2  X  10-" 

(80) 

(80) 

(91) 

7.  Metallic 

NR 

NR 

1.6 

NR 

m 

NR 

D’usts 

(80) 

8.  Coagulation  w/ 

NR 

NR 

2.2 

NR 

NR 

NR 

Silver  Ions 

(98) 

9.  Oxidation- 

NR 

NR 

NR 

NR 

10^ 

NR. 

Reduction 

(91) 

10.  Coagulation 

3.3  X  10^  3.3  X  10^ 

NR 

NR 

3-8 

NR 

with  Clay 

(99) 

(99) 

(99) 

11 .  Lime  Soda 

NR 

7-1 

NR 

NR 

4.2 

NR 

with  Clay 

(97) 

(97) 

12.  Foaming 

4.5 

6.7 

NR 

4.5 

3-9 

3.1 

(Flotation) 

(109) 

(109) 

(109) 

(109) 

(109) 

13 .  Slow  Sand 

NR 

NR 

2 

NR 

NR; 

NF 

Filtration 

(90) 

Notes . 

(1)  Percent 

Removal  e 

-  (sk 

.)ioo 

(2)  RN  -  not  reported  in  literature  reviewed 


contamination  of  the  Springfield,  Maes,  watershed  (as  reported  elsewhere  in 
this  report),  lahle  XI  was  prepared.  It  is  aijparent  from  this  table  that  ion 
exchange  is  the  only  single  decontamination  process  that  will  adeqitiatelv  remove 
the  minimum  concentration  reported  for  the  six  selected  isotopes  and  that  no 
process  will  adeq,iiately  remove  the  maximum  radiooontaminant ,  namely  I-I31. 
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